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This article presents an overview of current understanding of the interaction of low-energy 
positrons with molecules with emphasis on resonances, positron attachment and annihilation. 
Measurements of annihilation rates resolved as a function of positron energy reveal the presence 
of vibrational Feshbach resonances (VFR) for many polyatomic molecules. These resonances lead 
to strong enhancement of the annihilation rates. They also provide evidence that positrons bind to 
many molecular species. A quantitative theory of VFR-mediated attachment to small molecules is 
presented. It is tested successfully for selected molecules (e.g., methyl halides and methanol) where 
all modes couple to the positron continuum. Combination and overtone resonances are observed 
and their role is elucidated. Molecules that do not bind positrons, and hence do not exhibit such 
resonances, are discussed. In larger molecules, annihilation rates from VFR far exceed those expli- 
cable on the basis of single-mode resonances. These enhancements increase rapidly with the num- 
ber of vibrational degrees of freedom, approximately as the fourth power of the number of atoms in 
the molecule. While the details are as yet unclear, intramolecular vibrational energy redistribution 
(IVR) to states that do not couple directly to the positron continuum appears to be responsible for 
these enhanced annihilation rates. In connection with IVR, experimental evidence indicates that 
inelastic positron escape channels are relatively rare. Downshifts of the VFR from the vibrational 
mode energies, obtained by measuring annihilate rates as a function of incident positron energy, 
have provided binding energies for thirty species. Their dependence upon molecular parameters 
and their relationship to positron-atom and positron-molecule binding energy calculations are dis- 
cussed. Feshbach resonances and positron binding to molecules are compared with the analogous 
electron-molecule (negative ion) cases. The relationship of VFR-mediated annihilation to other 
phenomena such as Doppler-broadcning of the gamma-ray annihilation spectra, annihilation of 
thermalized positrons in gases, and annihilation-induced fragmentation of molecules is discussed. 
Possible areas for future theoretical and experimental investigation are also discussed. 
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I. INTRODUCTION AND OVERVIEW 

The subject of this review is the interaction of low- 
energy positrons with molecules. Positrons, the antipar- 
ticles of electrons, are important in many areas of science 
and technology. Much of their utility relies on the fact 
that, when an electron and positron interact, they can 
annihilate, producing a characteristic burst of gamma- 
rays. The lowest order process results in two back-to- 
back photons, each with the energy of the rest mass of 
the electron (or positron), 511 keV. 

The annihilation of low-energy (e.g., < 50 eV) 
positrons on atoms and molecules plays a particularly 
important role in many fields. In medicine, positron 
emission tomography (PET) exploits two- g amma annih i- 
lation to study human metabolic processes (jWahll . |2002|) . 
In materials science, there are numerous posit ron-based 
techniques to study the propertie s of matter (jColeman , 
2000 1 : i Dupasauier and Millsl Il995l iPuska and Nieminenl . 
1994 ISchultz and Lvnnl. Il988l ) ; inclu ding the Fermi sur- 

microscopic pores 
the free volume in 



faces in metals (IMaior et al. 



in solids (Gi dlev et al 




polymers (|Dlubek et all 19981). and the c omposition and 
structure of surfaces ( David et all l200l[ ) . In astronomy, 
511 keV annihilation radiation, the strongest gamma ray 
line of extraterrestrial origin, ha s proven useful i n elu- 
cidating astroph y sical p rocesses (IChurazov et all 120051: 



iGuessoum et all l2010t iRamatv et al , 19921) . A cur 



rent research goal is the creation of a Bose condensate 
of positronium (Ps) atoms (i.e., the electron-positron 
analog of the hydrogen atom) that offers promise for 
the development of an annihilation gamma-ray laser 



(ICassidv and Millsl . 120071: iMillsl . 120021 120071: iMills et all 
120041) . 

Typically, positrons from conventional sources (e.g., 
radioisotopes or electron accelerators) slow down from 
energies of kilovolts to hundreds of kilovolts to < 50 eV 
before annihilating. In the case of atoms or molecules, 
if the incident positron energy e is greater than the Ps- 
formation threshold E t h = Ei — Ep s , where Ei is the 
ionization energy of the target and Ep s — 6.8 eV is the 
binding energy of the ground-state Ps atom, then the 
dominant annihilation process is through Ps formation. 
The resulting Ps atom subsequently annihilates by emit- 
ting two or three gamma-ray quanta. 

In this review, attention is restricted to positron en- 
ergies below the Ps-formation threshold, < e < E t h, 
where the Ps channel is closed. Here annihilation occurs 
as a result of the overlap of the positron and electron 
densities during the collision. The basic rate in this case 
is the Dirac r ate A n for tw o- gamma annihilation in a free 
electron gas ( Dirad . Il930[ ) 



Ad = 7rrocn e , (1) 

where ro is the classical electron radius, tq — e 2 /mc 2 in 
cgs units, e and m are the electron charge and mass, c is 
the speed of light, and n e is the electron density. 

In his semin al discovery of the positronium atom, 
iDeutschl (Il951allbh found a curious effect. Although the 
annihilation rate for thermal positrons at 300 K in atomic 
and molecular gases was approximately in accord with 
Eq. (HJ) for some species (e.g., argon and nitrogen), the 
rate for dichlorofluoromethane CCI2F2 ("Freon-12") was 
much larger. Deutsch insightfully ascribed this effect 
to some type of resonant positron-molecule atta c hment 
process. A decade later, iPaul and Saint-Pierr 1 (119631) 
measured annihilation rates in gases of alkane molecules 
C n H2„+2, from methane to butane, n = 1-4. They found 
that the rate A was much greater than Ad and that the 
ratio A/Ad increased exponentially with molecular size. 

Annihilation rates in gases are conventionally normal- 
ized to the Dirac rate. The corresponding dimensionless 
quantity is the "effective number of electrons" 1 



•Zeff = 



A 



(2) 



where n is the density o f atoms or molecules ( Fraserl . 
[196111 

omeranchukl, [l949) . For a simple collision, and ne- 
glecting electron-positron correlations, one might expect 
that A ~ Ad, so that Z c g is comparable to Z = n e /n, the 
total number of electrons per target atom or molecule. 
However, values of Z e f{ are often much larger (e.g., for 
butane, Z eS /Z = 600). 



1 In chemical kinetics, Z e g corresponds to the (normalized) rate 
constant of the annihilation reaction. In positron physics this 
quantity is commonly referred to as the "annihilation rate" . 
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Positron annihilation in atoms and 
was subs equently studied for a wide 
speci es 



1980, 2002, 



(lAl-Qaradawi et all l2000t 



McNutt et al. 




Hevland etaR J1986 



molecules 
range o f 
ICharlton et al 



19851 : ISharma and McNuttl . Tl978t iTaol. 



1985 



Osmon ~ Il965allb1: ISharma et al. 



1965 



19701) 



Early experiments were done with thermal po sitrons in 
gases at atmospheric densi t ies, n ~ 1 amagat (iDeutscb , 
1953 [lGriffith and Hevland . 1 19781 IPaul and Saint-Piem . 
1963h . 2 Later, experiments were done at much lower 



densities usi n g posi t rons tr a pped and cooled to 300 K 
(llwata et all Il995l Il997al IMurphv and Surkol . Il99ll 
ISurko et aZj . ll988bl) . The experiments showed that the 
annihilation rates for many molecular species exceeded 
greatly the naive benchmark rate, Z e g ~ Z, and a 
number of chemical trends were identified. 

Since Deutsch's first results, these large anni- 
hilation rates were associated with some kind of 
resonance phenomeno n or attachment process. 



iGoldanskii and Savasovl ( 19641 ) discussed the possi- 
bility of resonance-enhanced annihilation due to a 
bound or virtu a l posi tron state close to zero energy. 
ISmith and Paull ( 1970h considered the possibility that 
the large annihilation rates in molecules were due to 
a vibrational resonanc e, and sever a l othe r explana- 
tions were proposed (iDzuba et al , 1996 ; Gribakinl . 



120001 : lLaricchia and Wilkinl . Il997t ISurko et a,l\ . Il988bh . 

However, progress was hampered greatly by the lack 
of data other than for positrons with thermal energy 
distributions at 300 K. The summary statement in 
1982 by Sir Harrie Massey was that annihilation studies 
were "compl etely myster i ous at present in a lmost all 
substances" ( Fraser et al l ll982UMassevLll982h . and this 
remained more or less correct for another twenty years. 

In the broader view, processes that are common- 
place in physics involving matter, such as low-energy, 
two-body scattering events, have frequently been found 
to be frustra tingly difficult to study when antiparticles 
are in v olve d (Charl ton and Humberstonl.l200ll:IColeman| . 
2000 : lEades and Hartmannl . Il999l ISchultz and Lynnl . 
19881 ) . The advent of efficient positron traps 



marked a t urning poin t ( Murphy and Surkol . Il992t 



ISurko et al 



of studies 



2005 1 Il988ah, e n abling a new generation 
Iwata et all, Il995l Il997at iKurz et all 1 19961: 



IMurphv and Surkol Il99lt ISurko et all Il988bl ). Experi- 
ments with trapped positrons cooled to 300 K permitted 
studies of test species at low densities (e.g., < 10~ 6 am- 
agat). This ensured that annihilation was strictly due 
to binary collisions, rather than many-particle effects 
(jlakubov and Khrapakl . Il982f ). and it enabled study of 
a broader range of chemical species, including low- 
vapor-pressure liquids and solids. Gamma-ray spec- 



tra w ere measured for many molecules ( Iwata et all , 
1997a). The Z e g jZ ratios were found to increase rapidly 
with molecular size up to species as large as naphtha- 
lene and hexadecane (C1 6 H34), reaching valu es > 10 4 
(IMurphv and Surkol Il99ll: ISurko et all . Il988bl) . 3 

A key to further progress was the development of a 
trap-based p ositron beam with a narrow energy s pread 
(~ 40 meV) ([Gilbert et aLl . ll997llKurz et alJ . ll998h . Us- 
ing this beam, annihilation rates for atoms and molecules 
were measured as a function of incident positron energy 
from 50 meV to many electron volts. The result was the 
discovery of resonances associated with the molecular vi- 
bration al modes, namely vib rational Feshbach resonances 
(VFR) (ICilbert rl aLl . l2002h . 

A crucial point is that VFRs generally require the ex- 
istence of a bound state of the positron and the molecule. 
They occur when the incident positron excites a vibra- 
tional mode and simultaneously makes a transition from 
the continuum into the bound state. The existence of 
both low-lying vibrational excitations and a positron 
bound state thus enables the formation of long-lived 
positron-molecule resonant complexes in a two-body col- 
lision. The lifetime of these quasibound states is limited 
by positron autodetachment accompanied by vibrational 
de-excitation. The upper limit on the lifetime is < 0.1 ns, 
set by the positron annihilation rate in the presence of 
atomic-density electrons. 

The annihilation rate as a function of positron energy, 
Z e g(e) (i.e., the "annihilation spectrum" ) for the four- 
carbo n alkane, butane, is shown in Fig. [T] ( Gi lbert et all . 
I2002D . While there is some qualitative correspondence 
between the Z c f{(e) and the infrared (IR) absorption 
spectrum of the molecul e, the shapes of the spectral fea- 
tures are quite different ( Barnes et all 120031 ) . 

These positron VFRs can be compared to resonances 
that play an important r ole in electron attachment 
to molecules and clusters ( Christophorou et all . Il984l 
iHotop et al .], 120031) . The electron collision results in 
the production of long-lived (metastable) parent an- 
ions, or molecular fragment negative ions via dissocia- 
tive attachment. A dominant mechanism of electron 
capture by molecules i s via negative-ion resonant states 
( Christophorou et all [l984). Dissociative attachment 
usually proceeds via electron shape resonances of ground 
or electronically-excited molecules. Such resonances are 
quite common in diatomic, triatomic and polyatomic 
species at energies in the range ~ 0-4 eV. The the- 
oretical description of them involves (complex) Born- 
Oppenheimer potential-energy s urfaces ( Bardslevl [i968a; 



_po_ 

iDomckd . Il98lt lO'Malleyl Il966l ). All the data indicate 
that positrons generally do not form shape resonances or 
electronic Feshbach resonances in low-energy collisions 



2 1 amagat = 2.69 X 10 19 cm -3 is the density of an ideal gas 
at standard temperature and pressure, 273.15 K and 101.3 kPa, 
respectively. 



3 The theoretical maximum for the magnitude of Z c ff is 
given by the unitar i ty li mit of the inelastic cross section 
JLandau and Lifshitd. 1 19771 ). Z eff < 10 7 for room-temperature 
positrons. 
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FIG. 1 The normalized annihilation rate Z e g(e) for butane 
C4H10 (•) as a function of the total incident positron energy 
e: (a) up to the Ps formation threshold, E t h = 3.8 eV, and (b) 
in the region of the molecular vibrations; dotted cur v e, the 
infrared absorption spectrum (|Linstrom and Mal lard. 2005) 
(logarithmic vertical scale, arbitrary units); solid curve, the 
vibrational mode density (in arbitrary units) , with the modes 
represented by Lorentzians with an arbitrary FWHM of 
10 meV; dashed line, mean energy of the C-H stretch fun- 
damentals. 

with molecules. Instead, energy-resolved annihilation 
studies point to the important role played by the VFR. 

These vibrational (or "nuclear-excited" ) Feshbach res- 
onances involve coupling of the electronic and the nu- 
clear motion beyond the Born-Oppenheimer approxima- 
tion. It cannot be described by potential-energy surfaces. 
This type of res onances was originally introduced by 
(|Bardslevl Il968bl) as an "indirect" mechanism for disso- 
ciative electron recombination and described using Brcit- 
Wigner theory. In the case of electrons, these VFRs 
lead to large attachment cross sections which typically 
reach their maximum values at thermal electron energies 
( Christophorou et all |l984). They are also responsible 
for the formation of long-lived parent negative ions for 
many complex polyatomic molecules. 

Referring to Fig. [1] the energy of the VFR correspond- 
ing to mode v is given by energy conservation, 

s v =u v - e b , (3) 

where is the positron-molecule binding energy, and 
uj v is the vibrational mode energy. The positron bind- 
ing energy (i.e., the positron affinity) can be measured 
by the downshift of the resonances from the energies of 
the vibrational modes. In Fig. Q] for butane, this is 
most easily seen as the shift in the C-H stretch vibra- 
tional resonance. The corresponding peak in Z e g occurs 



at 330 meV, as compared with the vibrational mode fre- 
quency of 365 meV, indicating that et = 35 meV. The 
resonances at lower energies are due to C-C modes and 
C-H bend modes and exhibit the same downshift. 

There are a number of important chemical 
trends asso ciated with reson ant annihilation on 
mole c ules (llwata et all, Il995t iMurphv and Surkol . 
Il99ll ; lYoung and Surkol l2008bl lcj). Examples are shown 
in Table HI Very small molecules, such as CO2, CH4, 
or H 2 0, have relatively small values of Z c g (e.g., 
Z c ff/Z < 10), and typically they do not exhibit resonant 
annihilation peaks. Positrons either do not bind to these 
species (i.e., < 0), or they bind extremely weakly. 
With the exception of methane, all of the alkanes exhibit 
VFRs, with values of e b increasing linearly with the 
number of carbon atoms n, and the magnitudes of Z c g 
increasing approximately exponentially with n. Most 
hydrocarbons, including aromatic molecules, alkenes and 
alcohols, exhibit similar resonant annihilation spectra. 



TABLE I Annihilation rates Z c g, and binding energies et for 
selected molecules. 



Class 


Molecule 


Z 


e b (meV)° 


ry a 

^cff 


Small inorganics 


H 2 


10 


< 


170 b 




NH 3 


9 


> 


300 b 


Methyl halides 


CH 3 F 


18 


> 


250 6 




CH 3 Br 


44 


40 


2000 6 


Alkanes 


CH4 


10 


< 


70 b 




C2H6 


18 


> 


900 c 




C 3 H 8 


26 


10 


10 500 c 




C6H14 


50 


80 


184 000 c 




C12H26 


98 


220 


9 800 000 c 


Alcohols 


CH3OH 


18 


> 


750' 




C2H5OH 


26 


45 


4500 6 


Aromatics 


CsH6 


42 


150 


47 000 c 




CioHg 


68 


300 


1 240 000 c 



"Valu es from energy-resolved measurements |Young and Surkol . 
2008b c); typical uncertainties in Z e ff and £5 are ±20% and 
±10 meV, respectively. 

'Maximum values for positron energies e > 50 meV. 

c Values of Z e g at the C-H resonance peak. 



Much progress has been made in the theoretical undcr- 
standing of resonant positron annihilation in molecules 
(lGribakkl2rooll200ll:lGribakin and Gilil2004D . A quan- 
titative theory has been developed for the case of isolated 
resonances of IR-active vibrational modes, such as those 
observed in experiments for selected small molecules. 
The prototypical example is that of the methyl halides, 
CH3X, where X is a F, CI or Br atom. Positron cou- 
pling to the IR-active modes is evaluated in the dipole 
approximation using data from IR absorption measure- 
ments. The only free parameter in the theory is the 
positron binding energy, which can be taken from ex- 
periment. This yields theoretical annihilation spectra for 
methyl halides that are in good agreement with the mea- 
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surements ( Gribakin and Led l2006bl ). 

A more stringent test of the theory relies on the fact 
that positron binding energies are expected to change 
little upon isotope substitution. For deuteration this 
was confirmed experimentally. The binding energies 
measured for CH3CI and CH^Br were used to predict 
Z e s for their deuterated analogs. The result is excel- 
lent agreement between theory and experim ent with no 
adjustable parameters (jYoung et all |2008() . In other 
small molecules, such as ethylene and methanol, IR- 
inactive modes and multimode vibrations are promi- 
nent, and must be included to explain the observations 
(|Young et aLl . l2008l ). 

This theoretical approach explains Z c s for small poly- 
atomics in which the positron coupling to the mode-based 
VFR, and possibly a few overtones, can be estimated 
(e.g., when they have dipole coupling). Their Z e g val- 
ues are between a few hundred and a few thousand. 4 
However, larger molecules with more than one or two 
carbons, have values of Z c g that cannot be explained 
by this theory (cf. Fig. [T] for butane). The current 
physical picture ascribes their large annihilation rates 
to large densi ties of vibrational res onances, known as 
"dark states" (|Gribakinl . l2000l l200ll) . that are not cou- 
pled directly to the positron continuum. The positron 
first attaches to the molecule via a vibrational "door- 
way state" (e.g., a dip ole-allowed mode- based VFR) 
(jGribakin and Gilll . l2004l ). The vibrational energy is then 
transferred to the "dark states" in a process known as 
intramolecular vibrational energy redistribution (IVR). 
Such IVR is important for many physical and chemical 
processes in molecules, including dissociative att achment 
(jNesbitt and Field! . Il996l lUzer and Milled . Il99lh . 

The magnitudes of resonant contributions to Z e g ex- 
hibit a relatively weak dependence on and on the in- 
cident positron energy e. It is of the form g = \J et/s 
and follows from rather general theoretical considera- 
tions. When this dependence is factored out, it is found 
experimentally that the resulting quantity, Z e fi/g, scales 
as ~ iV 4 , where N is the number of atoms in the molecule. 
This dependence on N is thought to reflect the rapid in- 
crease in the density of the molecular vibrational spec- 
trum with the number of vibrational modes. This depen- 
dence is interpreted as evidence that IVR does indeed 
play an important role in the annihilation process. 

Estimates of Z e g in large molecules, which assume that 
the IVR process is complete and all modes are populated 
statistically, predict Z e s values far in excess of those that 
are observed. Such estimates also fail to reproduce the 
energy dependence of Z c s, which is largely determined by 
the mode-based vibrational doorways. One hypothesis, 
as yet unconfirmed, is that the IVR process does not 



run to "completion". It appears that selective coupling 
of multimode vibrations leaves a large portion of them 
inactive. The calculation of Z e ff then requires a detailed 
knowledge of the vibrational mode couplings, and this 
has not yet been done. 

The energy-resolved annihilation experiments provide 
measures of positron-molecule binding energies, either di- 
rectly using Eq. (|3]), or for very weakly bound states, 
indirectly through the dependence of Z c s on g. To 
date, binding energies for about thirty molecules have 
been measured. They range from ~ 1 meV in small 
molecules such as CH3F, t o ~ 300 meV for large alkanes 
( Young and Surkd l2008bl ld\ A recent analysis indicates 
that these binding energies increase approximately lin- 
early with the molecular dipole polarizability and dipole 
momen t, and for aromatic molecules, the number of n 
bonds ( Danielson et al .1. 120091) . 



For atoms, comparatively accurate positron binding 
energies have been predicted theoretically for about ten 
species (|Mitrov et all l2002h . but there are no mea- 
surements. There have been a number of calcula- 



tions for positron binding to molecules ( 


Bressanini 


et al. 


1998: iBuenker and Liebermann. 2008 


: Bucnker et al. 


2005: lOarev et all I200S 


Choinacki and Strasburgcr 


2006: Danbv and Tennvsonl, 


1988: Gianturco et al. 


, 2006 


Kurtz and Jordanl 1978, 19811; Schrader and Moxom 


20011 ISchrader and Wanel.ll97fiUStrasburgeil 19991 2004 



large dipole moments which facilitate binding. In con- 
trast, the molecules for which the binding energies are 
known from experiment are either nonpolar or only 
weakly polar. Thus at present, there are no species for 
which experiment and theory can be compared, and so 
this is a critical area for future research. 

Presented here is a review of theoretical and experi- 
mental results for positron annihilation on molecules in 
the range of energies below the positronium formation 
threshold. Emphasis is placed upon the case in which 
positrons bind to the target and annihilation proceeds 
via the formation of vibrational Feshbach resonances. 
Current knowledge of positron- molecule binding energies, 
obtained from both experiment and theoretical calcula- 
tions, is summarized. These results are related to studies 
of positron-induced fragmentation of molecules, annihi- 
lation gamma ray spectra, annihilation in dense gases 
where nonlinear effects are observed, and to analogous 
electron interactions with molecules and clusters. 



II. THEORY 

A. Annihilation basics 



4 The heights of resonant peaks in the measured Z c ff spectra are 
related to the energy spread 8e of the incident positron beam. 
The value of Z c g ~ 10 3 corresponds to the typical 8e ~ 40 meV 
used to date (see Sec. IIIII . 



The process of electron-positron annihilation is de- 
scribed by quantum electrodynamics (QED). In the non- 
relativistic Born approximation, the cross section for an- 
nihilation into two photons averaged over the electron 
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and positron spins is ( Berestetskii et a/lll982f) 



= 7IT ~, 
V 



(4) 



where u is the relative velocity of the two particles. This 
cross section obeys the 1/v threshold law which describes 
inelastic collisions with fas t particles in the final state 
( Landau and LifshitzJ . [l977h . 

The two-photon annihilation described by Eq. (j4|) is al- 
lowed only when the total spin S of the electron-positron 
pair is zero. For S = 1 the smallest possible number of 
annihilation photons is t hree. The corresponding spin- 
averaged cross section is ( Berestetskii et all , 11982 



CT 3 7 = x(tt 2 ~ 9)ar 2 -, 

O D 



(5) 



where a = e 2 /he (in cgs units) is the fine structure con- 
stant, a sa 1/137. Since CT3 7 is 400 times smaller than 
(727, positron annihilation in many-electron systems is 
dominated by the two-gamma process. 

Numerically, the cross section in Eq. ((4]) is CT27 ~ 
10~ 8 c/t> a.u. 5 Hence the annihilation rate is usually 
much smaller than the rates for other atomic collision 
processes, even at low positron velocities (e.g., thermal, 
v ~ 0.05 a.u. at 300 K). When a fast positron, such 
as that emitted in a /3 + decay, moves through matter, it 
loses energy quickly through collisions, first by direct ion- 
ization, positronium formation and electronic excitation, 
and then by vibrational excitation and elastic collisions. 
As a result, the positrons typically slow to thermal ener- 
gies (i.e., ~ 25 meV for T = 300 K) before annihilation. 

At small velocities, v < 1 a.u., Eq. (Q| must be mod- 
ified to take into account the Coulomb interaction be- 
tween the electron and positron. The typical momenta 
exchanged in the annihilation process are p ~ mc. The 
corresponding separation, r ~ h/mc, is small compared 
to do, and in the nonrelativistic limit, the annihilation 
takes place when the electron and positron are at the 
same point. The cross section in Eq. (|4| must then 
be multiplied by the prob ability density at the origin 
(jLandau and Lifshitzl . Il977l ). 



|V(0)| 2 = 



2tt 



-2ir/v\ ' 



(6) 



where the wave function ip is normalized by — e ° 
at r 3> &0- This increases the annihilation cross section. 

The annihilation cross section for many- el ectron 
targets is t r aditio nally written as ( Fraseii Il968t 
|Pomeranchukl,ll949l) 



where Z e g represents the effective number of electrons 
that contribute to the annihilation. In the Born approx- 
imation, Z c g — Z, the total number of target electrons. 

However, at small positron energies (e.g., e < 1 eV), 
Z c g can be different from Z. First, there is a strong 
repulsion between the positron and the atomic nuclei. 
This prevents the positron from penetrating deep into 
the atoms, so that the annihilation involves predomi- 
nantly electrons in the valence and near-valence sub- 
shells, thereby reducing Z e g. On the other hand, the 
positron is attracted to the target by a long-range polar- 
ization potential — a,j/2r 4 , where ay is the target dipole 
polarizability, which enhances Z Q g . There is also a short- 
range enhancement of Z e g due to the Coulomb inter- 
action between the annihilating electron and positron, 
which has the same origin as the expression in Eq. ([5]). 
Finally, if the target binds the positron, the annihi- 
lation can be enhanced by positron capture into this 
bound state. The cross section for radiative capture (i.e., 
by emission of a photon ) is small, namely a c ~ a 3 a 2 , 
(Ber estetskii et al I Il982h . In collisions with molecules, 
the positron can transfer its energy to vibrations, form- 
ing a positron-molecule complex. This process is effective 
in enhancing the annihilation rate. It is the principal fo- 
cus of the present review. 

As follows from its definition by Eq. |[7J, Z e s is equal 
to the electron density at the positron, 



ieff 



r z 



r - rj)|* k (ri, . . . ,r z ,r)| 2 dri . . . dr z dr 



(8) 

where ^(fi, ... ,tz,t) is the total wave function of the 
system, with electron coordinates and positron coor- 
dinate r. This wave function describes the scattering of 
the positron with initial momentum k by the atomic or 
molecular target, and is normalized to the positron plane 
wave at large separations, 



*k(ri, 



,rz,r) 



$o(ri, 



k-r 



(9) 



where $0 is the wave function of the initial (e.g., ground) 
state of the target. For molecules, both and $0 also 
depend on the nuclear coordinates, which must be inte- 
grated over in Eq. ([5]). 

Equations and © determine the annihilation rate 
in binary positron-molecule collisions, 



A 



a a vn 



nrlcZ cS n, 



(10) 



c a = o"2 7 Z cff = 7rrQ-Z c 



(7) 



We make use of atomic units (a.u.), in which m = |e| = % = 1, 
c = a -1 Ri 137 a.u., and the Bohr radius ao = h 2 /me 2 (in cgs 
units) also equals unity. 



where n is the gas number density. To compare with 
experiment, this rate is averaged over the positron energy 
distribution. For thermal positrons this distribution is a 
Maxwellian, while in beam experiments, it is determined 
by the parameters of the beam. Empirically Eq. (|10p is 
also used to describe experiments at high densities where 
Z c g becomes density dependent (see Sec. IVIII.Cp . 
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B. Gamma-ray spectra and annihilation rates 

In the nonrelativistic limit, the two-photon QED an- 
nihilation amplitude can be expressed in terms of an ef- 
fective annihilation operator 



o (P) = / e-^>(r)^(r)dr, 



(11) 



where ip(r) and ip(r) are the electron and positron de- 
struction ope rators, 6 and P is the total momentum of 
the photo ns (|Dunlop and Gribakinl . l2006t iFerrelll . Il956t 
Led . Il957| ). The probability distribution of P in an anni- 
hilation event is given by 



W)(P) 



(f\O a (P)\i) 



(12) 



where |i) is the initial state with Z electrons and the 
positron (e.g., that with the wave function ^k), and |/) 
is the state of Z — 1 electrons after the annihilation. 
For P = 0, the two photons are emitted in opposite 



shifted by e n /2 relative to mc 2 . This shift is much smaller 
than the typical Doppler shift e due to the momentum 
of the bound electron, P ~ y / 2m\e n \, which corresponds 
to e ~ Pc ~ ^J\e n \mc 2 ^> |e n |. The resulting width and 
shape of the gamma spectrum contain important infor- 
mation about the bound electrons. 

In most experiments, the annihilation photons are not 
detected in coincidence with the final state /, and the ob- 
served spectrum is the sum over all final states, wie) = 
^2fWf(e). However, this spectrum still reveals contri- 
butions of different final states. For example, in par- 
tially fluorinated alkanes, annihilation with the tightly 
bound fluorine 2p electrons results in a broader spectral 
component than annihilation with the more diffuse C- 
H bond electrons. This allows one to deduce the rel- 
ative fraction of th e corresponding annihilation events 
(llwata et all ll997bT) (see Sec. lVllLAD . 

The total annihilation rate in the state i leading to the 
final state / is obtained by integration over the momenta, 



directions and have equal energies, E^i = E. 
mc 2 
e.g., 



72 



E~ 



A / =7rrgc/K/|(5 B (P)|i)| ! 



d 3 P 
(2tt) 3 ' 

For P ^ the photon energy is Doppler-shifted, and ^ toM annihilation rate in statc { is 



(16) 



E 7 i = E 7 + mc\V\ cos8, 



(13) 



f 



(i\n-(r)h+(r)\i)dr, (17) 



where V = P/2m is the center-of-mass velocity of the 
electron-positron pair, and 9 is the angle between V and 
the direction of the photon. Averaging the distribution 
of the Doppler shifts e = E 7 i — _E 7 = (Pc/2) cos 6*, over 
the direction of emission of the photons, gives the photon 
energy spectrum, 



w f (e) = i 



2|e|/c 



W f (P 



PdPdrip 

(2tt) 3 ' 



(14) 



where n_(r) = ^ (r)i/j(r) and n+(r) = (p'(r)ip(r) are the 
electron and positron density operators. The annihila- 
tion rate is thus given by the expectation value of the 
electron density at the positron. Equation (IT71) gives the 
two-photon annihilation rate in a system of one positron 
and one target atom or molecule. For a positron moving 
through a gas of density n, the annihilation rate takes the 
form of Eq. (ITU1) . Normalizing the initial state i to one 
positron per unit volume, as \l/k in Eq. (0), one obtains 



In Cartesian coordinates, 



Wf (P x ,P y ,2e/c)^^. (15) 



(i|n_(r)n+(r)|i)dr. 



(18) 



This form shows that the energy spectrum is proportional 
to the probability density for a component of P. This 
quantity can be measured either by sampling the Doppler 
spectrum of the gamma rays or by measuring the angular 
deviation of the two photons (see Sec. IIII.FI) . 

When a low-energy positron annihilates with a bound 
electron with energy e n , the mean photon energy E 1 is 



In Eq. the spin indices in V( r ) an d V?( r ) are suppressed, 

and summation over them is assumed. This form can be used 
in systems with paired electron spins or when averaging over the 
positron spin. The modulus-squared amplitude is then multiplied 
by the spin-averaged QED factor nr^c. In general, one should 
use the spin-singlet combination of the annihilation operators 
in Eq. Hilt . — $4.$t)i together with the two-photon 

annihilation factor AttTqC. 



In the coordinate representation, this yields Eq. ©. 

In the independent-particle approximation, the elec- 
tronic parts of the initial and final states are Slater de- 
terminants constructed from the electron orbitals (e.g., 
in the Hartree-Fock scheme). The incident positron is 
described by its own wave function (/?k(r), and the anni- 
hilation amplitude (f\O a (P)\i) takes the form 

A„ k (P) = J e - lP - r ^(rVk(r)dr, (19) 

where ift n ( r ) is the orbital of the annihilated electron. In 
this approximation 



z 

Z cS = J2 / lV>n(r)| 2 |^ k (r)| 2 dr, 

n=l 



(20) 



i.e., the average product of the electron and positron den- 
sities. 
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C. Positron-molecule wave function 



where p ep is the average electron density at the positron, 



The wave function for t he positron coll i ding with 
a molecule, can be written as ()Gribakinl . 120001 l200lh 



(o) 



(0)\ 



V 



e — e 



V 



(21) 



The first term on the right-hand side describes direct, or 
potential scattering of the positron by the target. The 
corresponding wave function ^H^p is determined by the 
positron interaction with the charge distribution of the 
ground-state target and electron-positron correlation ef- 
fects (e.g., target polarization and virtual Ps formation). 
It neglects the coupling V between the electron-positron 
and nuclear (vibrational) degrees of freedom. 

The second term describes positron capture into the vi- 
brational Feshbach resonances. It is present for molecules 
that can bind the positron. These resonances corre- 
spond to vibrationally-excited states ^ v of the positron- 
molecule complex, embedded in the positron continuum. 
They occur when the positron energy e = k 2 /2 is close to 
e„ = E„ — Eb, where is the positron binding energy, and 
E v is the vibrational excitation energy of the positron- 
molecule complex. Equation (|2ip has the appearance of 
a standard perturbation-theory formula, but the energies 
of the positron-molecule quasibound states >]/„ in the de- 
nominator are complex, e v — |r„, where T v = r^+r^+rj, 
is the total width of the resonance. In atomic units T v is 
equal to the decay rate of the resonant state. It contains 
contributions of positron annihilation and elastic escape, 
r° and r®, respectively, and possibly also the inelastic 
escape rate T % v . The latter describes positron autode- 
tachment accompanied by vibrational transitions to the 
states other than the initial state. 

Molecular rotations are, in general, not expected to af- 
fect positron annihilation. The rotational motion is slow 
compared to the motion of the positron or the vibrational 
motion. Accordingly, direct scattering can be considered 
for fixed molecular orientation, and the results averaged 
over the orientations. Positron capture in VFRs at low 
energies is dominated by the s wave, or at most a few 
lower partial waves. Hence in the capture process, the 
angular momentum of the molecule remains unchanged 
or changes little. 

The positron capture amplitude (^iVH^k '') deter- 
mines the elastic width in state v. 



rt = 2TT K*„|V|* 



(OK 



(2^)3 



(22) 



If the positron interaction with the vibrations cannot be 
described by perturbation theory, Eqs. (l2Tj) and (|22|) re- 
main valid provided the amplitudes (^^IVI^^) are re- 
placed by their nonperturbative values. 

According to Eq. (fTTj) . the annihilation rate of the 
positron-molecule state ^> v is given by 



r" 



r C Pep 



(23) 



Pep 



f Z 

J 226(T - Ti)\^ u (Ti, . . . ,T Z ,r)\ 2 dri . . .dTi 



dr. 



(24) 

with the integration extending to the nuclear coordinates 
in the wave function \E> '„. The amplitude of the nuclear 
motion is small, and p ep is expected to depend weakly on 
the degree of vibrational excitation in state v. 

To calculate Z c g, the wave function from Eq. (|2"Tj) is 
substituted into Eq. ©, which yields 



Z cS = (* k |^<5(r-r,)|^k) 



i=l 
Z 



= (4 0) iE^-^i4 H{ int !*r ce ) 

i=l I k J 



(e-e„-£r M )(e-e,, + ir i ,) 



. (25) 



where the capture amplitude A v is related to the elastic 
width by Y e v = 2tt\A v \ 2 . The terms on the right-hand 
side describe the contributions of direct and resonant an- 
nihilation and the interference between the two. We will 
now examine the two main contributions in detail. 

The separation of the wave function into the direct and 
resonant parts in Eq. (|2"Tj) is valid because the positron 
VFRs are narrow. This is a consequence of the weakness 
of coupling between the positron and the vibrational mo- 
tion (i.e., small capture widths T^, see Sec. lII.F| ). In spite 
of this, the resonant contribution to the annihilation rate 
for complex polyatomics exceeds the direct contributions 
by orders of magnitude. 



D. Direct annihilation: virtual and weakly bound positron 
states 

The potential scattering wave function satisfies 



the Schrodinger equation 

(T+U-E p 



(0) 
k 



(0) 



(26) 



where T is the kinetic energy operator for the electrons 
and positron, U is the sum of all Coulomb interactions 
between the particles (with the nuclei at their equilibrium 
positions), and E is the target ground-state energy. 

For positron energies below the Ps-formation thresh- 
old, annihilation occurs when the positron is within the 
range of the target ground-state electron cloud. At such 
distances, the interaction U between the particles is much 
greater than the positron energy e. Therefore, the e^ k ° 

term in Eq. (PIBl) can be neglected, and the solution 

at these small separations does not depend on e, except 

through a normalization factor. 
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When the positron is outside the target, con- 
tains contributions of the incident and scattered positron 
waves, 



,r z ,r) ~ $ (ri, ...,r z ) 



Akr 



+ /kk- 



where /kk' is the scattering amplitude, and k' = kr/r. 
Inside the target, 4^ is determined by matching it with 
Eq. (f27l) at the target boundary r = R, where R is the 
characteristic radius of the target. For small positron 
momenta, kR -C 1, the scattering is dominated by the s 
wave, and the amplitude /kk' can be replaced by the s- 
wave amplitude f n . As a result, the integrand in Eq. (Op 
for Z e g is proportional to |l+/o/i?| 2 (|Dzuba et aU [l993). 
This gives t he following esti mate for Z e g due to direct 



this gives the tollowmg esti m 
annihilation (|Gribakinl . 12000). 



* 4np e 8R (R 2 



2i?Re/ + ^l 

47T 



(28) 



where p e is the effective electron density in the region 
of annihilation, SR is the range of distances where the 
positron annihilates, and a e \ is the elastic cross section. 
At small positron energies, a c \ ~ 47r|/o| 2 , and in the zero- 
energy limit a c \ = Ana 2 , where a is the positron scatter- 
ing length, a — — /q at k = 0. 7 

A simple estimate of the factor Anp e 5R = F in Eq. (|2"8)l 
is obtained using the Ps density at the origin, p e ~ pp s — 
1/8tt, and SR ~ 1, which yields F ~ 0.5. Equation (|2"8|) 
then shows that the magnitude of Zgjj ^ s comparable 
to the geometrical cross section of the target (in atomic 
units), unless a c \ is much greater than R 2 . 

When the scattering cross section is large, the annihi- 
lation rate is greatly enhanced. This occurs when the 
positro n has a virtual or a bound st ate close to zero 
energy ( Goldanskii and Savasovl Il964l ). Such states are 
characterized by a small parameter k = 1/a, \k\ <C R ■ 



It is related to the energy of the bound state eo 
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(for k > 0), or virtual state, eo = k 2 /2 (for k < 0). 
This parameter determines the low-energy s-wave scat- 
tering amplitude /o = — (k + ifc) -1 and cross section 
cr c i ~ 4tt/(k 2 + A: 2 ) (|Landau and Lifshit j . Il977lh For 
small k, this cross section can be much greater than the 
geometrical size of the target. The last term in brack- 
ets in Eq. ([28} t hen dominates, and shows a simi- 



lar en hancement ( Dzuba et al. , 1993; Mit rov and Ivanovl . 
I2002D . 8 
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fc 2 ' 



(29) 



The applicability of Eq. (j29)) is shown in Fig. [2] 
It shows the Z e g values from the Schwinger mul- 
tichannel (SMC) calcula tion for C2H2 and C2H4 
(|do N. Varella et all [2002h . fitted using Eq. (gSJ) with a 
small vertical offset to account for the nonresonant Z e g 
background. According to the SMC calculation, both 
molecules possess virtual positron states. This results in 
the characteristic rise of Z c g at small positron momenta 
described by Eq. (|29]l. The virtual level in C 2 H 2 (fit- 
ted value K = 0.0041) lies closer to zero energy than in 
C2H4 (k = 0.0372), which manifests in the large Z e ff 
values for acetylene. The fitted factor F ps 0.25 for 
the two molecules is smaller than the estimate obtained 
from high-quality atomic calculations (see below). This 
is likely an indication of the lack of short-range correla- 
tion terms in the SMC calculation, which would enhance 
the electron density at the positron (see Sec. III.HTTT) . 



1000 



N 



100 



L0 




0.01 0.1 1 

Positron momentum (a.u.) 

FIG. 2 Comparison of the Z e a values calculated using the 
SMC method for C2H2 (squares) and C2H4 (circles) by 
Ido N. Varella et al\ (|2002ft (see Sec. III.H.1|) with the fit us- 
ing Eq. (|29l) with a constant vertical offset. Parameters of 
the fit: C 2 H 2 (solid curve), F = 0.261, K = 0.0041; C2H4 
(dashed curve), F = 0.230, « = 0.0372. 

Positron virtual states explain the large thermal Z e g 
values observed at room temp eratures in heavier noble 
gases dDzuba e£ aEl.ll996lll993h. Th e value of Z cS = 400 
observed for Xe ( Murphy and Surkolll990h is close to the 
maximum direct annihilation rate for thermal positrons 
at 300 K. It is estimated from Eq. (J2SJ to be 



Z ( c s r) < 10 3 . 



(30) 



Higher Z e g values observed in many polyatomics (see, 
e.g., Table Q} can be understood only by considering 
positron-molecule binding and resonances. 

The annihilation rate for the positron bound to an 
atom or molecule is 



7 If the target molecule has a permanent dipole moment fi, the 
long-range dipole potential fj, ■ r/r 3 dominates the low-energy 
scattering (Fabrikant, 1977). This makes <r G i infinite, while Z e g 
remains finite, making Eq. H28H invalid. 

8 The long-range polarization potential — a<j/2r 4 modifies the 
near- threshold form of <r c i and Z^ lr ' ( Gribaklnl. l200d : iMitrovL 
2002), but Eq. I I29II can still be used as an estimate. 



i=l 



P" : : "/;;<- / y^5(r-r i )\V (r 1 ,...,r z ,r)\ 2 dr 1 ...drzdr , 

(31) 

where ^0 is the wave function of the bound state. For a 
weakly bound state (e.g., <C 1 eV) T a can be estimated 
in a way similar to that used for Z^ 1 ^ above. When the 
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positron is outside the target (r > R), \l/o takes the form 



f (ri,---,rz,r) ~ $ (n, . . . , r z )— e 



(32) 



where A is the asymptotic normalization constant. 9 For 
weak binding (k <C R^ 1 ) the main contribution to the 
normalization integral 



|*o(ri, • • • , r z , r)| 2 dri . . . dr z dr = 1, 



comes from large positron separations where Eq. 
valid. This yields 



(33) 



31 is 



(34) 



By matching the wave function V&o in Eq. (|31l) at r = R 
with the asymptotic form in Eq. (|32[) . one obtains 



pa 



nr^c 4Trp e 5R\ A\ 2 = nr^cF—, 

Z7T 



(35) 



(|Gribakinl . l200lh . Hence the electron-positron contact 
density from Eq. ([2"4"[) is estimated by 



Pep - {F/2it)k. 



(36) 



Equation ([35)) shows that T a is proportional to K = 
■\72 gf, (i.e., to the square ro ot of the binding energy; 
see iMitrov and Ivanovl ( 20021 ) for an alternative deriva- 
tion). 10 This relationship between T a and k is confirmed 
by po sitron-atom bound state calculations ( Mitrov et all 
2002). Figure [3] shows values for six atoms with Ei > 
Ep s , namely Be, Mg, Cd, Cu, Zn and Ag, obtained us- 
ing high-quality configuration interaction and stochastic 
variational methods (see Sec. III.HTSI for details). Note 
that the datum for the LiH molecule also follows this 
trend , in sp ite of its large dipole moment [/j, — 5.9 D 
( Lidd . l2000f) ] and relatively strong binding. A linear fit 
through the atomic data points gives a value for the fac- 
tor F — 4irp e 5R in Eq. (|35|) . namely F « 0.66 a.u., 
which is close to the rough estimate given above. One 
can use this value to evaluate the annihilation rates for 
positron-molecule bound states from Eq. (|35|) . provided 
their binding energies are known. 



E. Resonant annihilation 

The effect of resonances on Z e g is described by the 
second and third terms in Eq. (1231) . It is dominated by 
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FIG. 3 Dependence of the annihilation rate F a for positron 
bound states on the para meter k = \Z2ea: solid circles, re- 
cent results for six atoms (iBromlev and Mitrovl . [2002. 2006, 
120101 ; IMitrov et all 12002], 120081); open circles, earlier re- 
sults for these atoms (IMitrov et all 120021 ) and LiH molecule 
(|Mitrov and Rvzhikhl l2000h ; dashed line is the fit F a = 
5.3k (in 10 9 s _1 ) which corresponds to F = 0.66 a.u. 



the diagonal part of the double sum in the last term. 
The off-diagonal and interference terms vanish upon av- 
eraging over the positron energy and can usually be 
neglected. The resonant contribution to the annihila- 
tion c ross section is described by the Breit-Wigner for- 
mula ( Landau and Lifshitzl . I1977I ). 



+ ±r 2 



(37) 



where b v is the degeneracy of the vth resonance. Equa- 
tion (jTj) , ([23]) . and ([37[) then give the resonant Z e s, 
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'Pep 



E 



(38) 



The contact density p ep can be estimated from Eq. (|36j) 
if the positron binding energy is known. To calculate 

Z^g s ^ , one also needs the energies and widths of the reso- 
nances. The former are determined by the positron bind- 
ing energy and the vibrational excitation energies of the 
positron-molecule complex. The elastic and total rates 
depend on the strength of coupling between the positron 
and the vibrational motion, and for overtones and combi- 
nation excitations, on the strength of anharmonic terms 
in the vibrational Hamiltonian. This makes an ab initio 
calculation of resonant Zeff a multifaceted problem. 



9 Equation (1 32 I t assumes that the ionization potential of the atomic 
system satisfies Ei > -Ep s . For Ei < Ep s the asymptotic form is 
that of Ps bound to the positive ion iMitrov et ail 12002ft . 
10 Equation (132 I t is valid if the positron-target interaction is short- 
range. It must be modified if the molecule has a dipole moment 
[see, e.g.. iFabrikanlj Jl977)]. However, Eq. (135 I t can be used as 
an estimate if the dipole force does not play a dominant role in 
the binding. 



F. Resonances due to infrared-active modes 

One case in which such a calculation is possible is that 
of is olated vibrational r e sonanc es of IR-active fundamen- 
tals ( Gribakin and Lee] , l2006b[) . Consider a small poly- 
atomic molecule which supports a bound positron state 
with a small binding energy = k 2 /2 <C 1. The wave 
function of the bound positron is very diffuse. Outside 
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the molecule it behaves as ifo(r) = Ar~ e~ Kr , with A 
given by Eq. (jM)) . 

Suppose that the vibrational modes in this molecule 
are not mixed with overtones or combination vibrations. 
Due to the weakness of the positron binding, the vibra- 
tional excitation energies of the positron-molecule com- 
plex are close to the vibrational fundamentals cj„ of the 
neutral molecule, E v ss u) v . In this case the sum in 
Eq. (|3"8"f is over the modes v, and the resonant energies 
are e v — uj v — Eb- Some (or all) of these modes can be IR 
active. Positron capture into the corresponding VFR is 
mediated by the long-range dipole coupling, and one can 
readily evaluate this contribution to Z^g . 

Consider a positron with momentum k incident upon a 
molecule in the vibrational ground state $o(R)> where R 
represents all of the molecular coordinates. If fc 2 /2 w e„, 
the positron can be captured into a VFR, and thus 
be bound to the molecule in a vibrationally excited 
state "^(R). The corresponding rate given by 
Eq. (|22|) . can be found usin g a method similar to the 
Born dipole approximation ( Land . [l980) with the cou- 
pling V — d • r /r 3 , where d is the dipole moment operator 
of the molecule. 12 This gives the amplitude 



(r) $: (R) ^ e* kr $ (R) drdR 




(39) 



where = ("J^dl^o), and 2F1 is the hypergeometric 
function. 13 The corresponding elastic rate is 
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(40) 



-1/2 



1F1 



(*■ 



• 5. 

' 2' 



1 — Sb/ijJu-, with a maximum 



where h(£) = £ 3/2 (l - £)~ 
dimensionless function of £ 
h « 0.75 at f w 0.89. 

Equation (|40p shows that the elastic rate for the reso- 
nance of an IR-active mode is determined largely by its 
frequency ui v and transition dipole amplitude d v . Their 
values are known for many species from IR abso rption 
measurements [e.g., sec iBishop and CheunsJ (|l982t) ]. 

This t heory has been succ e ssfully applied to the methyl 
halides ( Gribakin and Leel l2006bl). Energy- r esolve d 



measurements of Z e g for CH^F (jBarnes et all , l2003h . 



11 There is extensive experimental evidence that E v ss ui u for most 
resonances observed (cf. Sees. HVI andrVTl. Apparent exceptions, 
where shifts ~ 10-20 meV are observed are the C-H stretch mode 
of CH3F and the O-H stretch in methanol (cf. Sec. IIVI . 

12 Experiments show that the Born dipole approximation provides 
good estimates, or lower bounds, for the exci tation cross sections 
of IR-active modes by lo w-energy positrons llMarler et all |2006| ; 
iMarler and Surkd . I2005T) . 

13 Here, 2 Fi (§, 1; §; -z 2 ) = §z- 2 [(l + z^z" 1 arctan z - 1], 



CH3CI, and CH 3 Br (jBarnes et all 120061) show peaks 
close to the vibrational mode energies, pointing to a 
sizeable contribution of resonant annihilation in these 
molecules. The methyl halides have C% v symmetry, and 
all six vibrational modes are IR active (see Table [TT| for 
CH3CI). Methyl halides are also relatively small, so that 
IVR effects may not be important in the energy range 
of the fundamentals (see Sec. III.G[) . Thus for these 
molecules, Eqs. (j36j), (|38j) . and (l40l) allow one to cal- 

(res) 
eff • 



culate the contribution of all VFR to Z, 



TABLE II Characteristics of the vibrational modes of CH3CI. 



Mode 


Symmetry b v 


U3 v a 










(meV) 


(a.u.) 


(MeV) 


v\ 


A 1 1 


363 


0.0191 


57.2 


;/ 2 


At 1 


168 


0.0176 


22.9 


vz 


Ax 1 


91 


0.0442 


65.7 


7/4 


E 2 


373 


0.0099 


15.9 




E 2 


180 


0.0162 


20.9 


v$ 


E 2 


126 


0.0111 


6.4 


"Mode 


energies and 


transition 


amplitudes 


from 


IBishop and Cheund (|l982h: 
£5 = 25 meV. 


elastic widths 


from Ea.BOt for 



The only free parameter in the theory is the positron 
binding energy. It can be chosen by comparison with the 
experimental Z c s spectrum. To do this, the theoretical 
Z e g must be averaged over the energy distribution of the 
positron beam /b(e — e) from Eq. (|56j) . 



e)de. 



(41) 



This integral is simplified by the fact that the resonances 
are very narrow, since the total width T u = Tf, + is 
small compared to the energy spread of the beam. For 
example, the values in Table [TTI show that Tf, < 0.1 meV. 
The annihilation width is even smaller; for Sb = 25 meV, 
Eq. ([35]) yields T£ = 0.15 ^eV. Hence, 

h V e 



^eff 



(e) = 2tt p ep 



-A(e-e„), (42) 



where k v = \f2e v , and A(^) = ft>(—E) describes the 
shape of the resonances as measured in the positron-beam 
experiment (see Sec. IIILPp . 

The above estimates show that <§; T e v . In this case 
the total decay rate is dominated by the elastic rate, 
r„ « r^, and the contributions of individual resonances 
to the sum in Eq. (|4"2"|) are not sensitive to the precise 
values of Y e v . Therefore even relatively weak positron- 
vibrational coupling is sufficient to fully "turn on" the 
resonance contribution. This explains why the Z e g and 
IR absorption spectra of a molecule can be quite differ- 
ent, even when the resonant Z c r is determined by dipole 
coupling. 

Application of this theory to methyl halides and other 
small molecules is discussed in Sec. IIVI 
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G. Resonant annihilation in large molecules 

1. Vibrational level densities 

In general, a molecule with N atoms has 3 N — 6 vibra- 
tional degrees of freedom. If positron attachment pro- 
ceeds only via excitation of the single-mode VFRs, the 
resonant Z e g values will grow linearly with the size of 
the molecule. However the experimental Z e g data show 
a much faster increase (cf. the data for alkanes in Ta- 
ble HI . These large Z e g values can only be explained if 
the positrons can couple, at least indirectly, with multi- 
quantum vibrations. 

In large polyatomic species, the total vibrational level 
density increases rapidly with the excitation energy and 
is quite high, even in the energy range of the fundamen- 
tals. The spacing between the multimode VFRs in such 
spectra is small compared with the energy spread of the 
incident positrons. Averaging from Eq. (f3"5|) over 

an energy interval which contains many resonances near 
positron energy e, one obtains 



Z eS {e) 



2ir 2 Pep / r e (e) 



p(s + £ b ), 



(43) 



where p(e + £&) is the vibrational level density of the 
positron-molecule complex. From now on we omit the 
superscript in zjfg, since Z e {i in large molecules is al- 
most entirely due to resonant annihilation. In Eq. (|43[) 
it is assumed that the positron collides with a molecule 
in the ground vibrational state. Larger molecules have a 
significant thermal energy content at room temperature, 
which can be taken into account (see below). 

The Z c ff in Eqs. (|3"5|) and (1431) . and the contributions 
of individual resonances in Eq. (j42j) , are proportional to 
Pep/k. Since p ep oc ^fei [cf. Eq. (|36|)]. the resonant 
contribution to Z e g is pro porti onal to the dimensionless 
"kinematic" factor g = y^Sb/s. The magnitude and en- 
ergy dependence of Z e g beyond this factor are due to the 
dynamics of molecular vibrations and positron interac- 
tion with them, etc. Normalizing Z c ff by g has proven to 
be use ful in analyzing va r ious t rends in resonant annihi- 
lation (jYoung and Surkol . 120071 ) (see Sees. HVl and IV)). 

According to Eq. (j43|) . the resonant Z c g is propor- 
tional to the vibrational spectrum density p. This density 
can be evaluated easily in the harmonic approximation 
in which E v — ^2 k rikUJk, where rik are non- negative inte- 
gers, and u>k are the mode frequencies. For weakly bound 
positron-molecule complexes, these frequencies are close 
to those of t he neutral molecule, and can be taken from 
experiment ( Linstrom and Mallardl . 120051 ) or quantum 
chemi stry calculations (e.g., using Q-Chem. iKong et all 
f|2000h V To apply Eq. g3]) to the annihilation of thermal 
positrons at temperature T, the density must be averaged 
over the Maxwellian positron energy distribution, 



where Ub is the Boltzmann constant. Figure [4] shows 
p(eb) for alkanes with n = 3-9 carbons. These densities 
increase rapidly with the positron binding energy and 
with the size of the alkane molecule. 
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FIG. 4 Vibrational level densities in alkanes as functions of 
the positron binding energy: solid curves, calculated from 
Eq. (|44[) for room-temperature incident positrons; circles, es- 
timated from experimental Z c s values at 300 K for C3H8 to 
C9H20, using Eqs. (|36[) and (|43|) and assuming F e /r = 1, and 
plotted against the experimental binding energies (Table [Vj. 



The estimates made in Sec. III.FI show that r e /r m 1 
for the VFRs of IR-active modes. If the same were true 
for the multimode resonances, the increase of Z e g values 
along the alkane series would match the growth of their 
vibrational level densities, as per Eq. (|4"3"]l . However, ex- 
perimental room-temperature Z c g values increase at a 
much slower rate. As a result, the effective vibrational 
densities p estimated using Eqs. (|43l) and (|3^|) from exper- 
imental thermal Z c g values and binding energies (Table 
|V|) , and assuming r e /r = 1, are much lower than the cal- 
culated densities for all alkanes larger than propane (see 
Fig. 2]). Hence the above ass umption is generally inc or- 
rect, and instead, r e /r < 1 (jGribakin and Gilll200l . 

What is the physical reason for the suppressed r e /T 
ratio? First, the positron coupling to multimode VFR is 
likely much weaker than that of single-mode resonances. 
For example, the dipole coupling analyzed in Sec. III.FI 
can cause only single-quantum vibrational transitions (in 
the harmonic approximation). Combination vibrations 
and overtones can be excited due to anharmonic or Corio- 
lis terms in the vibrational Hamiltonian. In this case, the 
coupling strength of a single-mode excitation is divided 
between many multiquantum excitations. The value of 
r e for a multimode VFR is then only a small fraction 
of the typical single-mode T e . This is discussed in more 
detail in Sec. III.G.2I In addition, some vibrational ex- 
citations may be completely inaccessible to the incident 
positron (e.g., due to symmetry). This will reduce the av- 
erage (r e /r), or, effectively, reduce the vibrational den- 
sity in Eq. (|33")) . 

Another reason for the reduction of r e /T could be the 
contribution of vibrationally inelastic escape to the total 
width r. At present there is little direct experimental 
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evidence of inelastic positron escape (see Sec. E}. How- 
ever, if the system was in the regime of strong mixing 
between the single-mode and multimode vibrational ex- 
citations, then one cannot see why such channels would 
be closed. 14 

A simple estimate of this effect can be made assum- 
ing compl ete, statistical mixing o f all vibrational ex- 
citations (jGribakin and Le"el 120091 ). When a positron 
with energy e collides with a molecule with vibrational 
energy E v then all final vibrational states with ener- 
gies E' v < E v + e can be populated. Assuming that 
the positron coupling strengths to all vibrational exci- 
tations are similar, one has r e /r « 1/N(e + E v ), where 
N(e + E V ) = " p(E' v )dE' v is the number of open vi- 
brational escape channels. Equation (|4"3")) then becomes 



■> off 



2-K 2 p ep p(e ■ 



E v + Eb) 



N(e + E v ) 



(45) 



For a given molecule this expression contains only one 
free parameter, namely, the positron binding energy. 

The binding energies can be chosen by comparison 
with experimental room-temperature data for thermal 
positrons (Table |V| > averaging Z e s from Eq. (|45j) over 
the initial target states using the Boltzmann factors 
exp(— E^/Zc^T), and the Maxwellian positron energy dis- 
tribution. For alkanes with 3-8 carbons such a fit gives 
the binding energies Eb — 22, 42, 65, 90, 103 and 
122 meV, respectively. These values are in good agree- 
ment with those measured in the positron-beam experi- 
ments (cf. Table fV| . However, the dependence of Z e g on 
the positron energy predicted by Eq. (|45[) is in striking 
disagreement with the measured energy-resolved Z e g, as 
shown in Fig. [5] for butane and octane. The Z c g val- 
ues from Eq. (|45p decrease monotonically, the increase in 
N(e + E v ) being faster than that of the density p. In con- 
trast, the experimental Z g spectra show resonant peaks 
corresponding to the vibrational fundamentals (down- 
shifted by Eb, cf. Fig.Q]), and the peak Z D g values exceed 
the predictions of Eq. (|45l) by a factor of ~ 50. Thus 
the model that assumes complete statistical mixing with 
all vibrational excitations coupled to the positron con- 
tinuum, does not explain either the resonant structure or 
the large values of Z e g that are observed. 



2. Mode-based resonant doorway states 

The energy dependence of the measured Z c g in alka- 
nes shown in Fig. [5] is of the type expected for mode- 
based VFRs described in Sec. III.FI On the other hand, 
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FIG. 5 Comparison of experimental Z e g for butane (solid 
symbols) and octane (open symbols) with the predictions of 
the statistical model: diamonds, experimental Z^g values for 
thermal positrons at 300 K (cf . Table [Vj) ; circles, Z e s mea - 
sured as a function of positron energy (jBarnes et aU[2003T ); 
curves, Z e fr calculated from Eq. (|45[) . using Eb fitted to repro- 
duce thermal Z e g. 



the magnitude of Z e g and its increase with the size of 
the molecule, are much greater than the contribution of 
mode-based resonances. This suggests a two-step model 
of positron capture tha t involves mode-based v ibrational 
doorway resonances (jGribakin and Gilll . [20041 ) . 

In this model, the incident positron first forms a qua- 
sibound state with the molecule by transferring its ex- 
cess energy to a single mode with near-resonant energy 
Ld n ~ £ + Eb- This simple doorway state of the positron- 
molec ule complex [also terme d a zeroth-order "bright" 
state ( Nesbitt and Fieldl . ll996l) ] is embedded in the dense 
spectrum of multimode vibrations (or "dark" states, as 
they are not coupled directly to the positron continuum) . 
Due to vibrational state mixing caused by anharmonic 
or rotationally-induced coupling terms in the vibrational 
Hamiltonian, the doorway state can then "spread" into 
multimode vibrational states. This process of vibra- 
tional energy redistribution takes place on a time scale 



t ~ l/r spr , where T spl - is known as the spreading width 
in nuclear physics, or the IVR rate in molecular physics. 

To link the multimode- VFR and the doorway-state- 
resonance pictures, consider a perturbative expression for 
the clastic rate, 



r^ = 27rK^|v/|o, £ )| 5 



(46) 



where |0,e) describes the positron incident on the 
ground-state molecule. 16 The multimode eigenstate of 
the positron- molecule complex can be expanded in 



Vibrationally inelastic scattering following capture has been 
studied in electron-molecule collisions via energy-loss spectra. 
In addition to the excitation of modes and overtones, vibrational 
state "quasicontinua" have been observed in larger polya tomic 
molecules at electron-volt excitation energies llAllanl fl 984fl . 



The term "doorway resonance" originates in nuclear physics, 
where it means "a metastable state formed in the initial state 
of the reaction" , which "may decay partly into the open chan- 
nels (direct reactions), and partly through the coupling to the 
internal degrees of freedom" llBohr a nd Mottclson, 1998). 
See Eq. (1221 ; for simplicity, in Eq. I|46|l we assume that only one 
positron partial wave contributes to the rate and that its wave 
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the basis of the noninteracting (e.g., harmonic approxi- 
mation) multimode vibrational states 



4\ 



(47) 



where the coefficients c\ are obtained by diagonalizing 
the vibrational Hamiltonian in the basis of |3>»). 

Let us assume that of all only those which de- 

scribe a bound positron and a single-mode excitation 
(i.e., \n, £o)) where n indicates the mode) are coupled to 
|0, e). These states |n, £q) are the doorway states intro- 
duced above. The coefficients c\ u ' describe the mixing 
of the doorway states with the multimode eigenstates v 
(i.e., "spreading"). The corresponding probabilities can 
be approximated by a Breit-Wigner line shape, 



r 2 

spr 



/4 



(£„-£i) a + r» /4 



subject to normalization J^IC^I 2 = !• Here £^ and 
£/, are the energies of the eigen- and basis states of the 
positron-molecule complex, respectively, so that E v ~ 
e - e = e + e b and E t = uj n for = |n,e ). 

Using Eqs. (j43]), 06}, (gTJ and (g8j), one obtains Z eff , 
averaged over energy on the scale of closely spaced VFRs, 
as 



(48) 



Z, 



cff 



2n 2 p ep T s 



2 7 rr( £ )^( £ -a;„ + £6 ) 2 + ir| ] 



(49) 



pr 



where = 27r|(n,eo|^|0,e)| 2 is the elastic rate of the 
doorway of mode n. Equation (|49[) has the same form as 
Eq. (|38|) . except that it contains the elastic rates of the 
doorways, T^, and the sum is over the modes. 

Comparison with Eq. (|43|) shows that the resonant 
energy dependence of Z e g in Eq. (|4"9"jl is due to the mod- 
ulation of the elastic rate by the mode-based doorways. 
Evaluating Z e g from Eq. (|49l) at the doorway resonance 
energy (i.e., e — io n — Eb) and comparing it with Eq. (|4"3")l 
gives an estimate of the elastic rate of the VFR: 



r e (e) 



2H 



irp(e + e b )T s 



(50) 



The product p(e + £&)r spr is the number of vibra- 
tional eigenstates within the energy interval r spr . In 
the regime of strong level mixing this number is large, 
p(e + £&)r spr 3> 1. Equation ([50]) shows that the elas- 
tic rates of the multimode VFRs are much smaller than 
those of the doorways, as each VFR carries a small frac- 
tion of the positron coupling strength of the doorway. 

Equation (|49|) shows that the enhancement of Z c g for 
larger molecules is caused by (i) larger binding energies 



function is normalized to a (^-function in energy rather than to a 
plane wave. 
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FIG. 6 Comparison of the experimental Z e g for (•) bu- 
tane, C4H10 and (o) octane, CsHig (scaled by a factor 1/50) 
(jBarnes et aZ.L 12003 ). with Z cB from Eq. |Q5}. The theoret- 
ical Z e a is averaged over the positron beam energy distri- 
bution. The fit for butane uses e& = 35 meV, T%/T = 7.2 
for C-H stretch modes and T^/F = 1.2 for the rest; for oc- 
tane, Eb = 122 meV, T e n /Y = 84 for C-H stretch modes and 
F e „/r = 14 for the rest. 



and p e p oc y/eb, (ii) larger number of modes/doorways in 
the sum, and most importantly (iii) smaller total decay 
rates r(e) (i.e., longer lifetimes) of the VFRs. This de- 
crease in T(e) must be related to the suppression of T e 
due to higher level densities. The total rate is, however, 
always bounded from below by the annihilation contri- 
bution, T(e) > T a . The annihilation rates are typically 
smaller than the elastic rates of mode-based doorways 
by 2-3 orders of magnitude (see Sec. III.Fj) . Hence, one 
can expect a similar-sized increase in the contribution to 
Z e g from each mode-based doorway resonance in a large 
molecule, compared with that in a small molecule. 
Figure [5] shows the applicability of Eq. (141)1) . Here it is 



used to fit the exper imental energy-resol ved Z c g data for 
butane and octane ( Barnes et al\ . [2003h . For this com- 



parison, Z e g from Eq. ((49]) is averaged over the positron 
energy distribution (Sec. IIII.D|) . The result does not 
depend on the spreading width, as long as it is much 
smaller than the typical energy spread of the beam (i.e., 
r sp r *C 40 meV) . The binding energies of the two species 
are chosen to be e& = 35 and 122 meV. The remain- 
ing unknown parameter is the ratio r^/r(e). Given the 
differences between the C-H stretch and lower-energy res- 
onances in Fig. [51 two different values are used for these 
groups of modes, the former six times greater than the 
latter. The ratios T^/T for octane are approximately 
12 times larger than for butane. In the context of the 
model, this reflects the greater degree of mixing between 
multimode VFRs in the larger molecule. 

The results shown in Fig. |6] summarize the extent of 
our understanding of annihilation in large molecules. To 
explain the observed magnitudes of Z e g, positron cap- 
ture into multimode VFR must be invoked, mediated by 
a process such as IVR. However, if the enhancement is 
indeed due to IVR, it appears to be far from statistically 
complete. 
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3. Annihilation and the onset of IVR 



known, are also listed. 



Experimentally, resonant annihilation is observed in 
relatively small polyatomics with four or five atoms 
fSec. HV|) . as well as in much larger molecules (Sec.[\~). In 
the first case, the annihilation is due to the single-mode 
VFRs, possibly augmented by contributions of combina- 
tion and overtone resonances (see Sec. HV.Bj) . In large 
molecules enhanced Z e g values are assumed to be due to 
the spreading of the vibrational energy into multimode 
VFRs. The energy dependence of Z c g in large molecules 
is typically well represented by the spectrum of the fun- 
damentals, which act as doorways. Isolated overtone and 
combination VFRs cannot, in general, be identified as 
distinct doorways in the Z e R spectrum, though there are 
exceptions (e.g., benzene, Sec. |Vj). This can probably be 
explained by the much smaller ratios of /T(e) for the 
combination/overtone doorways [cf. Eq. (|49j) ]. 

The transition from the small- to the large-molecule 
behavior occurs with the onset of strong vibrational mix- 
ing (i.e., IVR) involving, in particular, the vibrational 
fundamentals. The phenomenon of IV R has been studied 
widel y using a number of techniques (jNesbitt and Field! 
1996). In particular, high-resolution (0.0005 cm -1 ) mea- 
surements of vibrational spectra of jet-cooled molecules 
allow direct observation of the splitting of single-mode 
transitions into clumps of vibrationally mixed multimode 
levels s pread over an energy interv al r spr ~ 0.02 cm -1 w 
2 ^eV |McIlrov and Nesbittl . Il990l) . 

Another technique used to study IVR in the 
range of the C-H stre t ch m odes, is IR fluorescence 
([Stewart and McDonald! Il983l) . Here, a fundamental vi- 
bration is excited by a short laser pulse. If this vibra- 
tion is mixed with multimode excitations, the amount 
of fluorescence at the fundamental frequency is reduced. 
The measured fluorescence, normalized to the known IR 
absorption strength, gives the "dilution factor" . Its re- 
ciprocal characterizes the number of eigenstates that are 
strongly coupled to the fundamental. If all vibrational 
states within the energy range r spr are mixed, the dilu- 
tion factor will be ~ (/cTspr)" 1 , where p is the vibrational 
spectrum density at the re levant energy. Note that asim- 
ilar factor enters Eq. ([50]) ■ IStewart and McDonald! (|l983l ) 
observed that the dilution factor drops rapidly when the 
level density is increased beyond the "threshold" value of 
p = 10-100/cm -1 . This is in agreement with the value 
r sp r ~ 0.02 cm -1 , as pT spl > 1 marks the onset of IVR. 

In order to see if a similar threshold governs the transi- 
tion from the small- to the large-molecule behavior in an- 
nihilation, vibrational densities have been evaluated for 
a number of molecules using the harmonic approxima- 
tion. The densities at the excitation energy of the C-H 
stretch mode (E — 2900 cm -1 = 0.36 eV) are shown 
in Table IIIII To assess the effect of finite molecular 
temperature, the densities were evaluated both at zero 
temperature, denoted p(E), and at room temperature, 
Pt{E) = E, P{E + E v )e- E v/ k * T /J2 v e- E *' k * T . In Ta- 
blc lllll values of the IR fluorescence dilution factor, where 



TABLE III Vibrational excitation densities for zero- and 
room-temperature molecules at the C-H stretch mode energy. 



Molecule 




P a 


pT 


Dilution 






1/cm -1 


1/cm -1 


factors' 


Acetylene 


C2H2 


0.05 


0.06 


- 


Methyl chloride 


CH3CI 


0.07 


0.07 




Ethylene 


C2H4 


0.08 


0.09 




Methanol 


CH3OH 


0.12 


0.13 




Ethane 


C 2 H 6 


0.42 


0.54 


0.8, 1 


Cyclopropane 


C3H6 


0.59 


0.81 


0.3, 0.7, 0.8 


Ethanol 


C 2 H 5 OH 


2.84 


6.22 




Propane 


C3H8 


5.19 


14.3 


0.2, 0.4 


Benzene 


C 6 H 6 


5.76 


22.3 


0.6, 0.7, 0.9 


Butane 


C4H10 


106 


921 


0.16, 0.05 c 



"Densities at E = 2900 cm- 1 for T = and T = 293 K, averaged 
over a Gaussian wi th 10 meV, FWHM. Mode fre quencies for most 
molecules are from |T~n strom and Mallard (2005V 

^Va lues for the C-H stretch modes {Stewart and M cDonald! . 
11983ft . 

c Values for 1-butyne jKim et all , Il987ft and isobutane 
( Ste wart and McDonald! . Il983ft . respectively. 



Energy-resolved Z c g measurements put the top four 
molecules in Table IIIII in the small-molecule category (cf . 
Sec. II VI) . Propane, benzene and butane at the bot- 
tom of the table, behave as large molecules (Sec. |VT) . 
as they exhibit Z e ff values that cannot be explained by 
the mode-based VFR. Ethane and cyclopropane appear 
to be borderline, with characteristics that place them 
in both categories (e.g., the C-H stretch peak in both 
molecules appears to be strongly enhanced). The val- 
ues of the densities in Table IIIII are broadly in agree- 
ment with this classification, with the threshold density, 
p ~ l/cm -1 . This value is lower than the threshold den- 
sity in the fluoresecence studies, perhaps due to differ- 
ences in the nature of the vibrational energy transfer in 
the two cases. Somewhat surprisingly, Z e f{ measurements 
for ethanol (cf. Sec. IIV.D|) indicate that it is a "small 
molecule" , in contrast with its vibrational density value. 
This suggests that the threshold value depends on the 
details of the vibrational Hamiltonian. 

Table IIIII shows that molecular temperature can 
strongly affect the density for larger molecules. The 
only experiments to investigate the effect of molecular 
temperature on Z e g in similar-sized molecules, heptane 
and pentane, show only a small effect at the C-H stretch 
peak (Sec. IV.E~2)) . While the analysis here indicates that 
molecules with p ~ l/cm -1 might change from exhibit- 
ing small- to large-molecule behavior as the temperature 
is increased, molecules with these density values exhibit 
only small changes in p with temperature. Thus the com- 
bination of this analysis and the experimental results on 
large molecules indicate that temperature likely plays a 
limited role in changing the details of the Z c g spectra. 
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When the level densities for larger molecules are ex- 
amined as a function of the excitation energy E, the 
densities in the C-H stretch energy range are an order-of- 
magnitude greater than those in the range of other funda- 
mentals, E < 1000 cm -1 . In principle, this could explain 
the stronger enhancement of Z e g in the C-H stretch peak, 
as compared to that in the low-energy mode peaks (see 
Fig. |5]and Sec. |V] for further discussion). 

To summarize, optics-based IVR studies show that the 
degree of IVR increases rapidly as a function of molec- 
ular size. The onset of strong IVR occurs for molecules 
that are similar in size to those for which mode-based, 
VFR-mediated annihilation begins to fail to explain the 
Z e g spectra (e.g., ethane, cyclopropane and propane). 
Thus, while indirect, these results fit well with the phys- 
ical picture that IVR is responsible for the very large 
annihilation rates observed in large molecules. 



H. Calculations of annihilation and binding 

I. Annihilation 

Most calculations of Z c g for molecules have been done 
for diatomics or small polyatomics with fixed nuclei, ig- 
noring the vibrational dynamics. This is a good approx- 
imation for molecules that do not possess VFRs, where 
direct annihilation is the dominant mechanism. Such a 
calculation is still far from trivial, given the large role of 
electron-positron correlations. 

The H2 molecule is an example. In this case one 
can construct a sufficiently flexible trial wave function 
for the positron and two electrons and use the gen- 
eralized Kcin method to solve the scattering problem 
(lArmourl Il984l) and calculate Z e g ( Armour and Bakerl 
I1985D . This calculation provides an accurate description 
of the positron elastic scattering cross section below 5 eV 
([Armour et aLl . ll990D . Calculation of Z e g requires the in- 
clusion of terms wi th explicit dependence on t he electron- 
positron distance ( Ar^uL^nTBake^, [HH]) . However, 
it still yields a thermal room-temperature annihilation 
rate of Z c g — 10.2, well bel o w the experimental value 
of 14.8 ± 0.2 (jMcNutt et all 11979ft . This discrepancy 
has been finally resolved usi ng the stochastic variational 
method ( Zhang et all [2009), which is one of the most 
powerful methods for studying few-body systems. 

Another ab initio method used to calculate positron- 
molecule scattering and ann ihilation is the Schwinger 
multichannel (SMC) method ([Germano and Limal I1993T 
Ida Silva et a/.l . [l994l ). In this scheme the (Z + 1) -particle 
wave function is expanded in a Cartesian Gaussian ba- 
sis set, with the functions centered on the atomic nuclei 
and on additional centers outside the molecule. The lat- 
ter are important for representing electron-positron cor- 
relation effects. This method has been a pplied to H2 
Lino et all Il998l). N? (|de Carvalho et all I2OOOI), C9H, 



tic scattering cross sections are in good agreement 
with the experimental data. In contrast, the calcu- 
lated Z s ff_j2foe^_Jc^room ;: temperature positron ener- 
gies (Ido N. Varella et all |2002) are well below the mea- 
sured thermal data, namely, Z c tf — 7.3 vs 14.8 for H 
dMcNutt et a/.l . ll979D . 9.3 ^ 30.5 for N 
[1982J), 73 w~L200 for C->Ra (llwata et all 
vs 3160 for C 2 H 2 (|lwata et all 11997 ' 



Hevl and et al. 
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19941) . and 145 
The discrep- 
ancy for H2 and N2 is most likely due to a lack of basis 
functions that describe short-range electron-positron cor- 
relations. The much larger gap between the theory and 
experiment for ethylene and acetylene has a different ori- 
gin. As discussed in Sec. IIV.C1 the energy-resolved Z e g 
data for both molecules show large contributions of res- 
onant annihilation that cannot be described by a fixed- 
nuclei calculation. 

In spite of the failure to reproduce the Z e g values, 
the SMC calculations provide an important clue about 
the physics of the positron interaction with ethylene and 
acetylene. They show that electron-positron correlations 
in these systems lead to strong positron-molecule attrac- 
tion that produces virtual levels close to ze ro energy 
(|de Carvalho et all 120031 : Ida Silva et ^lll996l ). As a re- 
sult, the calculated Z e g display a characteristic growth 
at low energies (cf. Fig. [2]). Observation of resonant 
annihilation for these molecules (Sec. IIV.CI) indicates 
that they do, in fact, support positron bound states with 
Eh < 10 mcV. This means that a relatively small (though, 
computationally challenging) improvement in the SMC 
calculations could prove the existence of these bound 
states. This would represent an important step towards 
ab initio theoretical description of resonant annihilation. 

Positron interaction with acetylene (and methane) 
was also examined by iNishimura and Gianturcol (|2003l ) 
using a body-fixed vibrational close-coupling (VCC) 
method. In this approach, the positron-molecule inter- 
action is modeled using a correlation-polarization poten- 
tial (CPP). These calculations revealed the existence of 
a positron virtual state for methane (k ~ —0.15 a.u.) 
and acetylene (k w —0.01 a.u.). As in the case of 
the SMC calculations, it appears that only a small in- 
crease in the potential is required to turn the virtual 
state into the observed bound state in acetylene. Such 
an increase may well be within the uncertainty of the 
CPP method. In fact, a fixed-nuclei calculation with 
the existing CPP, showed that the change from a virtual 
state to the bound states can be achieved by stretch- 
ing t he C-H bonds in C2H2, C 2 H4, an d C2H 6 by 15- 
30% ( Nishimura and Gianturcol . l2005bl) . In acetylene 
this change can also be induced by symmetric bendin g 
by about 16 degrees (|Nishimura and Gianturcol . l2004f) . 
However, this mechanism of bound-state formation ap- 



dda Silva et all Il996h . and C2H2 (|de Carvalho et all 
120031) . In all cases the differential and total elas- 



17 The SMC computer code used to calculate the Z B g values prior 
to 2001 contained an extra factor of Z due to a programming 
mistake ((do N. Varella et all I2002T) . which gave an illusion of 
agreement with experiment. 
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pears to be problematic, since the energy required to dis- 
tort the molecule far exceeds the positron binding energy 
and the thermal energy of molecules at ~ 300 K. 

The earliest calculation of positron annihila- 
tion on polyatomi c mole cules is probably that of 
IJain and Thompson! (|l983l ) for CH 4 and NH 3 . They 
described the interaction between the positron and the 
target by means of a potential V s (r) + V p (r), where V s (r) 
is the full electrostatic potential of the ground-state 
molecule, and V p (r) is a CPP. This potential has the 



correct form V p (r) 



-ad/2r at large distances. 



Solving the Schrodinger equation for the positron wave 
function fk(r) yields the scattering amplitude and 
cross section, and <Pk(r) is then used in Eq. (f20]t to 
calculate Z e s- IJain and Thompson! ( 19831 ) showed that 
it is important to go beyond this approximation and 
include the effect of distortion of the electron density 
by the positron. This led to an increase of Z e fi values 
for CH4 by a factor of two. Their final result for CH4, 
Z e g = 99.5 at £ = 0.025 eV, is reasonably close to, but 
still smaller than the experimental room -temperature 
value Z cS = 142 ± 1 (|Wright et all . Il983l) . They con- 
cluded that the remaining discrepancy could be removed 
in a better positron-molecule calculation. More signif- 
icantly, they also concluded that the somewhat large, 
measured Z g value was not the result of the formation 
of a positron-molecule complex. As discussed in Sec. 
III.D1 this is in complete agreement with the current 
understanding of the way in which positron-molecule 
virtual states can enhance low-energy Z c g values. 

An approach similar to that of IJain and Thompson! 
(|l983h has now been tested for a variety of organic and 
inorganic polyatomic molecules ( Gianturco et all . l200lb 
lOcchigrossi a nd Gianturco, I2003D . with ethane and ben- 
zene being the largest. The short-range part of the 
CPP employed in these calculations is based on density- 
functional treatments of the electron-positron correla- 
tions, while the long-range behavior is described an- 
alytically (e .g., as — Qrf/2r 4 , for the dipole polariza- 
tion term) (jjain and Gianturco! . Il99ll ). The positron 
wave function from these calculations is used to cal- 
culate Z e g from Eq. QUI) . More recent calculations 
( Franz and Gianturco! . l2006h 18 employed an enhance- 
ment factor in the integrand of Eq. ([5]). It depends on 
the electr on density and describes its local increase at th e 
positron ( Arponen . ll978l ; [Bororiski and Nieminenl . fl986l) . 

The main result of these fixed-nuclei calculations is 
that they fail to reproduce large experimental Z G ff values 
(> 10 2 ) for molecules such as C2H2, C 2H4, C2H6, and 
C 6 H 6 (jOcchigrossi and Gianturcol 120031) . This is an in- 
direct confirmation of the role of resonant annihilation 
involving nuclear vibrations in these molecules (see Sees. 



IIV.C1 IV.CI and IV.F[) . For smaller polyatomics, such as 
H2O and CH4, the computed room-temperature Z e g val- 
ues, 167 and 65, respectively, are a fact or of two or three 
smaller than the experimental values (jGianturco et all , 
2001). Here an adjustment in the CPP and/or the use of 
the enhancement factor in the calculation of Z e g could 
bring theory and experiment into agreement. Indeed, 
according to the energy-resolved Z c g measurements, nei- 
ther of these molecules shows clearly discernible contri- 
butions of resonant annihilation (Sec. HV.E"| . Similarly, 
for diatomics, such as H 2 , O2, N 2 , NO, and CO, the cal- 
culation with en hancement factors gives Z e g values be- 
tween 10 and 40 ( Franz and Gianturcol . [2006h . which are 
within a factor of two of experimental values for room- 
temperature positrons. 

There are few annihilation calculations that in- 
clude the dynamic interaction between the positron 
and molecular vibrations. The VCC calculations 



18 This work also corrected a normalization error that 
overestimated Z e g values for diatomics reported earlier 



(|Gianturco and Mukheried . I2OO0I) . 



of IGianturco and Mukheried dl999h for C0 2 , and 
IGianturco and MukherieeT (2000) for 2 , N 2 , NO, and 
CO, showed that vibrational coupling has a relatively 
small effect on the annihilation rates. This is to be ex- 
pected, since small molecules such as these most likely do 
not bind the positron (cf. Sec. lVI[) . and so the VFR mech- 
anism is "switched off" for them. The true magnitude 
of the effect of vibrational coupling for these molecules 
remains somewhat uncertain, as the reported theoretical 
Z e g values probably suffer from an uncertainty related to 
the normalization error (see above). It is also surprising 
that the Z e g values obtained in these calculations, us- 
ing a static potential, change little upon inclusion of the 
CPP, since practically all other calculations show that 
correlations have a large effect on positron scattering and 
annihilation. 

The zcro-range-potential (ZRP) model calculations for 
Kr2 (a weakly bound van der Waals dimer) demon- 
strated a numbe r of key features of r esonan t annihilation 
(|Gribakinl . l2002t iGribakin and Led . l2006al ). First, they 
showed that, while the potential representing each of the 
Kr atoms has no bound states, the dimer is capable of 
binding a positron. Second, they showed that positron 
binding has a relatively small effect on the vibrational 
frequency of the complex. It changes by only about 10% 
compared to the frequency of Kx2, in spite of the fact 
that the binding energy is relatively large (i.e., two times 
the vibrational quantum). Finally, this model showed 
the emergence of VFR in a dynamic positron-molecule 
calculation. The resonant contribution leads to a large 
increase in Z e g (e.g., for thermal positrons at 300 K, 
from Z c ff = 250 for direct annihilation to Z e ff = 700- 
950, depending on the details of the model). Unfortu- 
nately, van der Waals molecules such as Kr2 are diffi- 
cult to study experimentally, and the ZRP method is 
in general too crude to predict the binding energies or 
Z e g spectra for molecules for which the annihilation has 
been measured. That being said, this approach does pro- 
vide an easily solvable and instructive model for resonant 
positron- molecule phenomena (see Sec. |II.H"72]> . 
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Recently a 



t heory has been proposed 
(Id'A. Sanchez et all [2009), that described vibra- 
tionally enhanced annihilation using the Feshbach 
projection operator formalism. It assumed that the 
positron is captured into a resonant electronic state that 
determines the subsequent vibrational dynamics. The 
full implications of this theory have yet to be elucidated. 
A similar mechanism drives many electron-molecule 
attachment processes. However, in the case of the 
positron, there is no experimental evidence to date of 
the required resonant electronic states. 



2. Positron-molecule binding 

Calculation of positron binding has proven to be ex- 
ceedingly challenging. The electrostatic interaction be- 
tween positrons and neutral atoms or molecules (with- 
out large dipole moments) is dominated by the nu- 
clear repulsion. At large separations, the electric field 
of the positron gives rise to the attractive — ad/Zr 4 
potential. At short range, there i s an additional at- 
traction due virtual Ps form ation ( Dzuba et all Il995t 
iGribakin and Ludlowl . l2004f) . Together with polariza- 
tion, these forces can overcome the static repulsion and 
thus enable the formation of virtual levels or bound 
states. 

In the case of atoms, reliable calculations for positron 
bound states have been done for about ten species with 
one or two valence electrons (such as Li, Be, Na , Mg, 
Cu, Zn, etc.). The work by IMitrov et all (|2002l ) is a 
good review of the state of the field a few years ago. The 
calculations were done using a variety of methods: many- 
body theory and its combination with the configura- 



tion interaction (CI) (jDzuba et aZ. Ul999lll995h: stochas- 
tic variational method (S VM) ()Rvzhikh and Mitrovl . 
Il997t iRvzhikh et oil . |j998|); and the CI method with 
core polarization potential s ([Bromley and Mitrovl l2000t 
IMitrov and Rvzhikhl . Il999l ). 

These calculations provide useful insights into the 
physics of positron binding. One important parameter 
is the ionization potential of the atomic system Ei and 
its relation to the Ps binding energy Ep s . For systems 
with Ei > Ep s , the electrons are relatively tightly bound 
in the target. Since the positron is repelled by the atomic 
core, it then forms a loosely bound state and stays out- 
side the atom (i.e., represented asymptotically as A+e + ). 
For Ei < £p s however, the positro n can attract a valence 
electr on forming a "Ps cluster" ()Rvzhikh and Mitrovl . 
I1998D . In this case, the bound state is asymptotically a 
Ps atom orbiting the residual positive ion (A + + Ps). Fig- 
ure [7] shows the calculated binding energies Sb for atoms 
as a fu nction of Ej. A cal culation for a model "alkali 
atom" (|Mitrov et all Il999h shown by a dashed curve, 
suggests that Sb peaks at Ei = Ep s , and the calculations 
for real atoms generally support this picture. 

In the "alkali atom" model, a positron and a single va- 
lence electron move in the field of a fixed atomic core. As 
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FIG. 7 Binding energies of positron-atom complexes as a 
function of t heir ionization potential: squares , SVM and CI 



calcu lations (jBromlev and Mitrovl , 120071 . l20icil ; 



Mitro v et al. 



2002); dashed curve, model "alkali atom" (|Mitrov et al. 
19991 ). For positron binding to metastable states the latter 
are indicated in round brackets. The two bound states for Ca 
and Sr are shown in square brackets. 



the core potential is varied, both the ionization potential 
and the di pole polarizab i lity ch ange in such a way that 
ad oc E~ 2 ( Mitrov et all Il999i ). In real systems, and in 
particular, in molecules, ad and Ei can be regarded more 
akin to independent parameters. As shown later (see 
Sec. IVI|) . both are found experimentally to influence the 
observed binding energy. Theoretical studies of atoms 
and experimental studies of molecules indicate that their 
maximum binding energies are comparable, £b < 0.5 eV, 
while the values for the strongly polar alkali hydrides can 
be as large as 1 eV (see below). 

A major difference between positron binding to atoms 
and to molecules is that molecules can have permanent 
dipole moments. Theoretically, a static molecule with 
dipole moment > /i cr = 1.625 D = 0.639 a.u., pos- 
sesses an infinite nu mber of positroni c, as well as elec- 
tronic bound states (|Crawfordl . fl967h . For a molecule 
that is free to rotate, the critical dipole moment for bind- 
ing is greater than /i cr . This critical value increases as 
the molecular moment of inertia decre ases or the an gular 
momentum of the molecule increases (|Garretd . [l97lf) . 

At present, positron-molecule binding has been 
predicted theoretically for a few strongly polar 
molecules (lAdamson et all 120081: iBressanini et all 1998 : 



2006 
2006 



Buenker et al. , 20051 20061: IChoinacki and Strasburge 



Danbv and Tenny so n |. 
Kurtz and Jordaj Il981 



1988; iGianturco et al 



Pichl et al 



2005; 



Mell a et al 



Schrader and WangT 1976 



Schrader a nd Mo xoml. 
Strasburgerl 19961 



2001 



2001 



1999 



20011 . 12004 iTachikawa etall . l2003l l200ll ) . These calcula 



tions employed a variety of methods such as Hartree-Fock 
(HF), CI and diffusion Monte Carlo (DMC), as well as 
explicitly correlated Gaussians (ECG). Positron bound 
states have also been found using the CPP method for 
the nonpolar ca ge-like molec u le C2 0, which has a large 
a d = 25.4 A 3 (jCarev et all , l2008fh A representative 
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selection of recent results is given in Table IIV1 For 
simple molecules, such as LiH, the results obtained 
by different methods are generally in good agreement. 
The binding energies obtained at the static HF level 
increase considerably when correlations are included 
[e.g., see £& values for HCN, urea (NH2)2CO, and acetone 
(CH3)2CO]. In fact, the current CI values for the larger 
polyatomics may still considerably underestimate the 
true £h values due to incomplete CI expansions. 



TABLE IV Calculated positron-molecule binding energies. 



density follows the increase in as described by Eq. (|36|) 
for bound states with nonpolar species. 



Molecule 


» a 


Eb 


Method 


Reference 




(D) 


(meV) 






LiH 


5.88 


1000 


DMC 


Mella et al. (20001 






909 


FCSVM 6 


MR C 






1005 d 


ECG 


Strasburger (2001) 






1043 e 


ECG 


BA C 






1626 


NEO 6 


Adamson et al. (2008) 


LiF 


6.33 


599 


DMC 


Mella et al. (2001) 


BeO 


6.26 


680 


DMC 


Mella et al. (2001) 






375 


MRD-CI 6 


Buenker et al. (2007) 


NaH 


6.40 


1031 


MRD-CI 


Gianturco et al. (2006) 


RbH 


9.03 


1639 


MRD-CI 


Gianturco et al. (2006) 


MgO 


6.42 


472 


MRD-CI 


BL C 


LiO 


6.84 


304 


MRD-CI 


BL C 


CH 2 


2.33 


19 


CI 


Strasbureer (2004) 


HCN 


2.98 


2 


HF 


CS C 






35 


CI 


CS C 






38 


DMC 


Kita et al. (2009) 


Urea 


3.99 


6 


HF 


Tachikawa et al. (2003) 






13 


CI 


Tachikawa et al. (2003) 


Acetone 


2.88 


1 


HF 


Tachikawa et al. (2003) 






4 


CI 


Tachikawa et al. (2003) 


C20 




780 f 


CPP 


Carev et al. (2008) 






230-250 / 


CPP 


Carev et al. (20081 



"Dip ole m oment s from iBuenker et al ] J2005T I: iGutsev et al\ 

l|l997l) ; lLidel 1120001 ) or as cited. 

6 FCSVM, fixed-core SVM; NEO, nuclear-electronic orbital 
m ethod; MRD-CI. mul tir eferen ce single- a nd double-excitation CI. 

iMitrov and Rvzhikbl (l2000h (MR); iBubin and Adamowic"3 
j2004t) (BA): IBuenker and Liebermannl l|2008l ) (BL); 
IChoinacki and Strasbureed ||2006| ) (CS). 

d Adiabatic positron affinity. 

e Non-Born-Oppenheimer variational calculation. 
•'The C % isomer of C20 is predicted to have a deeply bound s-type 
state and three weakly bound p-type states. 



The positron density in bound states with po- 
lar molecules is asymmetric, with a strong pile-up 
outside the ne gatively-charged e nd of the molecule 
(IBuenker et all l2007t IStrasburgerL [l999l l200ll) . This 



is 



shown by Fig. [8] for LiH. This figure also shows that cor- 
relations (e.g., included through ECG) produce a large 
increase in the positron density at the molecule, com- 
pared to the static HF calculation. This increase in the 
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FIG. 8 Positron density in e + LiH along the molecular axis: 
dashed curve, HF (et = 0.16 eV); solid curve, ECG (eb = 
0.94 e V) ; Li atom is at X = —1.51 a.u., and H at X — 1.51 a.u. 
From IStrasburgeil(|l999l ). 



Both the electron and positron densities change 
little with vibrationa l excitation of the molecule 
([Gianturco et aLl . l2006l) . In the alkali hydrides, positron 
binding noticeably increases the bond lengths and soft- 
ens the vibrational modes. For example, the energies 
of the first vibrational excitations in LiH and e + LiH 
are 168 and 109 meV, respectively ()Gianturco et all 
l2006t iMella et all |2000|) . This is related to the fact 
that the structure of this molecule is closer to Li + HPs 
than to a loosely bound positron orbiting the neutral 
molecule. In contrast, in BeO and MgO, positron binding 
changes the bond leng ths and the vibrational fr e quen- 
cies by less than 1% ( Buenker and Liebermannl . 120081: 
IBuenker et~al\ . 120071 ). 

At present there are no ab initio calculations of 
positron binding to alkanes or other nonpolar or weakly 
polar molecules (except for C2o). Positron binding in 
such systems is exclusively due to electron-positron cor- 
relation effects. To obtain a bound state, the calculation 
must include them accurately, since binding does not ex- 
ist at the static (e.g., HF) level. As an alternative, and 
give n the availability o f experimenta l data (Sees. Ivl and 
EU), iGribakin and Led (|2006aL l2009h explored positron- 
molecule binding to alkanes using ZRP. The ZRP is the 
simplest form of a mo del potential, suited to study- 
ing low-energy processes ( Demkov and OstrovskvL [l988h . 
The idea of the model potential approach is to fit the po- 
tential to experimental data (e.g., the binding energy for 
a given molecule) and then use it to study binding for a 
range of similar molecules. 

In the ZRP method, the bound-state wave function of 
the positr on in the field of N centers p laced at Ri has 
the form ( Demkov and OstrovskvL fl988l ). 



*(r) 



A g-Klr-R.,1 

l—l ' 



(51) 
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where n > is related to the bound-state energy by £o = 
—k 2 /2. The interaction with each center is parameterized 
by KQi through the boundary condition, 



const x (|r — RJ 



«0i)- 



(52) 



Subjecting ^ from Eq. (j5TI) to N conditions (|52|) yields a 
set of linear homogeneous equations for A», whose solv- 
ability determines the allowed values of «. 

The alkanes, C n H2„+2, were modeled by a planar zig- 
zag chain of n ZRPs, each representing the CH3 or CH2 
group. The distance between the neighboring ZRPs is 
given by the length of the C-C bond 2.91 a.u., and the 
angle between adjacent bonds is equal to 113°. The pa- 
rameter KQi = —0.69 a.u. was chosen to reproduce the 
binding energy for dodecane (n = 12, e — —220 meV) 
(|Gribakin and Leei 120011 . Figure [5] compares the results 
of this calculation with the measured binding energies for 
alkanes up to n — 16 (cf. Table [Vj. 



close to the center of the molecule. On the density plot 
(Fig.[TUl right), this corresponds to an area of low density 
near the midpoint. The actual positron wave functions 
are expected to differ from that given by Eq. (|51l) in that 
the latter does not exclude the positron from the regions 
inside the atomic cores (which are of "zero range" in the 
model). However, the atomic cores are relatively small 
compared to the extent of the positron wave function, 
in keeping with the main assumption of the ZRP model. 
Thus the model captures the main features of positron- 
molecule bound states. 
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FIG. 9 Positron binding energies for alkanes from experiment 
(crosses) (Table W\ and ZRP model calculation (circles). The 
ZRP model is fit to Eb = 220 meV for dodecane. 

Figure |9] shows that the model gives a good over- 
all description of positron binding to alkanes, including 
the prediction of the second bound state. However, it 
fails to capture in a quantitative way some of the de- 
tails. The model predicts binding for n > 4, whereas 
ethane (n = 2) is observed to bind positrons. The 
model also predicts that a second bound state emerges 
for n — 1 3, while experimentally this state is observed 
at n = 12 ({Barnes et all 120061 : [Young and Surkol l2008bl ) 
(cf. Fig.ES). 

To visualize the bound states, the two-dimensional 
density 



\V(x,y,z)\ 2 dz, 



(53) 



where x and y are in the plane of the carbon chain, is 
shown in Fig. [TU]for the first and second bound states of 
tetradecane (n — 14). Both states are quite diffuse, with 
the positron spread over the whole molecule. The wave 
function of the second bound state must be orthogonal to 
the ground state, and so it changes sign on a nodal surface 



FIG. 10 Two-dimensional density of the positron wave func- 
tions in the ZRP model for the 1st (left) and 2nd (right) bound 
states in tetradecane (n = 14). In these plots, the first car- 
bon atom is at the origin, and the C-C bonds are alternately 
parallel and at 67° degrees to the x axis. 



III. EXPERIMENTAL TOOLS AND PROCEDURES 

In this section, experimental techniques are described 
that have been used to investigate low-energy positron- 
molecule annihilation and related effects. The method- 
ology has varied considerably over the past half century 
of these studies, and no attempt is made to be complete. 
Emphasis is placed upon typical and/or best practices. 
The reader is referred to the original papers for further 
details. 



A. Annihilation-rate measurements with thermalized 
positrons in atmospheric pressure gases 

Positrons from conventional sources such as ra- 
dioisotopes or electron accelerators typically have 
energies ranging from ~ 1 keV to ~ 0.5 MeV. Thus, 
to study positron interactions at low energies, some 
method must be used to slow the positrons. Early 
annihilation rate measurements were done using the 
test species themselves as "moderators" and measur- 



production and annihilation signals 


(Dcutsch, Il951allb 


Griffith and Hcvlandl. 1978; 


Hevland et all 


1982 


Osmon. 1965a; Shearer and Deutsda.ll 949c 1 Wright et al. 



Il983l ). In a typical experiment, the 1.28 MeV gamma 
ray that accompanies positron emission from a 22 Na 
positron source provides a start signal, and the detection 
of a single, 511 keV gamma ray from a two-gamma 
annihilation event is used as a stop signal. Achieving 
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an adequate signal-to-noise ratio requires working at gas 
densities > 0.1 amagat. 

Annihilation rates are obtained from the dependence 
of the exponential time decay of the 511 keV annihilation 
signal (i.e., after the positrons come to thermal equilib- 
rium) on the test-gas density. Using fast electronics, it 
was also possi ble to measure the s l owing down of the 
fast positrons ( Sharma and McNuttl . 1 19781) . Where nec- 
essary, the rate of thermalization was increased by adding 
a light spe cies with a small value such as molecular 
hydrogen (jWright et all . Il985t l , and the long-time ortho- 
positronium component was quenched using a small ad- 
mixtur e of a gas with u npaired electron spins such as NO 
or 2 (|Deutschl . [l95laM rJ). 



B. Buffer-gas positron traps as tailored sources of positrons 



Buffer-gas (BG] trap s ([Murphy and Surkol . Il992t 
ISurko et all . Il999l Il988bl ) proved to be a useful tool to 
tailor positron gas es, plasmas and beams for positron an- 
nihilation studies ( Surko and Greavesl I2004T ) . Positrons 
from a sealed 22 Na radioactive source are slowed to 
electron-volt energies usi ng a solid-neon refl e ction mod- 
erator [efficiency 1-2% ([Greaves and Surkol Il996l) ]. A 
50 mCi 22 Na source and neon moderator produces ~ 
5—10 x 10 6 slow positrons per second. The slow positrons 
are then guided magnetically into a b uffer- gas, Penning- 
Malmberg trap ( Murphy and Surkol . Il992l ). It consists 
of a uniform magnetic field (B ~ 0.1 T) coaxial with 
a set of cylindrical electrodes biased to form a stepped, 
three-stage, potential well. The stages contain a nitrogen 
buffer gas with successively lower pressures. Positrons 
become trapped by losing energy through electronic ex- 
citations of the N2 molecules. They then cool to the 
ambient (i.e., room) temperature in the third stage by 
additional collisions with the N2. More recently, a small 
amount of CF4 was added to increase the coolin g rate 
dGreaves and Surkol 1200(1 ISullivan efal\ . l2002bTl . The 
result is a magnetized, thermal gas (or plasma) of as 
many as 10 8 positrons at a temperature of 300 K. The 
trapped positrons are in a background gas pressure < 
10~ 6 torr, which can then be pumped out, depending 
upon the experiment. 



C. Annihilation-rate measurements in positron traps 

Positron traps have enabled an improved method 
to study the interaction of t hermal posi t rons with a 
large variety of te s t spe c ies (llwata et al , 19951 120001 ; 
iMurphv and Surkol Il99ll ISurko et all . Il988bf ). In this 



case, the positrons are trapped, the buffer gas is pumped 
out, and a test gas species is introduced into the trap 
at a low pressure. Low pressures ensure that the an- 
nihilation events are due to two-body positron-molecule 
interactions and thus that three-body processes are neg- 
ligible. This technique also permitted study of low vapor 



pressure targets. Furthermore, the positron temperature 
(e.g., 300 K) could be measured to verify that they are in 
thermal equilibrium with the test species. This was done 
using standard plasma techniques, dumping the positrons 
and measuring thei r energy distribution using a retarding 
potential analyzer ( Murphy and Surkol 1 19921 ). 

The apparatus for these therm al annihilation-rat e 
measurements is shown in Fig. [TT] (llwata et all , Il995f ). 
In this case, the trapped positrons were shuttled to a 
separate confinement stage surrounded by a cryogenically 
cooled surface to reduce the level of impurities in the vac- 
uum system. It operated with liquid nitrogen (77 K) or 
an ethanol-water mixture (— 7°C) depending upon the 
test species. The annihilation was monitored by hold- 
ing the positrons for a given time, then dumping those 
that had not annihilated onto a collector plate and mea- 
suring the gamma ray signal using a Nal(Tl) scintilla- 
tor and a photomultiplier. Experiments were also con- 
ducted in which the thermal positrons were heated by 
short bursts of radio-frequency noise to provide a mea- 
sure of the a nnihilation r ate as a function of positron 
temperature ( Iwata et all . 12000). 
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FIG. 11 Schematic diagram of the apparatus used to study 
the interaction of trapped, thermal positrons with low- 
pressure gases. The positrons are contained in a fourth trap- 
ping stage surrounded by a cold surface to minimi ze the effects 
of imp urities present in the vacuum system. From lTwata et all 
(jl995h . 

For test species that are gases at 300 K, the pressure 
is controlled to micro-torr precision by feedback using a 
capacitance manometer and a piezoelectric valve. A liq- 
uid delivery system, while offering less control, permitted 
studies of a broader range of test species. It can be used 
with species that have relatively low boiling points (i.e., 
close to room temper ature) and pu lverized solid samples, 
such as naphthalene ( Youngl . [20071 ). For liquid samples, a 
freeze-pump-thaw procedure is used to eliminate volatile 
contaminants. The sample is placed in a temperature- 
regulated bath, and a needle valve is used to leak va- 
por into the annihilation cell. While there was no feed- 
back control of the pressure, it was still reasonably stable 
(|lwata et aLl . ll995l ). 

Some vapors must be run at pressures below the sen- 
sitivity of the capacitance manometer to avoid detector 



22 



saturation. In this case, the pressure of the test species 
is measured with an ion gauge calibrated against the 
manometer at higher pressures. A few species studied 
(e.g., naphthalene) have such low vapor pressu res that 
accurate pressure calib r ation w as not possible (Young], 
120071 : lYoung and SurkoL feoOSbl) . 

The major uncertainty in these Z e g measurements is 
estimated to be a ±20% systematic error in the measure- 
ment of test-gas pressures. 



D. Trap-based cold positron beams 

To obtain detailed information about positron inter- 
actions with matter, the tool of choice is a beam of 
positrons with a well defined and tunable energy. In early 
experiments, the beam energy resolution was limited by 
the energy spread of positrons emerging from a modera- 
tor (e.g., a fraction of an e lectron volt). While the use of 
a cold primary moderator (iBrown et all [l986) or remod- 
erator ( Gullikson and MillsT 1987|) could reduce the en- 
ergy spread to 30-40 meV, these techniques were rarely, 
if ever used in atomic physics experiments. The advent of 
the BG positron accumulator enabled an efficient, pulsed 
and tunable, low-energy positron beam wi th a compara- 
ble, and potentially smaller ener gy spr ead ( Gilbert) . l2000t 
iGilbert et al. |, Il997t iKurz et aZ.Ul998| ) . 

In this technique, ~ 10 4 positrons are ac cumulated in a 
time ~ 0.1 s and cooled fo r a similar time ( Gilbert et all 
l2002HSullivan et aLl . l200"i . Then the exit-gate electrode 
is lowered to a potential Ve , which sets the beam energy. 
The bottom of confining potential well is then raised 
(i.e., to ~ 0.25 V above Ve) in a time ~ 5 /is to pro- 
duce a positron pulse with a similar time duration. The 
beam-transport energy (typically ~ 2 eV) is set below 
the threshold for positronium formation to avoid positron 
loss and extraneous gamma-ray signals. All this is done 
in the BG-trap magnetic field (i.e., B ~ 0.15 T), produc- 
ing pulses of 1—3 x 10 4 positrons at a few Hz rate. The 
pulses are magnetically guided to the annihilation cell in 
fields ranging from 0.03 to 0.1 T. 



The beam energy distribution is measured using the 
cylindrical electrode o f the gas cell as a retarding poten- 
tial analyzer (RPA) (jYoungl . 120071 ). The mean energy 
of the beam in the annihilation cell can be verified in- 
dependently by measuring the time-of-flight delay of the 
positrons passing through the cel l as a function of in- 
creasin g cell potential Vc ( Marleii 120051: ISullivan et all 
2002b). Typical uncertainties in the mean beam energy 
40 meV. 



are 



The parallel e nergy distribu tion can be modeled by a 
single Gaussian ( YouneJ . [2007|) . At a more accurate level, 
there is typically a high-energy tail containing ~ 10% of 
the beam particles, depending upon the beam-formation 
protocol and the relative magnitudes of the magnetic 



field in the buffer- gas trap and the measurement region. 19 
Parallel energy spreads as small as 18 meV (FWHM) ca n 
be achieved using this technique ( Gilbert et all Il997th 
when the beam is kept in the same strength magnetic 
field as the buffer-gas positron accumulator. Typical val- 
ues of parallel energy spread in the annihilation experi- 
ments described here are somewhat larger, 25—30 meV, 
due to the fact that the annihilatio n cell was at a lowe r 
magnetic field than that of the trap ( Barnes et a/.l . l200"3h . 



The spread in transverse energies in the BG trap is set 
by the ambient temperature Tj_ = 25 meV (i.e., corre- 
sponding to an electrode temperature of ~ 300 K). 

The positron energy distribution in the beam is mod- 
elled as 
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(54) 

where en and e± and the parallel and transverse positron 
energies, a is the root-mean-squared width of the paral- 
lel energy distribution (i.e., corresponding to a FWHM 
of oV81n2), and e is the mean parallel energy of the 
positron beam. 

This distribution has been verified using the resonant 
annihilation peak for the C-H stretch modes in propane 
(cf. Sec. [V}. The line width of the VFR peak is assumed 
to be negligible (i.e., < 1 meV; cf. Sec. HH, as is the 
spread of the propane C-H stretch frequencies. The dis- 
tribution in Eq. (|54[) is convolved with a delta function, 
yielding t he distribution of ( total) positron energies of 
the beam (Gribaki n and Led . l2006bl ). 
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where $(x) is the error function. Note that this energy 
distribution depends on the difference between the total 
energy e and mean parallel energy e. 

The resulting fit, shown in Fig. [12] is in excellent agree- 
ment with the data. Due to the spread in perpendicular 
energies, the position of the peak observed as a function 
of the mean parallel energy (which is set and measured 
by varying Vc) is about 12 meV below the true energy of 
the peak (i.e., as a function of the total positron energy). 
The annihilation spectra in this review are presented as 
a function of the total positron energy, which is taken to 



19 The ratio e±/B and the total positron energy e = Sm + e± arc 
both conserved in this magnetic beam-transport system, where 
eii is the parallel beam energy, and e±_ is the transverse energy 
of the particles' gyromotion in the plane perpendicular to the 
magnetic field B fearnes ef "oZI . 12003 ; fSullivan et n?ll2002bt l. 
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be 12 meV higher than the mean parallel energy set in 
the experiment. 
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FIG. 12 Measurement of the distribution of total positron 
energies in the trap-based beam using a vibrational- Feshbach 
annihilation resonance: solid curve, prediction of Eq. (|56l) . 
normalized arbitrarily, shifted in energy, and fitted to the 
energy-reversed, normalized, C-H stretch peak in Z e s for 
propane (•), with T± — 26 meV and a pa rallel energy spread 
(FWHM) of 27 meV. From lYounel (|2007l h 



E. Energy-resolved annihilation measurements 

Central to this review are measurements of positron 
annihilation on molecules, resolved as a function of 
incident pos i tron energy dBarnes et al. 



„J2003J, 120061 

iGilbert et al\ . 120021: lYoung and Surkol. l2008bl h The anni- 



hilation cell is shown schematically in Fig. [13] It consists 
of a cylindrical, gold-plated electrode 4.4 cm in diameter 
and 17 cm long. The gamma-ray detector and associated 
shielding restrict the detector field-of-view to a region 
< 15 cm in length along the axis of the cell. Magnet 
coils outside the cell impose a field of ~ 0.075-0.095 T, 
with the lowest value in region viewed by the detector. 
Metal baffles shield the gas cell and the detector from 
spurious gamma-ray decays. The system for handling 
gases and vapors is described in Sec. IIII.CI above. 
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FIG. 13 Schematic diagram of the gas cell, shielding, and 
detection apparatus used for energy-resol ved posi t ron a nnihi- 
lation measurements (not to scale). From] Young (20(53). 

Single gamma rays are detected using a Csl crystal 



and a photodiode, followed by a single-channel analyzer 
centered on the 511 keV annihilation gamma-ray line. 
The absolute detector efficiency and the sensitivity along 
the path of the positron beam are measured using a cal- 
ibrated gamma-ray test source ( Gilbert et all [2003) . 

Pulses of positrons pass through the gas cell several 
times while annihilation events are recorded, with the to- 
tal scattering kept below 15%. A typical time window fo r 
this measurement is ~ 15 fis (lYoung and Surkol 12008b). 
Positrons are kept in flight while the annihilation events 
are recorded to avoid spurious gamma-ray background 
signals. To avoid detector saturation, the average sig- 
nal level is adjusted to ~ 1 count per 10 positron pulses. 
Typical test-gas pressures range from 0.1 to 100 /xtorr, 
due to the large variations in annihilation rate for dif- 
ferent chemical species. Background signals can be as 
low as one count per 10 9 positrons cycled through the 
annihilation cell. 

A typical spectrum consists of ~ 10-25 pulse-train 
measurements at each energy, taken at 10-15 meV in- 
tervals over the relevant range of positron energies (e.g., 
< 500 meV). This measurement is then repeated a few 
hundred times. Complications due to scattering restrict 
measurements to energies > 50 meV. Absolute values 
of Z e ff are obtained from measurement of the detector 
efficiency, the pulse strength, the detector-sensitivity- 
averaged path length, and the test-gas pressure. 

Uncertainties in these parameters are estimated to re- 
sult in a 20% overall systematic uncertainty in the ab- 
solute magnitude of Z g. The error bars shown in this 
review indicate the statistical uncertainty due to the fi- 
nite number of counts at a given energy. In many cases, 
these statistical errors are smaller than the size of the 
data points. The linearity of the signal with test-gas 
pressure is checked to ensure that scattering and three- 
body effects are negligible. 

A separate annihilation cell was used for studies of an- 
nihilatio n as a function of the te mperature of the test 
species ( Young and Surkol . [^OOSal . It is capable of cool- 
ing target gases down to 100 K. To ensure that the test 
species are actually thermalized, the annihilation cell was 
operated at temperatures below the sticking temperature 
of the test species on the cell walls. 



F. Gamma-ray spectral measurements 

The center-of-mass momentum of the annihilating 
electron-positron pair contains important information 
about the momentum distribution of the bound elec- 
trons in an atomic or molecular target. It produces both 
Doppler shifts of the photon energies and changes in their 
directions (see Sec. III.B[> . Using such measurements, one 
can distinguish, for example, annihilation on different 
atoms, or annihilation on valence and inne r-shell elec- 
trons in a given species ( Iwata et aZ.I . Il997bh . 

In a two-gamma event, the photon momentum distri- 
bution can be determined by measuring either the an- 
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gular deviation of the gamma rays [the so-called angular 
correlation of annihila tion radiation (ACAR) technique 
( Coleman et aZ.I . fl994ft ] . or the Doppler broadening of the 
annihilation gamma-ray line. For the tenuous samples 
studied here, Doppler broadening is the method of choice. 

The a pparatus for these measurements is shown in 
Fig. [□] ((iwata et all Il997al) . The gamma rays are de- 
tected using a high-purity Ge detector with an energy 
resolution of 1.16 keV (FWHM). Positrons are trapped 
and cooled, then the N2 trapping gas is pumped out, 
and the test gas is introduced. Gamma-ray data are 
recorded, and then this cycle is repeated. The errors in 
the measurement of gamma-ray energies are estimated to 
be ±0.02 keV. The errors in spectral intensity are pre- 
dominantly statistic al due to the finite number of counts 
(llwata et all Il997al ). 



G. Annihilation-induced fragmentation 

Positron annihilation on atomic and molecular species 
produces a spectrum of positive ions that can be mea- 
sured using time-of-flight techniques. The first experi- 
ment of this kind arranged for positrons in a Penning trap 
to interact with molecular species for a short time. Then 
the resulting ions were dumped and the m ass spectrum 
measured using a t i me-of -fiight technique ( Glish et~al\ , 
I1994 IPassner et all . ll989l ). Subsequently, more detailed 
studies were conducted using i mproved techniques to 



measure the ion mass spectra (iDonohue et al , 199C : 
Hulett et artll99lll996aHMoxom et aZ.L i2000; Xu et al. 



1994 , 19931 . Il995l ). In this case, positrons from an elec 



tron LINAC were moderated and then accumulated in a 
Penning trap, where they were allowed to interact with 
the test species. The product ions were accelerated and 
detected using a microchannel plate. Significantly higher 
mass resolution was achieved using a spatially varying 
(quadratic) potential to arrange the same arrival time 
for ions of the same mass starting at different initial po- 
sitions in the trap. 



IV. ANNIHILATION ON SMALL MOLECULES 

Research on positron annihilation on atomic and 
molecular targets has typically been focused in two ar- 
eas: atomic and small molecular species where mod- 
est annihilation rates are observed, and large hydrocar- 
bons that are characterized by very large annihilation 
rates. Som e small polyatomic molec ules occupy a mid- 
dle ground ( Gribakin and Leel , l2006bl ) . Here annihilation 
proceeds in a different manner than in other small targets 
due to the fact that these molecules can temporarily bind 
a positron via VFRs. Further, in contrast to larger hy- 
drocarbons, the theoretical description of annihilation in 
these targets is relatively simple. The focus of this section 
is VFR-mediated annihilation in small molecules and the 
contrasting case of annihilation in small molecules that 



do not bind positrons. 

Energy-resolved Z c g spectra, such as those shown in 
Fig. [TJ reveal peaks corresponding to the resonant trans- 
fer of energy from the positron to specific molecular vi- 
brations. In each of these resonances the positron be- 
comes temporarily attached to the molecule, resulting in 
a greatly enhanced annihilation rate. The energies of the 
resonances are given by Eq. ([3]), namely, e v = E v — 
where E v is the energy of the excited vibration v and 
£b is the positron-molecule binding energy. Fundamental 
vibrations for which E v — uj v , can produce strong annihi- 
lation. As a result, the molecular annihilation spectra are 
somewhat similar to infrared-absorption spectra. How- 
ever, the magnitudes of the annihilation resonances are 
not proportional to the IR absorption strengths. They 
follow a different scaling with molecular size. While IR- 
active modes dominate the spectra of many molecules, 
there are cases in which nominally IR-inactive modes, as 
well as combinations and overtones, also appear to pro- 
duce annihilation resonances. 



A. Halomethanes as a benchmark example of VFR 

The singly halogenated methanes are a near-perfect 
set of molecules to test theories of positron VFRs. Each 
molecule has only six vibrational degrees of freedom, 
all of which are dipole-active. There are three energy- 
separated pairs of fundamental vibrations: the C-H 
stretch modes, the C-H bend modes, and the C-X modes 
where X is the halogen. The energy-resolved annihi- 
lation spectra and the infrared absorption spectra for 
these molecules are shown in Fig. [TU In CH3CI and 
CH3B1', one can discern VFR from all of the infrared- 
active modes. The high-energy peak is due to the C-H 
stretch mode, and the broad low-energy feature is due to 
the C-H bend and C-X modes. There is no evidence of 
multimode VFRs. 

As shown in Fig. [T3] the position of the C-H stretch 
resonance shifts downward in energy as the size of the 
halogen atom is increased, reflecting an increase in the 
positron binding energy. The binding energies range from 
near zero in CH3F to about 40 meV in CHaBr. Since vi- 
brational Feshbach resonances cannot occur unless the 
positron is bound to the molecule, the small positive en- 
ergy shift of the C-H stretch peak in CH3F is likely the 
result of a very small binding energy and a small positive 
shift in the mode energy. 

According to Eq. (|4"2"|) , the annihilation rate in small 
molecules can be described by a sum of Breit-Wigner 
resonances for each mode, v, convolved with the instru- 
mental positron energy resolution function /b(£). All of 
the modes are dipole active in the halomethanes. Thus 
the elastic capture rates, which are roughly proportional 
to the IR strengths, are expected to be much larger than 
the annihilation rates (see Sec. III. Fl) . As a result, T ~ T e , 
and Eq. (|4"2"|) is greatly simplified. Using Eq. (|3"6"|) . one 
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FIG. 14 Annihilation rates Z e ff (•) and IR absorption (solid 
curves) for methyl halides: (a) CH3F, (b) CH3CI, and (c) 
CHaBr ijBarnes et all 120061 '). Vertical bars below each plot 
indicate the vibrational energies from IR and Raman mea- 
surements, and selected mode frequencies ( "all modes" ) from 
(jLinstrom and Mallard l2005h . 
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FIG. 15 Comparsion of experimental and theoretical Z c a for 
methyl halides CH3X (• and solid curves) and CD3X (o and 
dashed curves) for (a) X=F, Eb = 0.3 m eV; (b) -Y=C 1 , £b = 
25 meV; and (c) ^=Br e b = 40 meV dBames et all |2006| ; 
iGribakin and Led . [2006bl ; lYoung et al. I l2008h . Dotted curves 
show the contributions of direct annihilation. 
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where g v = ye^/ev- Since F is assumed to be constant 
(see Sec. III.Dj) . the relative magnitudes of the dipolc- 
active resonances in a given small molecule are deter- 
mined by the g v factors. The only adjustable parameter 
is the binding energy, which can be determined by com- 
parison with experiment. The application of Eq. (|57j) 
to the halomethanes and their deuterated counterparts 
is shown in Fig. 1151 The small contribution from non- 
resonant direct annihilation, described by Eq. (|2T)1) . is 
included. 

For the two larger halogens, the binding energy was ob- 
tained from the position of the C-H stretch peaks, while 
for CH3F, the (small) binding ener gy was determined 
by fitt ing to the magnitude of Z e g (jGribakin and Led . 
l2006bl ). The resulting binding energies for CH 3 F, CH 3 C1, 
and CHaBr are 0.3, 25, and 40 meV, respectively, and the 
agreement between the theoretical and experimental Z c r 
spectra is remarkably good. As the binding increases, 
each annihilation peak shifts downward in energy and 



increases in magnitude in accord with the predicted g- 
factor scaling. The theory also predicts successfully the 
absolute magnitudes of each of these features. 

For a wide variety of hydrocarbons, experiments have 
established that the positron binding to the molecule is 
unchanged when the hydrogen atoms are replaced with 
deuterium atoms. This is plausible, since the binding 
energy is expected to be primarily a function of the 
electr onic rather than th e vibrational, degrees of free- 
dom ( Barnes et all , 120031 ). The only difference is that 
the deuterated species have lower vibrational energies 
due to the larger reduced masses associated with the C- 
D modes. The binding energies for the hydrogenated 
species were used to predict the Z g spectra for the deu- 
trated species, thus providing a stringent test of the the- 
ory. The results of this comparison are shown in Fig. 1151 
Note that the g-factors and hence the magnitudes of 
the peaks, are larger for the deuterated species, because 
the resonances occur at smaller positron impact energies. 
The theory, now with no free parameters, works well for 
the deuterated halomethanes just as it did for the hydro- 
genated species. To our knowledge, no other theory has 
demonstrated such close agreement, in both magnitude 
and shape, with the observed positron-molecule annihi- 
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lation resonances. 



B. Methanol: a case of multimode VFR 

In methyl halides, each vibrational mode produces a 
measurable VFR with a relative magnitude given by the 
g factor, and there is no evidence of multimode excita- 
tions. Experiments show that this "selection rule" must 
be relaxed for a variety of other molecular species. As 
in electron-molecule collisions, positrons can be expected 
to excite vibrations that are nominally dipole forbidden. 
The nature of annihilation VFRs, as described by the 
theory, will be to even out fairly large variations in the 
capture rate as long as T e v T a . As shown in Fig. [TBI 
there is evidence of such higher-order vibrations in the 
Z e ff spectrum of methanol CH3OH. This molecule is iso- 
electronic with CH3F. Its vibrational spectrum is also 
similar except for the additional O-H stretch vibration. 
However, the Z c g spectrum of methanol is quite different 
than that of CH3F and the other methyl halides. There 
is a significant increase in magnitude of the high- and 
low-energy peaks in methanol relative to those of CH3F. 
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FIG. 16 Comparsion of experimental Z e g (•) for methanol 
(CH3OH) with theory: dotted curve, contribution of direct 
annihilation '; dashed curve, total Z e g due to VFR of 
vibrational fundamentals for Eb = 2 meV; solid curve, total 
Z e ft due to resonant annihilation inv olving the 12 modes and 
9 overtones and combinations. See lYoung et al\ (|2008l ) for 
details. 



The positions of the C-H stretch peak and the peaks 
at lower energies indicate that the binding energy in 
methanol is small. The fitted binding energy of CH3F 
is 0.3 meV; in methanol it could be an order of magni- 
tude larger but still remain within the experimental en- 
ergy uncertainty. There is evidence of an additional peak 
above the C-H stretch modes, presumably due to the O- 
H stretch mode which has the energy woh = 456 meV. 
If this interpretation is correct, this VFR is downshifted 
relative to o>oh by an amount somewhat greater than the 
positron binding energy. 

While one can vary the binding energy to describe bet- 
ter some of the enhancement in Z e g seen in Fig. 1161 there 



are features in the spectrum that cannot be explained by 
mode-based VFRs. In particular, the higher-energy peak 
in methanol is much broader than the corresponding C- 
H stretch peaks in the halomethanes and the prediction 
based on the IR-active modes. It is also much closer in 
magnitude to the low-energy peak. This discrepancy can 
only be resolved by considering additional resonances. 

Figure [TB] shows the results of two calculations using 
a binding energy of 2 meV. The dashed curve is the 
calculation which includes only the fundamental vibra- 
tions, all of which are dipolemactive. It falls significantly 
short of explaining the spectrum. The IR absorption 
measurements in methanol reveal a number of relatively 
weak overtones and combination vib rations (likely som e 
of which are Fermi resonances) ( Bertie and Zhanglll997h . 
Using these data allows one to estimate the elastic rates 
Y e v . They are generally smaller than those of the funda- 
me ntals, but still satisf y the relation T e u T a [see Table I 
by (jYoung et all I2008T )] . The result of adding these nine 
two-quantum overtones and combinations to the twelve 
fundamentals is shown by the solid curve in Fig. 1161 This 
prediction is clearly in better agreement with experiment. 
This comparison provides strong evidence that multi- 
mode vibrations can make significant contributions to the 
Z e g spectra of small molecules. There are two remaining 
discrepancies. One is the significantly downshifted O-H 
stretch peak. The second is the higher experimental Z e s 
values below 100 meV. In methanol there is a torsion 
mode at ^40 meV. Its overtones and combinations could 
provide the missing spectral weight in this region, but 
there is at present no estimate of this effect. 



C. VFR from dipole-forbidden vibrations 

The theoretical analysis of methanol indicates that 
multimode vibrations can contribute significant spectral 
weight to an annihilation spectrum. However, there is 
also strong evidence that modes with very weak (or nomi- 
nally zero) dipole coupling can also produce VFRs. They 
could, for example, arise from higher-order nondipole 
coupling (e.g., electric-quadrupole-active modes), but 
there is no simple way to assess their possible contri- 
butions. These nondipole features are exemplified by 
the experimental and theoretical Z g spectra of ethylene 
shown in Fig. [T7] (a). This molecule has five IR-active 
modes and six IR-inactive modes. The shift of the C-H 
stretch peak indicates a binding energy of about 10 meV. 

Including only IR-active {B u symmetry) modes under- 
estimates Z e H by a factor of 2 at the peaks located at 
100 and 350 m eV. This calculation , which uses the IR 
strengths from iBishop and Cheung! ( 1982t ). completely 
misses contributions in the interval of energies between 
these modes. However, the form of Eq. (|4"2")l does not de- 
pend on the nature of positron- vibrational coupling, and 
all vibrations with sa T e u 3> T a contribute equally to 
Z e g (i.e., to within the g factor), as per Eq. (fST)) . As 
shown in Fig. 1171 including the remaining modes (i.e., 
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those with A g and B g symmetry, corresponding to the 
capture of s-, p-, and <i-wave positrons) greatly improves 
the agreement. The further inclusion of the fourteen IR- 
active combination vibrations (not shown) results in Z c s 
values exceeding those observed experimentally. How- 
ever, as shown in Fig. I17[ if one weighs their contribu- 
tions by an empirical factor T1/T u = 1/n, where n is 
the number of vibrational quanta involved, the predic- 
tion agrees well with the experimentally measured spec- 
trum. This analysis shows that nearly all of the spectral 
weight between the two largest peaks in ethylene can, at 
least in principle, be attributed to multimode VFR. How- 
ever, there is presently only an ad hoc approach available 
to decide whether "borderline" vibrational capture chan- 
nels sit above or below the cutoff value of the coupling 
strength set by T a . 



trum matches the magnitude of that observed. Still, the 
peaks in the theoretical Z e g spectrum are much more 
prominent . A similar an alysis has been performed for 
ammonia (|Gribakinl[20T0l) . It also highlights the role of 
overtones and combinations, and suggests that rotational 
broadening of VFR and rotational Feshbach resonances 
may be required to explain the Z c g spectrum in this small 
molecule. 

The modeling discussed here indicates that, to explain 
the Z e g spectra of molecules such as C2H4 and C2H2, 
one needs to include the VFR of IR-inactive modes, over- 
tones, and combinations. However, determining Y e v and 
r„ for these resonances requires calculations of the full 
vibrational dynamics, and so at present the theory is in- 
complete. 
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FIG. 17 Comparsion of the experimental Z e g (•) with the- 
ory for (a) ethylene (C2H4, ej, = 10 meV), and (b) acety- 
lene (C2H2, Eb = 5 meV). Dotted curves, direct annihilation 
Z^g 11 ' '; solid curves, total Z e s due to all IR active modes. In 
(a) dashed curve, total Z e ff due to A g , B g , and B u modes 
that capture s-, p-, and d-wave positrons; chain curve, same 
with the addition of 14 IR-act ive overtones and combinations 
listed in (|Georges et ai, 1999). In (b) long-dashed curve, all 
modes; chain curve, all modes with E 9 , u , II 9lU , and A 9 sym- 
metries. The contributions of all overtones and combinations 
are ad hoc weighted by a factor 1/n, where n is the number 
of vibrational quanta. 



Shown in Fig. [IT] (b) is the spectrum of acetylene 
(C2H2). Its binding energy is too small to measure di- 
rectly, and hence it is obtained by fitting to the spectrum. 
As in the case of ethylene, the theoretical Z c g which in- 
cludes only IR-active fundamentals, does not agree with 
the measured spectrum. Adding all of the other modes 
provides little improvement. It is only after adding the 
overtones and combinations that the calculated Z c g spec- 



D. Effect of molecular size on the magnitudes of VFR 

As discussed in the next section, larger molecules ex- 
hibit annihilation VFRs that cannot be explained by 
positron coupling to fundamentals or combination and 
overtone vibrations. In these molecules, an additional 
enhancement mechanism appears to be operative that 
causes the magnitudes of the fundamental resonances to 
grow rapidly with molecule size. The smallest molecule 
to show evidence of such enhanced VFR is ethane. As 
shown in Fig. 1181 the Z c r spectrum for this molecule has 
a distinctly different spectral shape than that of the other 
molecules discussed in this chapter. The high-energy C- 
H stretch peak is three times larger than the low-energy 
C-H bend peaks, which is inconsistent with the simple 
g scaling of VFR magnitudes. The calculated Z e g for 
ethane, shown in Fig. [T5J indicates that the VFRs of 
IR-active and dipole-forbidden modes, populated by the 
positron s-, p-, and d-wave capture, can explain the spec- 
tral weight at lower energies. However, they do not ac- 
count for the magnitude of the C-H stretch peak, nor for 
the magnitude of Z c g between the low- and high-energy 
peaks, which may be due to overtones and combinations. 

Similar analyses for propane and cyclopropane (see 
Sec. lV.F[) show that these molecules exhibit even stronger 
enhancements of the high- and low-energy peaks, while 
the "background" of multimode excitations between 
these resonances does not appear to be similarly en- 
hanced. However, molecular size alone does not appear to 
be an accurate predictor of the scaling of Z c g magnitudes. 
As shown in Fig. [T51 the analysis for ethanol (C2H5OH; 
one more atom than ethane) indicates that there is no 
evidence of such an enhancement for this molecule. Here 
the agreement with experimental Z e g of the simple IR- 
active-mode VFR theory without overtones or combina- 
tions is remarkably good. 

In summary, the theoretical framework of 
iGribakin and Led {2006b) allows us to make a clearer 
distinction between those molecules that exhibit "small 
molecule" behavior and those that exhibit "large 
molecule" behavior (i.e., where Z e g values can be orders 
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FIG. 18 Experimental Z c g spectrum (•) and theory for 
ethane (C2H6). The calculations use et = 1 meV and in- 
clude: dotted curve, direct Z e s\ solid curve, same with VFR 
of IR-active modes; dashed curve, same with the addition of 
Ai g modes; chain curve, same with the addition of E g modes. 
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FIG. 19 Comparsion of experimental Z c g (•), IR absorption 
spectrum (solid curve), and theoretical Z e a (dashed curve) 
for ethanol (C2H5OH). The calculation [cf. Eq. (|42|) ] uses 
Eb = 45 meV an d the mode frequen cies and transition dipole 
amplitudes from IShaw et ali [j 19901 ). 



of magnitude larger than those predicted by the mode 
VFR theory). However, the physics that is responsible 
for this threshold is still poorly understood. 



E. Nonresonant annihilation in small molecules 

Many molecules do not exhibit resonant annihilation 
but display instead relatively smooth featureless Z c g 
spectra. To better understand these molecules, recall 
that the minimum requirements to observe annihilation- 
mediated VFR are the existence of a positron-molecule 
bound state and a vibrational mode that couples to the 
positron. 

The existence of weak binding with Eh ~ k 2 /2 is linked 
to a large positive value of the positron-molecule scat- 
tering length k^ 1 (Sec. III.DI) . The magnitudes of the 
resonances are then determined by g — njk. If k is nega- 
tive, the bound state is replaced by a virtual state in the 
continuum, and VFR are absent. In both cases, however, 
one expects a nonresonant background due to direct an- 



nihilation, proportional to (k 2 + n 2 )^ 1 [cf. Eq. (|29|) ]. 

Small nonpolar or weakly polar molecules are far 
less likely to bind positrons, and so it is not surpris- 
ing that many of these molecules lack VFRs. One 
such molecule is C02- It has a r elatively flat spec 
trum with Z e g « 35 above 150 meV (lYoung and Surko 
2008d) . and the thermal Z cS of 54.7 ([Wright et al 



19851 ). Since this molecule has 22 electrons, these val- 
ues are not far from the uncorrelated electron gas pre- 
diction. Using a vibrational cl ose coupling formalism, 
iGianturco and Mukheried (1999) predicted a resonance- 
free spectrum with a nearly constant Z c g ss 50, which is 
in reasonably good agreement with the measurements. 

The spectrum of methane CH4, shown in Fig. |2"TJI 
is also relatively featureless, and similar in both mag- 
nitude and e nergy dependence to the Z G ff for CF4 
(not shown) ([Barnes et al .], I2003D . This is consistent 
with these molecules not supporting the positron bound 
states. It refutes an earlier conjecture (jGribakinl , l2000t 
llwata et al. , 2000), which was based on the analysis 
of room-temperature Z e g for methane and its fluoro- 
substitutes. However, the difference between molecules 
close to the border between VFR-active and VFR- 
inactive species can be stark. Figure [20] shows a com- 
parison of CH4 and CH3F. The analysis presented in 
Fig. [T5l indicates that CH3F has a very small binding en- 
ergy (~ 0.3 meV). Yet this produces a distinct resonance 
at 150 meV, and the spectral weight in the vicinity of the 
C-H stretch mode is nearly doubled. 
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FIG. 20 Measured Z cS spectra for CH 4 (•) and CH3F (o), 
and IR absorption spectrum of methane (solid curve). 

As shown i n Fig. [21] the water mol ecule also lacks 
distinct VFR ([Young and Surkol . l2008d ). Its Z cS rises 
smoothly with decreasing energ y, reaching a value of 319 
for thermal positrons at 300 K ( Iwata et all Il995). The 
spectrum is well represented by the expression for direct 
annihilation [Eq. (|2"9"1) ]. 20 The fitting parameter k is con- 
sistent with a virtual state at ~ 1 meV. A constant Z c g 



20 In principle, Eq. ||29|I should be modified to account for the per- 
manent dipole moment of this molecule. 
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FIG. 21 Experimental Z c g spectrum (o) for water, and a fit 
based on the direct annihilation model: solid curve, Z e s from 
Eq. ((29j) with k 2 /2 = 0.3 meV plus a constant offset of 20. 



offset of 20 was also included in this fit [cf. the first term 
in Eq. (@8])]. 

The Z c g data for these small polyatomics can 
be com pared with the direct annihilation calcula- 
tions of IGianturco et al. I (120011) . In that work, the 
positron-molecule interaction is described by a local 
correlation-polarization potential, and Z c g includes the 
electron-positron contact den sity enhancement factor 



( Boronski and Nieminenl . il 9861) . The calculations predict 
a steady rise in Z e g with decreasing positron energy be- 
low 0.5-1 eV. The calculated annihilation rates for CH4, 
NH3, and H2O (including those at thermal positron en- 
ergies) are about a factor of 2 lower than the experimen- 
tal values, while for CF4, the calculated Z P ff is about 
two t imes larger. With the exception of NH3 ( Gribakinl . 
l2010i ). the calculations confirm that VFR are not neces- 
sary to explain the Z e g spectra in these molecules. 

Some molecules exhibit other spectral features. One 
example is the sawtoothlike oscillation centered at ~ 
380 meV in CH 4 (Fig. [20]). While the IR spectrum in- 
dicates strong absorption at this energy, the magnitude 
and shape of this feature in Z e g are not consistent with 
the VFR observed in other molecules. Such sawtooth 
featur es are also observed in th e Z e s spectra of CO2 and 
ITO (lYoung and Surkol l2008d) . 

The origin of these features is at present unclear. 
iNishimura and Gianturcol ( 2005al) predicted that water 
should have fairly strong vibrational excitation cross sec- 
tions, with sharp onsets. Channel coupling could result in 
additional structure in the Z a spectrum near the vibra- 

? — 1 I K 

tional excitati on thresholds (lYoung and Surkol 12008a) . 
Alternatively, lYoung and SurkoM 2008c ) suggested that 
these features could represent interference between di- 
rect and resonant annihilation. However, the latter is 
not compatible with the evidence that CO2, H2O, and 
CH4 do not to bind positrons. 

Other non-VFR-type features are either predicted 
or allowed by theory but have yet to be observed. 
One example is a shape resonance which could occur 
if the positron became temporarily trapped inside a 
positive-energy potential barrier. However, low-energy 



positron scattering and annihilation is usually domi- 
nated by the s-wave component of the incident positron 
wave function, which has no centrifugal barrier. Fur- 
ther, the atomic cores are repulsive, so a shape res- 
onance is unlikely. Exceptions in other systems in- 
clude the prediction of a p-wave positron shape res- 
onance in the Mg atom (Mit rov and Bromlevl . l2007t 
iMitrov et all [2008) and cage-stat e shape resonances 
in cubane (CsH s ), C2n and Cm dCarev et all 12005 



19991: IGianturco et all I200E 



IGianturco and Lucchese.. . _ . 

Finally. INishimura and Gianturcol (|2003l ) hypothesized 
that the presence of a virtual state, by itself, can lead 
to a long-lived intermediate state (e.g., following a vi- 
brational de-excitation collision), and this is expected to 
produce a broad spectral feature. Experimental investi- 
gation of these predictions is warranted. 



F. Small molecule summary 



The theory of iGribakin and Led (l2006bl) (Sec. HPT) 
provides a remarkably useful framework for understand- 
ing resonant annihilation in small molecules. The rel- 
ative contribution of each VFR is a competition be- 
tween resonant elastic scattering and annihilation. As 
long as the coupling to the positron continuum is strong 
enough, the magnitude of each resonance is proportional 
to g = •v/eb/e, and the annihilation spectrum is a sum of 
distinct resonances. 

In the case of the methyl halides and ethanol, the con- 
ditions for the simple g scaling were confirmed by direct 
calculation of the capture rates using the Born-dipole- 
type approximation (Sec. III.F|) . For methyl halides, the 
binding energies were extracted from the measurements 
on the protonated species. They were then used to make 
a prediction of Z e ff for the deuterated species, thus pro- 
viding a stringent test of the theory. Agreement between 
experiment and theory is excellent. 

In other cases, however, the VFR of the IR- active 
modes are insufficient to explain the observed annihi- 
lation spectra. In methanol, it is necessary to include 
dipole-active multimode vibrations; and in ethylene and 
acetylene, one must also include IR-inactive modes and 
n-quantum overtones and combinations to match the ex- 
perimentally measured Z e g. In these cases, the VFR 
of overtone and combination vibrations had to be ad- 
justed arbitrarily (i.e., by factors 1/n) to match the ex- 
perimental Z c ff. Thus, while the theory of Sec. [IFF] pro- 
vides useful insights into the vibrations that are likely in- 
volved, it lacks quantitative predictive accuracy for many 
molecules. 

This theory has also helped to elucidate the boundary 
separating species that have enhanced Z c g (i.e., beyond 
that predicted by single resonances) from those that do 
not. The factors responsible for this transition appear 
to be molecular size and the density of vibrational states 
near a resonance (Sec. III.Gj) . This enhancement is likely 
due to IVR. There is some evidence that it is already 
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operative in molecules as small as ethane and is dominant 
in propane. 



V. IVR-ENHANCED RESONANT ANNIHILATION IN 
LARGER MOLECULES 

A. Overview 

As indicated in Table U annihilation rates for large 
molecules grow rapidly with molecular size. Similar to 
the case of small molecules, these large rates are under- 
stood to be the result of positron attachment via VFR. 
As shown in Fig. Q] for butane, the Z e R spectra have dis- 
tinct peaks, the positions of which are strongly corre- 
lated with those of the IR peaks but are shifted down- 
ward by the binding energy. As was the typical case for 
small molecules, the VFR in larger molecules also occur 
predominantly at energies corresponding to those of the 
fundamental vibrations. The low-energy plateau in the 
alkane spectra (e.g., below 140 meV in Fig. [1} is due to 
C-H bend modes and C-C modes, while the high-energy 
peaks are due to C-H stretch modes. In the case of large 
molecules, however, the Z e ff values in these peaks are en- 
hanced by orders of magnitude above those predicted for 
VFR in small molecules. 

To validate the identification of the resonant peaks 
with the vibrational modes, shown in Fig. [22] is the 
spectrum of butane and nonane and the corr espond 
ing fully deutera t ed com pounds (IBarnes et al 



respond- 
L 2003; 



lYoung and Surkol l2008bl ). When the spectra of the 
deuterated compounds are corrected for the change in 
vibrational mode frequencies, the deuterated and hydro- 
genated data are in good agreement. 

As discussed in Sec. [TT1 a plausible explanation for the 
large magnitudes of Zeff hi large molecules is that, by 
the excitation of a vibrational fundamental, a large num- 
ber of otherwise inaccessible multimode VFR (so called 
dark states) contribute to the annihilation, mediated by 
the process of intramolecular vibrational energy redistr i- 
bution (IVR) (|Gribakinl |2000t iGribakin and G ill 12004ft . 
In this paradigm, the incident positron is first captured 
into a mode-based doorway state (e.g., involving infrared 
active modes) that then couples to a bath of quasidegen- 
erate multimode dark states. This increases the multi- 
plicity of the final capture states and causes an approxi- 
mately proportionate enhancement in the resonant anni- 
hilation rate. As discussed below, it appears that the fun- 
damental vibrations act as doorways in large molecules. 
This in turn leads to the excitation of some, but not 
all of the nearby multimode states in the IVR process, 
which results in enhanced annihilation. A similar cou- 
pling to dark states occurs in m olecular photoabsorption 
([Stewart and McDonald Il983l ). 

This doorway-state model provides a common thread 
in the discussion of annihilation for large molecules. A 
particularly important link in this model is the connec- 
tion between the observed rapid rise in the resonant an- 
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FIG. 22 Z e g spect r a of hydrogenated butane (solid 
circles) llBarnes et al\, 12003ft and nonane (open circles) 
(|Young and Surkol l2008bft compared with (the appropriately 
scaled) fully deutera ted analogs (dashed and chain curves) 
(|Barnes et al\ , [2003). Mode energies for the deuterated 
species have been scaled by the appropriate reduced mass 
factor to match those in the hydrogenated compounds, while 
the positron binding energy was assumed to be independent 
of deuteration. The magnitudes of Z e g have been scaled by 
the appropriate g factor at each incident positron energy. 



nihilation rate with molecular size and IVR. A number 
of experiments and analyses clarify and, in some cases 
quantify, this physical picture. 



B. The alkane molecule paradigm 

Alkane molecules have been studied more extensively 
than any other molecular species (cf. Figs. HI and l2"2"j) . 
Figure [53] shows six examples. These spectra resemble 
closely the spectra of the fundamental vibrations, al- 
beit downshifted by the positron-molecule binding en- 
ergy. The binding energies increase by 20-25 meV for 
each carbon-based monomer added to the alkane chain. 
Thus, the binding energy for propane is 10 meV, while 
for nonane it is 145 meV. The magnitudes of the resonant 
annihilation peaks grow rapidly with molecular size, in- 
dicating that IVR, albeit likely incomplete, is operative. 

The theory of resonant annihilation for small molecules 
fSec. III.F]) prescribes a restricted role for binding energy 
in determining Z c g spectra, namely, 



Z c s oc b v g{e) = b v \Je h je 



(58) 



where b v is the multiplicity of the excited resonant cap- 
ture states. The vibrational coupling appears only via 
the positron capture rate T e , which usually cancels the 
total rate T for all but very weak resonances. 

As shown in Fig. [M] the Z e g spectra for alkanes larger 
than ethane are self-similar when they are first scaled by 
the factor g(e) , then the resulting magnitudes of Z e g /g(e) 
are normalized at the C-H stretch peaks, and, finally, 
the spectra are shifted up ward by their binding energies 
(lYoung and Surkol l2008bft . A surprising result from this 
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FIG. 23 Z c a spectra as a function of incident positron en- 
ergy for a variety of alkane molecules, showing the systematic 
shift of the spectra to lower positron energy as molecular size 
is increased. Values of Z e s are absolute except for C14H30, 
which is in arbitrary units due to difficulties measuring its 
vapor pressure. 
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FIG. 24 Normalized and energy-shifted Z e s/g spectra for 
alkanes with n = 3-8 carbons. For comparison, room- 
temperature Z e g/g data for alkanes are also shown (circles 
with plus sign) at energies Eb + £t, where et = \ksT = 
37.5 meV. Each room temperature datum is labeled by the 
number of carbons n in the molecule. See text for details. 



analysis, not presently understood, is that the relative 
magnitudes of the high-energy and lower-energy peaks 
in the alkanes remain the same over a factor of 3 in 
molecular size and a factor of 10 3 in the magnitude of 
Z c g. The resulting self-similar spectra shift downward 
with increasing molecular size as the positron binding 
energy increases. As discussed in Sec. I VIII this has the 
consequence that the i? e ff spectrum, measured with ther- 
mal positrons, can be assigned to the vibrational modes 
populated through the corresponding VFRs in the ther- 
mal energy range. This connects in a quantitative way 
the energy-resolved measurements and thermal measure- 
ments of Z e g. 

As the size of the alkane is increased beyond twelve 
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FIG. 25 Z e ff spectra for (o) dodecane, (diamonds) tetrade- 
cane, and (squares) hexadecane. The vertical arrows indicate 
the positions of the C-H stretch mode VFR peaks for the sec- 
ond bound state (i.e., the first positronically excited state) in 
each molecule. The spectra for tetradecane and hexadecane 
have been normalized arbitrarily since their vapor pressures 
were too low to measure reliably. The large peaks at lower 
energy are the C-H stretch mode VFR for the first bound 
states (the positronic ground states). 



carbons (dodecane), a new feature appears in the Z e g 
spectrum at an energy close to that of the C-H stretch 
peak. As shown in Fig. [55j the magnitude of this peak 
grows as the alkane size increases, and it shifts downward 
in energy, just as the original C-H stretch peak. In hex- 
adecane (the 16-carbon alkane), the C-H peak occurs at 
a mere 55 meV incident positron energy, corresponding 
to a binding energy of 310 meV. 

This new resonance is attributed to a second positron 
bound state (i.e., the first positronically excited state) on 
the molecule. It is populated by a C-H stretch peak VFR 
in a manner similar to the larger ground-s tate peak that 
occur s at smaller incident positron energy ( Barnes et all 
120061) . The small peak at 365 meV in dodecane is iden- 
tified as the C-H stretch peak of the first, positronically 
excited bound state with a binding energy of a few milli- 
electron volts, while the larger peak at 150 meV is due to 
the positron in its ground state. The positions of these 
peaks are in good agreement with a model calculation 
described in Sec. |HJT2J (cf. Fig.©. 

An important feature of these positronically excited 
resonances is their magnitude relative to the correspond- 
ing ground-state resonance in the same molecule. For 
small molecules, it is expected that the contribution of 
a resonance will be proportional to g — y^eb/e. For 
both tetradecane and hexadecane, the ratios of the mag- 
nitudes of the first and second bound-state C-H stretch 
resonances are equal to the ratios of the g factors for 
these resonances. In particular, since the positron over- 
lap density is proportional to y/eb (Sec. III.H.2|) . it is ex- 
pected to be smaller for the second bound state, due to 
the smaller binding energy. This scaling of Z e g for the 
positronically excited states and the scaling with g of the 
spectra shown in Fig. [M] both demonstrate the impor- 
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tant role of the factor g in determining the magnitudes 
of annihilation peaks in large molecules. 



C. Dependence of Z cH on molecular size 

As shown in Fig. [26l positron-molecule annihilation 
has been studied for a variety of chemical species. 
While the vibrational modes and energy levels in these 
molecules differ in various ways, all of the hydrocarbons 
studied contain strong C-H stretch vibrational modes 
that result in prominent annihilation resonances. Conse- 
quently, they provide a convenient benchmark to deter- 
mine positron-molecule binding energies and the relative 
magnitudes of the annihilation rates. 
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FIG. 26 Z e fj at the C-H stretch peak vs binding energy Eb 
for alkanes, C„H2 n +2 (D), with the number of carbons n 
indicated; rings (hexagons); halomethanes (©); ethylene (•); 
methanol (+); 1-chlorohexane (A); fluoroalkanes (A); deuter- 
ated species (v)- 

Relatively small changes in molecular structure can 
have significant effects on both the positron binding ener- 
gies and the overall annihilation rates. For example, the 
C-H stretch peak in 1-chlorohcxane is shifted downward 
relative to that of hexane by nearly 100 meV, while the 
magnitude of the C -H stretch resonance is in creased by 
nearly a factor of 3 ( Young and Surkol 12008 bfh 

Based upon the alkane data, one might conclude that 
binding energy and the magnitudes of the annihilation 
resona nces are strongly correlated ( Young and Surkol 
l2008bD . However, as shown in Fig. [26l this hypothesis 
is not borne out by the data. With the exception of the 
alkanes, the binding energies and the magnitudes of the 
C-H stretch resonances appear to be only weakly corre- 
lated. For example, the Z e R values for 1-chlorohexane, 
benzene, and naphthalene are an order of magnitude or 
more smaller than those for alkanes with similar binding 
energies. 

This lack of correlation of Z e g with binding energy 
motivated further analysis of the Z c g magnitudes. Shown 
in Fig. [27] and listed in Table [V] are data for Z c s the C- 
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FIG. 27 Z e g at the C-H stretch peak, normalized by the fac- 
tor g = ^eb/s, vs the number of atoms in the molecule. The 
symbols are as in Fig. 1261 The solid line is the fit given by 
Eq. (|59[) . The error bars on 1-fluorohexane indicate the un- 
usually large uncertainty for this datum. 



H stretch peak, normalized by <?, as a function of the 
number of atoms N. With the exception of the partially 
fluorinated compounds, which will be discussed below, 
the magnitudes of the annihilation rates for the relatively 
wide variety of molecules studied lie close to a universal 
curve. 1-chlorohexane and benzene are no longer outliers 
to the extent they were in Fig. |21)| The empirical scaling 
of Z e s with N is found to be ( Young and Surkol 120071 , 
l20Q8bh . 



Z eS /g = 2.3JV 9 , 



(59) 



with q = 4.1 ± 0.1. This scaling likely reflects the depen- 
dence of the total number of accessible positron-molecule 
vibrational states (i.e., including dark states) on the num- 
ber of vibrational degrees of freedom. 

It is not surprising that Z e ff depends linearly on g = 
n/k, so long as the weak binding picture (i.e., s C 1 
a.u.) is valid. Furthermore, the 1/k factor arises from 
the normalization of the incident positron wave func- 
tion. The surprising aspect of this g scaling is that it 
appears to incorporate the only dependence of Z e g on 
£6. While the overall vibrational density of states is cor- 
related with the number of atoms N, it can certainly 
change without N changing (e.g., by chemical substitu- 
tion). Further, small- molecule theory assumes that the 
lifetime of a vibrationally excited positron-molecule com- 
plex is short compared to the annihilation time. How- 
ever, if the states accessible by IVR are especially long 
lived, the magnitudes of the resonances might then de- 
pend upon the capture rate rather than g. The annihi- 
lation rate would then be expected to saturate, growing 
only linearly with molecular size. The fact that the scal- 
ing of Eq. (|5"9"|) remains valid for the largest molecules 
studied to date is evidence that this saturation limit has 
not yet been reached. 
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TABLE V Parameters for annihilation on various large 
molecules: the number of atoms N; positron binding energies; 
Z e g at the C-H stretch mode and for thermal distributions of 
positrons at 300 K; and the normalized rate, /g. The 

Z e s values have an overall uncertainty of 20% (cf. Sec. lIII.E|l . 



Species 


N 


(meV) 


7 (CH)l 




Z l ? t H) /g 

oft ' 3 






Alkancs 






IVlcthanc 


5 


< 




142 




Ethane 


8 


> 


900 


660 




Propane 


11 


10 


10 500 


3 500 


63 000 


Butane 


14 


35 


21000 


11 300 


65 000 


Pcntanc 


17 


60 


80 000 


37 800 


180 000 


Hcxanc 


20 


80 


184 000 


120 000 


350 000 


Heptane 


23 


105 


590 000 


242 000 


930 000 


Octane 


26 


115 


1 090 000 


585 000 


1 610 000 


Nonanc 


29 


145 


2 000 000"* 


643 000 


2 500 000 


D o d c c an c 


38 


220 


9 800 000 


1 780 000 


8 000 000 


2nd B.S. 


38 


> 


1200 000 






'i/ol radecanc 


44 


200 


llx" 




6.8a; e 


2nd B.S. 


44 


50 


2.8x" 




7.0x e 


Hcxadccanc 


50 


310 




2 230 000 


6.4y e 


2nd B.S. 


50 


100 


4.0y° 




6.5y 






Alkanc isomers 






Isopcntanc 


17 


60 


80 000 

Rings 


50 500 


180 000 


Cyclopropane 


9 


10 


3 600 




21 500 


Cy clo hcxanc 


18 


80 


94 000 


20 000 


180 000 


Benzene 


12 


150 


47 000 


15 000 


58 000 


Pn apht halcnc 


18 


300 


1 240 000 


494 000 


640 000 






Haloalkancs 






l-fiuoropropanc 


11 


30 


1520 




5100 


2,2-difiuoropropanc 


11 


25 


900 


8130 


3 300 


1-fiuorobutanc 


14 


70 


5 600 




11500 


1-fiuorohcxanc 


20 


105 


60±30xl0 3 


269 000 


94 ± 47 x 10 3 


1-chlorohcxanc 


20 


175 


520 000 




540 000 


1-fiuorononanc 


29 


100 


930 000 




1050 000 






Dcutcratcd 






d-benzene 


12 


150 


61000 


36 900 


57 500 


d-butanc 


14 


35 


28 500 




75 000 


d-nonanc 


29 


145 


2 400 000 


641 000 


2 300 000 


d-naphthalcnc 


18 


~ 300 









a Determincd from the en ergy shift of the C-H peak in Z e g . 

6 At the C-H stretch peak feaxnes et ^ , l200ll200lil ; lGilbert et all 
120021 ; lYoung and Surkol 120071 , l2008bl1 . 

To r thermal positrons at 300 K l lH evland et all fl982l ; llwatal , 
119971 : llwata et all |1995| , 11994 Il997al ; Wright et ad , ll983f) . 

d zfff H ^ for nonane i s from lYoung and Surkol {2008b) as opposed 
to lBarnes et al. as the latter had a narrower positron energy 

distribution for this molecule. 

E Absolute Z e ff could not be determined, so the values are multi- 
plied by arbitrary factors x and y. 



D. Toward a model of annihilation in large molecules 

Our present theoretical understanding of the annihi- 
lation spectra of large molecules is incomplete. As dis- 
cussed in Sec. HH it is not possible to explain the large 
values of Z e g that are observed using only VFR involving 
the vibrational fundamentals. This difficulty could pos- 
sibly be overcome by considering mode-based resonances 
as doorways for positron capture in complex multimode 



VFRs. However, estimates assuming complete IVR show 
that Z e g would increase much faster with molecular size 
than is observed (Sec. III.G.1|) . Further, when all com- 
bination and overtone vibrations are assumed to cou- 
ple to the positron continuum, the resulting spectrum 
is predicted to be featureless, bearing little resemblance 
to the experimentally measured spectra. Nevertheless, 
the presently available theoretical models do provide a 
useful framework with which to interpret some of the 
experimental results. In the following, additional experi- 
ments and analysis are discussed that elucidate details of 
VFR-mcdiated annihilation in large molecules and place 
constraints on viable theoretical models of this process. 



E. Inelastic autodetachment 

Inelastic escape channels are a potentially important 
mechanism that can limit the magnitudes of resonant an- 
nihilation peaks (Sec. III. Gj) . They can occur when a reso- 
nantly captured positron is released from the molecule by 
the de-excitation of a vibration other than that produced 
during the initial capture. Such vibrations can either be 
excited through IVR following resonant capture or by 
thermal excitation. This process is expected to lead to a 
reduction in the annihilation rate, since the positron will 
spend less time on the molecule. 

An important consideration in such an inelastic pro- 
cess is how the resonant Z e g is affected by the binding 
energy. If the positron is in a weakly bound state, many 
vibrational modes will have sufficient energy to eject it; 
while for deeply bound states, fewer modes are able to do 
this. Thus, if inelastic escape channels were present, one 
would expect an additional dependence of the Z e g on e^, 
beyond the scaling with g. However, this contradicts a 
number of experimental results. For example, the C-H 
stretch-peak magnitudes for the first and second bound 
states of tetradecane and hexadecane strictly follow the 
g scaling with no additional dependence on Sb- Further- 
more, the Z e g at the C-H stretch peaks for nearly all 
molecules studied adhere strictly the Z e g/g oc iV 41 scal- 
ing, with relatively little deviation due to their disparate 
binding energies. This appears to rule out inelastic es- 
cape channels as being generally important in determin- 
ing Z e g values. 



1. Fluorine-substituted alkanes 

There is, however, one notable exception. As shown 
in Fig. [27l partially fiuorinated alkanes, such as the 1- 
fluoroalkanes, deviate significantly from the scaling of 
Eq. (|59|) . The data indicate that the substitution of 
a single fluorine atom for a terminal hydrogen in a given 
alkane reduces the height of the C-H stretch peak in 
Z c g by as much as an order of magnitude. The C-F 
stretch mode might be expected to play a significant role 
in this apparently inelastic process since the cross sec- 
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FIG. 28 The Z e g spectra for ( o) pro pane and (•) 1- 
fluoropropane (| Young and Surkol. l2008bl ). and (A) 2,2- 
difiuoropropane ( Gilbert et all 120021 ). The vertical dot- 
dashed lines indicate the expected energies of C-F stretch res- 
onances and escape thresho lds in each molecul e, based on C-F 
stretc h mode energies from iDurig et al\ (|l98ll ) ; iGuirgis et al\ 
(|l999h . See text for details. 



tion for positron-impact excitation o f the C-F stre t ch vi- 
bration in CF 4 is un usually large ( Marler et all l2006t 
ISullivan et aZI . l2002af) . 

Such an inelastic channel can significantly reduce the 
magnitude of Z c e by increasing the total post-capture 
escape rate. In this case [cf. Eq. (02)], 



Z cS (v) oc 



rj. 



(60) 



where Y % v is the inelastic escape rate (e.g., via de- 
excitation of mode n with energy so that increasing 
Tl reduces Z c g. This process requires that the multi- 
mode state v contains quasidegenerate components in 
which mode n is excited. This restriction imposes a 
threshold at s = ui n — £b, above which Z e g is reduced. 

As shown in Fig. [28l the Z e s spectra for 1- and 2,2- 
fluoropropane display just such a suppression at larger 
incident positron energies as compared with the analo- 
gous hydrogenated compounds. In 1-fluoropropane, the 
C-F stretch annihilation resonance is expected to occur at 
~ 90 meV, assuming Eb = 30 meV and ujcf ~ 120 meV 
( Guirgis et all Il999[) . In 2,2-difluoropropane, the C-F 
stretch peak is expected to occur at ~ 125 meV, assum- 
ing E b = 25 meV and wcf ~ 150 meV (jDurig et all 
Il98ll ). In both cases, the suppression of Z e g occurs near 
these threshold energies (shown by dot-dashed lines in 
Fig. 

As discussed by lYoung and Surkol ()2008bl ). there are 
similar decreases in the magnitudes of the Z c g at higher 
incident positron energies in the larger partially fluo- 
rinated alkanes, 1-fiuorobutane, 1-fiuorohexane, and 1- 
fluorononane, relative to the analogous hydrogenated 
molecules. However, there is also an increase in the mag- 
nitudes of the Z e g spectra at smaller positron energies, 
making the thresholds for escape-channel behavior less 
distinct. The behavior of 1-fiuorononane is more difficult 



to explain. The positron binding energy of 160 meV in 
this molecule would app ear to preclude the C-F stretc h 
mode [w„ < 127 meV (|Crowder and Lightfootl . Il983l) ] 
from acting as an inelastic escape channel. However, a 
decrease in annihilation is observed at higher energies 
nonetheless ( Young and Surkol l2008bl ). This is currently 
not understood. 



2. Effects of molecular temperature on Z c g 

For molecules at finite temperatures, the energy for 
positron detachment which suppresses Z e g can be sup- 
plied by pre-existing thermally excited modes rather 
than those excited in the attachment process. Thus 
Z e fi might be expected to increase significantly with 
decreasing molecular temperature. At one point, it 
was thought this could explain the empirical obser- 
vation tha t Z e g oc (2n + 2) exp^/fee? 1 ) in alkane 
molecules ( Barnesl |2004[ ). A qualitatively different ef- 
fec t of varying molecular temperatu re was suggested 
by iNishimura and Gianturcol ( 2005b| ). They proposed 
that vibrational excitation and the associated change in 
molecular geometry could be required to induce or in- 
crease positron binding to molecules, and this, in turn, 
could affect the activation of annihilation resonances. In 
this model, one would expect that increasing the molec- 
ular temperature should result in an increase in Z e g. 

Experiments were done to test these ideas using a 
specially constructed cold cell, so that Z c g spectra 
could be measured at dif ferent molecular temperatures 
(lYoung and Surkol l2008ah . The apparatus is described 
in Sec. IIII.EI Care was taken to ensure that the test gas 
in this flowing-gas system actually cooled to the ambient 
temperature of the cold cell. In addition, the test-gas 
pressure was maintained a safe margin below the equi- 
librium vapor pressure at each temperature in order to 
avoid condensation on surfaces near and inside the cell. 
The spectra for pentane at 153 and 300 K are shown 
in Fig. [29l There is only a small (~10%) increase in 
the magnitude of the C-H stretch resonance with the 
change in molecular temperature, while at lower inci- 
dent positron energies, there is a somewhat larger in- 
crease in Z e ff (i.e., ^30%). Similar results were ob- 
tained for heptane at 195 and 300 K, but the increase 
in the low-energ y portion of the spectrum was somewhat 
larger (~50%) ([Young and Surkol l2008al) . The spectra 
of both molecules indicate that their binding energies do 
not change with changes in the molecular temperature. 

These experimental results indicate that there is 
clearly no Boltzma nn- facto r- like d ependence of Z c g, as 
that considered by iBarne 1 (120041) . If there were, the 
pentane C-H stretch peak would have grown by a fac- 
tor of 10. These results also tend to rule out the sug- 
gestion, made by INishimura and Gianturcol ( 2005bl ) for 
C2H2, C2H4, and C2Hg that thermal deformation of the 
molecular bonds can enhance significantly the binding 
energy and hence increases the rate of VFR-mediated 
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FIG. 29 Energy-resolved measurements of Z e s for pentane 
(C5H12) at 300 K (o) and 153 K (•) using a cold cell 
(jYoung and Surkoll2008al ). 



positron attachment and annihilation. The data show 
the opposite trend, namely, that increasing the molecu- 
lar temperature leads to smaller values of Z e &. 

These findings also confirm other results indicating the 
absence of thermally excited escape channels. For ther- 
mally excited modes to provide effective escape, their 
energies must exceed the positron binding energy. As 
the binding energy increases (e.g., with alkane size), one 
would expect these channels to be switched off. Hence 
thermally activated detachment would produce an addi- 
tional dependence of the resonant Z c g on Eb beyond the 
g factor. Such an effect is not apparent in the alkanes, 
where the empirical scaling is described by Eq. (1591) for 
both the positron ground and first excited states. 

A question remains as to why the temperature effect 
appears to be stronger at low impact energies than at 
high energies, yet the ^-normalized spectra for alkanes at 
300 K (cf. Fig. [2~4"|) remain self-similar. 



F. Other IVR-related phenomena 

While the iV 41 scaling provides strong evidence of 
IVR, there are still significant questions as to how the 
IVR, induced by positron capture, proceeds in large 
molecules. We discuss here a few of the outstanding is- 
sues. 

The partially fluorinated alkanes (Sec. IV.E.1[) pro- 
vide a clear example of post-capture vibrational energy 
transfer, which suppresses the Z e g spectra due to in- 
elastic detachment. The so-called "intermediate" mul- 
timode state containing the excited C-F stretch mode is 
suggestive of a tiered IVR model, in which vibrational 
energy redistribution occurs i ncrementally through a n 
ever-growing set of dark states ( Nesbitt and Field] , 1 1996(1 . 
This could explain why only a few multimode excita- 
tions involving the C-F stretch (out of all possible exci- 
tations) have such a disproportionate influence. A simi- 
lar tiered model has been used to describe laser- e xcited 
vibrational dynamics in phenols ( Yamada et all |2007[ ) 
and to calculate IVR rates for acetylenic stretch modes 
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FIG. 30 Energy-resolved -Z e ff (•) and infrared absorption 
(solid curve) for benzene. The Z c ft spectrum has been 
shifted upward by the binding energy (Eb = 150 meV) for 
direct comparison. Th e norm alization of the IR absorption 
(jLinstrom and Mallard! . 120051 ) is arbitrary. Vertical lines in- 
dicate the positions of the vibrational modes. 



(|Stuchebrukhov and Marcus! . [l993h . 

In general, the haloalkanes have larger values of ther- 
mal Z c g at 300 K than their hydrogenated counterparts. 
There is also evidence that thermal values of Z e g for 
partially deuterated benzenes and other substitu ted ben- 
zenes are similarly enhanced ( Iwata et all Il995l ) . These 
observed increases in Z e g may be due to the increase in 
the density of vibrational "dark" states. Similar physics 
may account for the large thermal Z e ff values in CCI4 
and CBr 4 (9000 and 40 000, respectively). 

It is not clear at present what makes a "good" vibra- 
tional doorway state. In most large hydrocarbons, it is 
only the fundamental vibrations that appear to produce 
VFR; however, there are exceptions. Shown in Fig. [30] 
is the Z e g spectrum of benzene, shifted upward by its 
binding energy and compared with the infrared absorp- 
tion spectrum. Note the distinct peak at ~235 meV in 
the shifted Z c g spectrum. While there are no nearby fun- 
damental vibrations, there are two IR-active combination 
vibrations at 227 and 244 meV in the IR spectrum. Thus, 
the additional peak in Z c g appears to be evidence of un- 
usually strong capture into multimode doorways, likely 
enhanced by IVR. 

In elucidating the role of IVR it is of interest to com- 
pare propane and cyclopropane. Figure [31] shows that 
transforming the propane molecule into a ring reduces 
the C-H stretch peak by a factor of 3, approximately in 
agreement with the iV 4 scaling in Eq. (|59"1) . At the same 
time the plateau at lower energies is narrowed so that 
only one broad peak remains. However, its magnitude in 
cyclopropane is practically the same as in propane. Thus 
the cyclopropane spectrum does not obey the self-similar 
scaling observed in the alkanes. It is possible that the 
low-energy peak contains mode-based VFRs, with little 
enhancement due to IVR. The reduction in the magni- 
tude of the C-H peak may be due to the absence of low- 
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FIG. 31 En ergy-resolved Z e g (• ) for (a) propane and (b) cy- 
clopropane (|Barnes et all I2006D . The solid curves are the 
normalized IR absorption spectra, an d the vertical bars below 
each plot show the vibrational modes (|Linstrom and Mallard! . 
2005). Shown as insets are the molecular structures. 



frequency modes in cyclopropane. This is consistent with 
Table IIIII which shows that cyclopropane has a markedly 
lower vibrational density at the C-H stretch energy as 
compared with propane. 

As shown in Fig. [3lJ there is an identifiable feature in 
cyclopropane at ~ 250 meV (i.e., in the gap between the 
C-H stretch and lower-energy modes) that does not ap- 
pear in propane or larger hydrocarbons. It occurs in the 
energy range where there are peaks in the IR spectrum 
and it is likely due to combination and overtone vibra- 
tions (similar to that observed in benzene). This again 
points to the possibility that appreciable IR coupling is a 
predictor of the strength of multimodc doorways in large 
molecules. In cyclopropane this peak is not particularly 
enhanced (e.g., relative to the C-H stretch resonance), 
and looks closer to the effects due to combination and 
overtone VFRs observed in small molecules (see Sec. IIV1 
e.g., ethylene, Fig. IT71) . This phenomenon also bears fur- 
ther scrutiny. 



G. Large molecule summary 

There are several defining characteristics of the Z e g 
spectra of large molecules. They exhibit a series of peaks, 
the positions of which bear close resemblance to those of 



the fundamental vibrations, downshifted by the positron- 
molecule binding energy. The amplitudes of these reso- 
nances grow rapidly with molecular size but exhibit a 
relatively weak dependence on binding energy and inci- 
dent positron energy via the factor g. The amplitudes 
of the C-H stretch resonances in hydrocarbons obey the 
empirical scaling of Eq. (fBTJ)) . This scaling likely reflects 
the rapid increase in the number of vibrational degrees 
of freedom of the molecule with molecular size. It sug- 
gests that the VFR are enhanced by the IVR process. 
However, the extent of this IVR appears to be limited 
and does not reach the statistically complete final state. 
Inelastic escape channels appear to be relatively unim- 
portant, at least in the hydrocarbons studied thus far. 

Our theoretical understanding of annihilation pro- 
cesses in large molecules is not nearly as well developed 
as it is for small molecules. Qualitatively, a positron ap- 
pears to excite a vibrational fundamental (or in some 
cases, such as benzene, a combination or overtone) and 
populates a doorway resonance. The vibrational energy 
imparted to the molecule can then flow to some set of 
quasidegenerate dark states. However, if this picture is 
correct, there must be an operative escape channel. If 
not, all VFR will give approximately the same contribu- 
tion to the annihilation rate. Considering the available 
experimental evidence, it seems plausible that the reso- 
nant Z e f{ magnitudes are determined by a balance be- 
tween diffusion to some (limited) set of dark states, and 
positron ejection through either the doorway state that 
it entered or other nearby doorway states in what might 
be termed a "quasielastic" escape process. 

Beyond this qualitative picture, there are an assort- 
ment of relatively poorly understood effects. The most 
prominent of them is arguably the self-similar spectra ob- 
served in alkanes. There is also the occasional occurrence 
of multimode VFR, and the inelastic detachment ob- 
served in partially fluorinated alkanes. These effects pro- 
vide tantalizing hints of the many types of complex inter- 
nal dynamics that can be responsible for VFR-enhanced 
annihilation depending upon the particular situation. At 
a minimum, further investigations are warranted. For 
example, experiments with higher positron energy reso- 
lution would be useful to determine additional details of 
the processes involved. 



VI. POSITRON-MOLECULE BINDING ENERGIES 

A. Relation to molecular properties 

As described in Sees. IIVI and [Vj positron-molecule 
binding energies have now been measured or esti- 
mated for about 30 molecules (Barnes et all 2003. 2006; 
iGilbert et all 120021 : lYoung and Surkd. l2008bll^Y~ These 
data were analyzed recently by iDanielson et all ( 2009f ) 
with a view toward determining which molecular param- 
eters affect Eh . Table I VII lists the available binding en- 
ergies together with molecular dipole polarizabilies ad, 
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dipole moments /x, and ionization energies Ei 



energies continue to grow with n. 



TABLE VI Positron binding energies, normalized thermal 
annihilation rates Z^g /Z, and physical parameters for selected 
molecules. 



Formula 



e b b Z Z$ { /Z c 



Ei " 









(moV) 


(mcV) 






(A 3 ) 


(D) 


(eV) 






Alkancs and 


related 


molecules 








Ethane 


c 2 


Ho 


> 


-14 


18 


37 


4.4 





11.5 


Propane 


c 3 


H 8 


10 


10 


26 


13.5 


6.3 


0.1 


11.1 


Butane 


C4 


H10 


35 


31 


34 


330 


8.1 





10.6 


Pentane 


c 5 


H12 


60 


54 


42 


000 


10.0 





10.4 


Hcxanc 


Ce 


H14 


80 


77 


50 


2400 


11.8 





10.2 


Heptane 


C 7 


Hie 


105 


100 


58 


4200 


13.7 





9.9 


Octane 


C s 


His 


115 


123 


66 


8800 


15.5 





10.0 


Nonanc 


c 9 


H20 


145 


146 


74 


8700 


17.4 





10.0 


Dodccane 


Ci 


2H26 


220 


214 


98 


18000 


22.9 





9.9 


Tetradccanc 


Ci 


4H30 


260 


2G1 


114 


- 


26.6 





9.9 


Hcxadccanc 


Ci 


6H34 


310 


300 


130 




30.3 





9.9 


Butanc-dlO 


c 4 


Did 


35 


31 


34 




8.1 





- 


Nonane-d20 


c 9 


D20 


145 


146 


74 


8700 


17.4 





- 


Acetylene 


c 2 


H 2 


> 


-28 


14 


230 


3.3 





11.4 


Ethylene 


c 2 


H 4 


20 


-17 


16 


75 


4.2 





10.5 


Isopentanc 


c 5 


H12 


60 


57 


42 


1200 


10.0 


0.1 


10.3 


Cyclopropane 


c 3 


H 6 


10 


0.7 


24 




5.7 





9.9 


Cyclohcxanc 


Ce 


H12 


80 


68 


48 


420 


11.1 





9.9 








Aromatics 












Benzene 


c 6 


Ho 


150 


149 


42 


360 


10.4 





9.3 


Bcnzcnc-d6 


c 6 


D 6 


150 


149 


42 


730 


10.4 





9.3 


Naphthalene 


Ci 


oHs 


300 


290 


68 


7300 


16.6 





8.2 








Alcohols 












Methanol 


CH30H 


2 


5 


18 


84 


3.3 


1.7 


10.9 


Ethanol 


c 2 


H5OH 


45 


27 


26 




5.1 


1.7 


10.5 




Partially halogcnatcd hydrocarbons 








Methyl fiuor. 


CH 3 F 


> 


-3 


18 


77 


2.4 


1.85 


12.9 


Methyl chlor. 


CH3CI 


25 


23 


26 


580 


4.4 


1.9 


11.2 


Methyl brom. 


CH 3 Br 


40 


35 


44 




5.6 


1.8 


10.5 


1-fl. propane 


c 3 


H 7 F 


30 


45 


34 




6.0 


2.0 


11.3 


2,2-difl.prop. 


c 3 


H 6 F 2 


25 


51 


42 


190 


5.9 


2.4 


11.4 


1-fl. butane 


C 4 


HgF 


70 


27" 


42 




7.8 






1-fi. hcxanc 


c 6 


H13F 


105 


73" 


58 


46000 


11.5 






1-fi. nonanc 


c 9 


H19F 


160 


141 e 


82 




17.0 






1-chl.hcxane 


c 6 


H13C1 


175 


138 


66 




13.6 


2.0 


10.3 



"Measured in energ y-resolved annihilation experiments 
j Young and Surkol . l2008bl ld) . 
'Predicted by Eq. JBll . 

c Z^f i s measured at room temp erature with thermalized positrons 
ijlwatal . H997I : llwata et al.l. 11993 ). 

d Po l arizab i lities, dipo l e moments and ion ization energies from 
iLidd l|200Ch ; lMcCIellarJ l|l963t ): iMillerl l|l990l 1. 

e Lacking values of fi for these species, the predictions of Eq. I|61|l 
are lower bounds on et- 



Figure [32] shows the experimental binding energies for 
alkanes, aromatic molecules, and methyl halides and cal- 
culated binding energies for several atoms, as a function 
of their ionization potential. The magnitudes of Eb calcu- 
lated for atoms and measured for molecules are quite sim- 
ilar (i.e., Eb < 0.5 eV). The binding energies for smaller 
alkanes, aro matics, and methyl h alides follow a "model 
atom" curve (|Mitrov et aUll999h . However, the ioniza- 
tion energies for alkanes with more than n = 7 carbons 
remain practically constant, while the positron binding 
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FIG. 32 Binding energies of positron-atom and positron- 
molecule complexes as a function of thei r ionization potential . 
Calculations: squa r es, va rious atoms (|Bromlev and Mitrovl . 
120101; iMitrov et all 120021 ); dashed curve, model alkali atom 
(jMitrov et aU 1 19991 ). Experiment: circles, alkanes with n car- 
bons; diamonds, the aromatic molecules, benzene and naph- 
thalene; triangles, methyl halides (see Sees. HVl IVl and I VI[) . 



The dipole polarizability ad characterizes the strength 
of the long-range positron-molecule attraction. Figure [33l 
shows positron binding energies for molecules as a func- 
tion of ad- Also shown a re the theoret i cal va lues for three 
atoms, Be, Zn, and Cd (|Mitrov et aLl . l2002l) . This figure 
suggests an approximately linear relationship between Sb 
and ad for gro ups of homologous spec ies. Based on this, 
the analysis of iDanielson et all (|2009[ ) began with a lin- 
ear fit of Eb to ad for alkanes. However, as seen in Fig. [331 
this fit generally underestimates the binding energies for 
other classes of molecules, most notably those with per- 
manent dipole moments and aromatic molecules. This 
motivated the inclusion of two additional parameters in 
the analysis, namely the molecular dipole moment \x and 
the number of 7r bonds for aromatic molecules . Such 
use of the bonds is similar to the approach used by iMillerl 
(1990) to parametrize the molecular polarizability. 

The resulting fitting function, using the numerical val- 
ues of the molecular pa rameters as listed in Table [VTl is 
([Danielson et aUl2009[ ). 



e b = 12.4(a + 1.6^ + 2.4^ -5.6) [meV], (61) 

where £& is in meV, ad is in units of A 3 , and /i is 
in units of debye. This expression can be viewed as 
a lowest-order Taylor expansion of Eb as a function of 
the variables ad, fJ-, and N w . This fit to the bind- 
ing energy data is shown in Fig. [34j Generally, the 
agreement between the predictions of Eq. (fSTj) and the 
measurements is quite good. The most significant dis- 
crepancies are 1-chlorohexane, acetylene (C2H2), ethy- 
lene (C2H4), and 2,2-difluoropropane. For all but 2,2- 
difluoropropane, Eq. (I61[) underestimates the binding en- 
ergy. These molecules probably possess some additional 
attraction beyond that described by Eq. (|6"Tj) . With the 



38 



l-Cl-Hexane 




1-F-Nonane 
M-F-Hexane 
1-F-Bulane 



nn-24 
cm J 



35 



FIG. 33 Measured positron binding energies Sb as a function 
of the dipole polarizability q^: large circles, alkane molecules 
used in the linear fit shown by the solid line; smaller circles, 
alkane-related molecules; triangles, molecules with permanent 
dipole moments; squares, aromatics with ir bonds. Dotted 
lines, guides to show the linearity of et for the different se- 
ries; diamonds, calc ulated et for the atoms Be, Zn and Cd 
(|Mitrov et al.l 12002), shown for comparison. 



exception of 1-chlorohexane and 2,2-dinuoropropane, the 
outliers all have double and triple bonds. The addition of 
the N„ term in Eq. (|6ip . which improves agreement for 
the aromatic species, overestimates £{, for these smaller 
molecules. Thus, while a similar effect may well be op- 
erating in small molecules with ir bonds, its magnitude 
appears to be smaller than in the aromatics. 




Napthat 




l-Cl-Hexane 

Benzene 



1-F -Propane 
2,2-F-Propane 



-50 50 100 150 200 250 300 350 
12.4 (a + 1.6u + 2.4N, - 5.6) [meV] 



FIG. 34 Binding energy from the fit Eq. (|61|l (solid line) using 
the polarizability ad, dipole moment /i, and the number of tt 
bonds Nt, for aromatic molecules. Symbols as in Fig. 1331 

According to Eq. (|oTj) . binding is assured if either 
ad > 5.4 A 3 or /i > 3.6 D. The first of these conditions 
is similar to that for positron bin ding to a mo d el "aJ - 
kali atom", namely a d > 3.5 A 3 (|Mitrov et all Il999l) . 
The second condition can be compared with the the- 



oretical critical value fi c = 1.625 D, beyond which 
the dipole supp orts an infinite series of bound states 
( Crawfor dL Il967h . Probably more relevant is a "practi- 
cal" value /i c « 2.5 D required to obtain binding energies 
> 1 me V in negative ions (i.e., electron-molecu le bound 
states) (jAbdoul-Carime and DesfrancoisL Il998l ): see be- 
low. These comparisons show that the threshold values 
of ad and fi from Eq. (|6ip appear to be quite reasonable. 

Another feature of Eq. (|6"Tj) , which is supported by ob- 
servations (Table |VT|). is that it predicts that the binding 
energies for fully deuterated hydrocarbons are close to 
their hydrogenated analogs (cf. deuterated butane, ben- 
zene, and nonane) and that isomers have similar values 
of Sb (e.g., pentane vis-a-vis isopentane). 

For the species studied to date, the largest binding 
energies for a given number of atoms are seen in the aro- 
matic molecules. A study of larger species (i.e., polycyclic 
aromatic hydrocarbons), such as the three- and four- ring 
variants anthracene and pyrene, will be of interest. If ej, 
increases beyond the energy of the C-H stretch mode (i.e., 
the highest-energy mode), the resonant energy thresh- 
old in these molecules becomes negative and the VFR 
of the vibrational fundamentals become inaccessible to 
the positron. This has already been observe d in deuter- 
ated naphthalene ( Young and Surkol l2008bl ) . Neverthe- 
less, the annihilation spectra for these molecules may 
still exhibit VFRs associated with positronically excited 
states, or those where the positron excites higher-lying 
overtone and combination vibrations. 

For molecules with similar chemical composition, a d 
increases with molecular size. However, the linear in- 
crease in Sb with ad predicted by Eq. (I61[) is expected 
to saturate at some point (e.g., when the positron de 
Broglie wavelength becomes smaller than the size of the 
molecule). The data in Figs. |3"31 and IMl indicate that the 
experiments have not yet reached this limit. 

The predictions of Eq. (|61l) can also be compared with 
available calculations for positron binding to molecules. 
All of the molecules listed in Table IVIII are predicted to 
bind positrons in agreement with the theoretical calcu- 
lations. The absolute values of Eb predicted by Eq. (joTj) 
for HCN and formaldehyde (H2CO) agree to within a 
factor of 2 with the calculations. For urea and acetone 
Eq. (|6"Tj) predicts much larger values than the calculated 
£h. This confirms the expectation ( Strasburgerl 2004 ) 
that t he binding energies calculated by iTachikawa et all 
(2003) :br urea and acetone are underestimated. 

For LiH, the calculated Sb is ten times greater than 
the prediction of Eq. (|6 1 [) . This is also true for other al- 
kali hydrides, Na H, KH and RbH, where the calculations 



give Eh > 1 eV ( Buenker et all l2005t iGianturco et all 



2006). This discrepancy is due to the fact that Eq. (|6"T|) 
is based on binding energies for molecules with rela- 
tively large ionization energies. The physical picture 
of such bo und states ga i ned fr om positron calculations 
for atoms ( Mitrov et all 120021 ). describes them as the 
positron moving in the attractive potential of the neu- 
tral molecule. In contrast, the ionization energy of LiH 
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TABLE VII Comparisons of theoretical predictions for 
positron- molecule binding energies with those of Eq. (I61|l . 



Name 


Formula 


ad 


1" 


st a 


s b b 






(A 3 ) 


(D) 


(meV) 


(meV) 


Hydr. cyanide 


HCN 


2.5 


3.0 


21 


35 


Formaldehyde 


H 2 CO 


2.8 


2.3 


12 


19 


Urea 


(NH 2 ) 2 CO 


9.7 


4.6 


140 


13 


Acetone 


(CH 3 ) 2 CO 


6.4 


2.9 


66 


4 


Lith. hydride 


LiH 


3.8 


5.9 


94 


1000 



"Predicted from Eq. 16U . 
''Values from quantum -ch emist ry 
llChoinacki and S trasburgerl l2006t), 



calculationsj HCN 

H9C O llStrasburgerl , 



l2004h. urea and acet o ne ([Tachikawa et all 120031). and LiH 
llBubin and Adamowic3. |2004| ; iMella et all feoOOl ; IStrasburgerl , 
120011) . 



and the other alkali hydrides is close to the binding en- 
ergy of the Ps atom (e.g., Ei — Ep s = 1.1 eV for LiH). In 
this case, the calculations indicate that the relevant phys- 
ical picture is closer to that of a P sH complex attached 
to th e positively charged metal ion (jBuenker et al. , 2005. 
2006) ; cf. Fig. [5] Thus it is not surprising that there is a 
discrepancy between the predictions of the Eq. (|oTj) and 
calculations of Eb for LiH. 

Equation (j6l"j) can be used to identify candidate 
molecules for further theoretical and experimental stud- 
ies. Theoretical calculations favor molecules with small 
numbers of atoms and simple electronic structures. Ex- 
perimental studies require species with vapor number 
densities > 10 -7 amagat at moderate temperatures (e.g., 
T < 500 °C) that are stable at this temperature and not 
deleterious to the vacuum system. Most convenient for 
study are small molecules with relatively large binding 
energies (e. g., Eb > 20 meV) . 

Recently iDanielson et al ] (|2010h measured £& for sev- 
eral molecules of this type. For carbon disulfide CS2 
(ay = 8.8 A 3 and fi — 0), the measured value Eb — 
75 meV is a factor of 2 greater than 40 meV, predicted 
by Eq. (frJTj) neglecting w bonds. Other molecules stud- 
ied, which were chosen for their simplicity and relatively 
large dipole moments, include acetaldehyde CH 3 CHO 
(ad = 4.6 A 3 , fi — 2.75 D, and Eb = 90 meV), acetone 
(CH 3 ) 2 CO (a d = 6.4 A 3 , n = 2.9 D, and e b = 173 meV), 
and acetonitrile CH 3 CN (a d = 44 A 3 , = 3.9 D, and 
Eb = 180 meV). All three have binding energies signifi- 
cantly larger than those predicted by Eq. (|61[) . In this 
previously unexplored regime in which /i>2D, the de- 
pendence on both yit and ad is much stronger than that 
given by Eq. (|61[) . These data provide a significant new 
opportunity to make quantitative comparisons between 
theory and experiment. It is hoped that they will stim- 
ulate new theoretical calcul ations. Thus a re c ent C I cal- 
culation for acetonitrile by iTachikawa et al. I (|2010h pre- 
dicted Eb = 135 meV, which is within 30% of the experi- 
mental value. 



Equation (ftJTj) could also be useful in describing the be- 
havior of low-energy positrons in a wide range of chemical 
environments, including biological systems (e.g., in PET 
analysis) and in materials science in conjunction with 
studies of insulators using techniques such as ACAR and 
positron annihilation lifetime spectroscopy. 



B. Comparison with negative ions of molecules and clusters 

Unlike positron-molecule complexes, negative molec- 
ular ions have been studied extensively. They can be 
separated into two classes. There is a class of stable 
molecular anions with large el ectron affinities ~ 1-3 eV 
( Rien stra-Kiracofe et al ., 2002). In these anions, the ex- 
cess electron occupies an unfilled orbital in the valence 
shell of the molecule, and they have no positron analog. 

A second arguably less well-studied class of negative 
ions is closely related to the positron-molecule com- 



1998 



1994 



plexe s discussed here (lAbdoul-Carime an d Desfrancois 



2002t Desfrancois et al 



.Abdoul-Carime et al 
l2004h ." In these "dipole-bound anions," the elec- 
tron is only weakly bound to the target by a combi- 
nation of dipolar, quadrapolar, and polarization forces 
(lAbdoul-Carime and Desfrancoisl Il998t iGutowski et all 
19981 ). Because of the small binding energy, the wave 
function of the excess electron is diffuse, residing predom- 
inantly outside the molecule due to the Pauli exclusion 
principle. In the positron-molecule case, the analogous 
short-range repulsion is due to the positive charges of 
atomic cores. The minimum dipole moment required for 
the formation of stable anions of common closed-shell 
molec ules has been determine d experimentally to be 2- 
2.5 D ( Desfrancois et al] . \l994) . in qualitative agreement 
with the critical dipole moment of 3.6 D, predicted by 

Eq. (J6T1). 

Table I Villi lists the measured and calculated 
binding energies and relevant molecular parame- 
ters for a selection of these dipole -bound anions 
( Abdoul-Carime and Desfrancoisl Il998() . This table 
also lists the positron binding energies from quantum- 
chemistry calculations (where available) and estimated 
from Eq. (|6ip . Positron binding energies are gener- 
ally greater than their electron counterparts. This 
could be related to an electron-positron correlation effect 
known as virtual Ps formation, which contributes sig- 



nificantly to the positron-atom at traction (|Dzuba et all . 
119951 : iGribakin and Ludlowl . I2004D . One consequence of 
this is that positrons are predicted to bind to atoms such 
as M g, Zn, and Cd, wh i ch do not form sta ble negative 
ions (|Dzuba et all Il995fc iMitrov et al. 1. 120021 ). 

Another interesting electron analog of positron- 
molecule bound states is the case of ne gative ions of 
small molecular clus ters, such as (N^O)™ ( Weber et all 
Il999h and (C0 2 )„ (|Leber et all l2000h. Using laser 
assist ed photoelectron attachment 
Il994f ) have observed prominent VFRs 
fragment negative ions [i.e., (N 2 0) g 



Usi ng 
Klar et al. I (|!992al lbl. 

in the yields of 
0~ with q < n, and 
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TABLE VIII Measured and calculated electron binding ener- 
gies and positron binding energies (in meV). 



recent, energy-resolved d ata for Z c r described above 
(Bar nes et all 2003, 2006) and the current understand- 



Molecule 


Formula 


M° 




Electron 


Positron 






(D) 


(A 3 ) 


e b b 








Formaldehyde 


CH 2 


2.33 


2.8 




0.02-0.05 


18 


12 


Acctaldchydc 


CH 3 CHO 


2.75 


4.6 


0.65 


0.95 




42 


Acetone 


(CH 3 ) 2 CO 


2.88 


6.4 


2.8 


1.6 


4 


66 


Hydr. cyanide 


HCN 


2.98 


2.6 




3.3-5 


35 


21 


Nitrobenzene 


C 6 H 5 N0 2 


4.2 


13.5 


28 


30.4 




270 



"Values used by Abdoul-Carimc and Dcsfrancois (1998). 
'Experimental data llAbdcml-Carime and De sfrancoist, Il998h . 

c Electrostatic m odel calculations 

l|Abdoul-Carime and Dcsfrancois. 1998|)- 
^Calculated values, see Table FVlTI 
e Predictions of Eq. (16 1 D - 



(C02)g with q < n]. These experiments are particu- 
larly relevant to the resonant processes discussed in this 
review. First, these VFRs were clearly identified with in- 
dividual molecular vibrational modes. Second, these res- 
onances displayed downshifts with increasing cluster size 
(especially clear for the CO2 clusters), which provided a 
measure of the electron-cluster binding energies. Finally, 
these cluster anions have the structure of a weakly bound, 
diffuse excess electron attached to an essentially unper- 
turbed neutral cluster, similar to the positron-molecule 
bound states described here. 

For the CO2 clusters, the measured electron binding 
energies are comparable to the positron eb for alkanes 
with similar numbers of carbon atoms. They also in- 
crease approximately linearly for n = 4-20, with some 
evidence of saturation at larger n. This behavior has 
been successfully modeled theoretically by combining the 
— ad/2r 4 polarization outside the cluster w i th a c onstant 
short-range potential inside ( Leber et all l2000h . It is 
likely that including the long-range — cvd/2r 4 potential 
would also improve the modeling of positron binding to 
alkanes (cf. Fig.©. 

A theory has been constructed to describe electron 
collisions with van der Waals clusters, such as (C02) n 
(jFabrikantl . 120051: iFabrikant and Hotopl . 120051) . It allows 
one to calculate the attachment cross sections and de- 
scribe the VFR that are observed. In this theory the elec- 
tron interaction with vibrational degrees of freedom of 
the CO2 monomers is described in the dipole approxima- 
tion. This is analogous to the approach used in Sec. III.FI 
to describe positron-molecule VFR. 



VII. ANALYSIS OF ANNIHILATION RATES MEASURED 
WITH THERMALIZED POSITRONS 

There is an extensive body of experimental data on 
the annihilatio n of thermalized positrons in mo l ecular 
gases at 300 K dHevland et aUll982t llwata et all 1 19951: 



ing of resonant annihilation. In particular, it is now pos- 
sible to consider the relationship between the thermal 
Z c s values and those due to VFR and IVR. Figure [24] 
shows a comparison of energy-resolved Z e g spectra for 
alkanes with n = 3-8 carbons. For each molecule, its 
Z c g was normalized by g = -^/et/e, shifted upward in 
energy by the binding energy and then normalized to 
unity at the C-H stretch peak. In this representation, it is 
clear that the spectral shapes of the alkanes change rela- 
tively little with molecular size. As the size of an alkanc 
increases, its Z e g spectrum shifts to lower energies. A 
mode accessed at 50 meV incident positron energy in the 
butane (n — 4) spectrum will be accessed at 5 meV in 
hexane [n = 6). Thus, the annihilation rates Z^g mea- 
sured with thermalized positrons at 300 K are intimately 
related to the corresponding values of the energy-resolved 
spectra, albeit shifted downward by the binding energy. 

As shown in Fig. [24] the 300 K data (plus signs in cir- 
cles) align well with the energy-resolved data when they 
are normalized by g = y^e^/sr (and the C-H stretch 
peak heights) and assigned a corrected energy of Sb + £t, 
where £t = iksT = 37.5 meV is the average thermal 
positron energy at 300 K. Also note that, due to their 
larger binding energies, alkanes with nine or more car- 
bons begin to sample the trough in the spectrum between 
the plateau and the C-H stretch peaks. This explain s 
the trend, pointed out previously ( Barnes et all 120031) . 
that the ratio of Z^g to Z e s at the C-H stretch peak de- 
creases by a factor of 2 when the number of carbon atoms 
in the alkane is increased from n < 9 to n > 9. This 
analysis connects in a quantitative way the thermal data 
with the energy-resolved measurements taken at higher 
positron energies. Thermal positrons at 300 K annihi- 
late by the same IVR-enhanced resonant mechanism as 
higher-energy positrons; the only difference is the specific 
vibrational resonances involved. 

The values of measured at 300 K can also be used 
to test Eq. (|6Tj) which estimates the pos itron binding en- 
ergy in terms of molecular parameters ( Danielson et all 
2009). A selection of molecules and their thermal anni- 



hilation rates Z^ from I Danielson et al. I (|2009D are listed 



in Table IIXI While these Z^g values do not provide di- 
rect evidence of VFR and hence positron binding, they 
can indicate whether positrons do or do not bind to 
the molecular species ( Gribakin. I2000ll200ll:flwata et all 



iPaul and Saint-Pierrel Il963t ISurko et all ll988bT . It is 
useful to examine these results in light of the more 



l200Ct lYoung and Surkd . l2008bh . We take Z% > 10* 
(Sec. III.D[) to be an indicator of resonant annihilation 
and hence positron binding. This is by no means a nec- 
essary condition, and so in the analysis below, we also 
use the somewhat arbitrary condition, Z* s /Z = 10, as 
the boundary separating the two groups of molecules. 
Namely, Z^/Z > 10 for a given molecule is taken to 
mean that positrons bind to this target and vice versa. 
The predictions of Eq. (frJTj) are listed in Table 11X1 to- 
gether with the molecular parameters used to calculate 
them. Generally, molecules with relatively large values 
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of ZfgjZ are predicted correctly by Eq. (|6Tj) to have 
positive binding energies. One borderline exception is 
methane which has a value of Z^g/Z — 14 but does not 
bind positrons. 

TABLE IX Values of Z%/Z for a variety of chemical species, 
the predictions of Eq. ()61[) f° r their binding energies, the 
molecular parameters relevant to this analysis, and the molec- 
ular ionization energies Ei. 



Name 


Formula 


(meV) 


Z 


Z%/Z b 


(A 3 ) 


(D) 


(eV) 




Small molecules 










Carbon dioxide 


co 2 


-36 


22 


2.5 


2.7 





13.8 


Sulfur hexafl. 


SF 6 


-14 


70 


1.2 


4.5 





15.3 


Water 


H 2 


-15 


10 


32 


1.5 


1.9 


12.6 


Nitrous oxide 


N 2 


-29 


22 


3.5 


3.0 


0.2 


12.9 


Nitrogen dioxide 


N0 2 


-26 


23 


47 


3.0 


0.3 


9.8 


Ammonia 


NH 3 


-12 


10 


160 


2.3 


1.5 


10.2 


Methane 


CH 4 


-37 


10 


14 


2.6 





12.7 




Alkcncs and alkyncs 








1-hcxcnc 


CeHi 2 


81 


48 


3900 


11.6 


0.3 


9.5 


trans 3-hcxcnc 


CeHi 2 


74 


48 


4100 


11.6 





8.9 


1,3-hcxadicnc 


CeHio 


72 d 


46 


8500 


11.4 


_ 


8.5 


1,3,5-hexatricnc 


CeHg 


69 


44 


9400 


11.2 





8.3 




Pcrhalogcnatcd 


alkancs 








Carbon tetrafl. 


CF 4 


-34 


42 


1.2 


2.9 





16.2 


Hcxafiuorocthanc 


C 2 F 6 


-10 


66 


2.3 


4.8 





14.6 


Pcrfiuoropropanc 


C 3 F 8 


13 d 


90 


1.7 


6.7 




_ 


Pcrfiuorohcxanc 


C 6 F 14 


84 


102 


3.3 


12.4 





12.8 


Pcrfiuoro octane 


C 8 F 18 


131 


210 


5.1 


16.2 





12.6 


Carbon tctrachl. 


CC1 4 


58 


74 


130 


10.3 





11.3 


Carbon tctrabrom. 


CBi-4 


120 


146 


270 


15.3 





10.3 


Carbon tctraiod. 


CI 4 


235 


218 


37 


24.5 









Partially fiuorinatcd alkancs 








Din. methane 


CH 2 F 2 


0.4 


26 


31 


2.5 


1.8 


12.6 


Trifiuoromcthanc 


CHF 3 


-4 


34 


7.3 


2.7 


1.7 


14.8 


Fluorocthanc 


C 2 H 5 F 


21 


26 


120 


4.2 


2.0 


12.4 


1,1, 1-trifi. ethane 


C 2 HsF3 


29 


42 


38 


4.2 


2.3 


13.3 




Oxygen-containing 


molecules 








1-propanol 


C 3 H s O 


50 


34 


590 


7.0 


1.7 


10.2 


Acetone 


C 3 H 6 


06 


32 


3100 


6.3 


2.9 


9.7 


Other 


aromatics and substitute 


d benzenes 






Anthracene 


C14H10 


422 


94 


46000 


22.8 





7.5 


Dccahydronaphth. 


CioHis 


151 


78 


5000 


17.7 





9.4 


o-xylcnc 


CsHio 


208 


58 


3100 


14.1 


0.6 


8.6 


Toluene 


C7H8 


179 


50 


3800 


12.3 


0.4 


8.8 


Hcxafiuorobcnzcnc 


C 6 F 6 


141 


90 


13 


9.8 





9.9 


Octafluorotolucnc 


C 7 F 8 


lG5 d 


114 


11 


11.7 




9.9 


Octafluoronaphth. 


CioF 8 


272 


132 


23 


15.5 





8.9 


Nitrobenzene 


CeHsNOs 


254 


64 


6700 


12.1 


4.2 


9.9 


Chlorobcnzcnc 


C6H5CI 


204 


58 


1250 


12.1 


1.6 


9.1 


Bromobcnzcnc 


C 6 H 5 Br 


215 


70 


2300 


13 


1.7 


9.0 


Fluorobcnzcnc 


C 6 H 5 F 


176 


50 


900 


10.0 


1.6 


9.2 


1 ,2-difluorobcnzcnc 


C 6 H 4 F 2 


189 


58 


570 


9.8 


2.4 


9.3 



"Predicted by Eq. J6T1 . 

b Z*% i s measured at room temp erature with thermalized positrons 
l llwatal . 1 19971 ; llwata et a/1. 11995T) . 
^Parameters from lLidd feuutfl ; iMcClellanl [119631 ); iMilleil dl99Ct) . 

^Lacking values of fj, for these species, the predictions of Eq. I|61|l 
are lower bounds on ej. 

As shown in Table IIX1 the criterion from Eq. (|ST|) that 
molecules will bind for ay > 5.4 A 3 has the consequence 



that most large molecules can bind positrons whether 
they have a permanent dipole moment or not. The alka- 
nes are an example of this, where only methane does 
not bind. For nonaromatic molecules the ay term dom- 
inates for all but relatively small molecules with corre- 
spondingly small values of a ( i- A related trend is seen 
in halogen substitution, where Eb rises rapidly as the size 
(and hence the polarizability) of the halogen is increased. 
For example, CF4 does not bind positrons, while CBr4 is 
predicted to have a binding energy in excess of 120 meV. 

Restricting comparison to cases where the model pre- 
dicts \Eb\ > 15 meV in deference to the likely error bars 
for the model, negative binding energies are predicted 
for almost all small molecules. In particular, the di- 
atomic molecules, H2, D2, N2, O2, CO, and NO (not 
shown in Table HX|) . with Z%jZ values ranging from 
2.3 to 7.4 (i.e., < 10) all have negative binding ener- 
gies £& < —40 meV ([Danielson et aLl . l2009t ). As shown in 
Table ITXl neither methane nor carbon tetrafluoride bind, 
i n agreement with the interpretation of their Z^g values 
( Iwata et aLl . Il995h . As indicated in Table IVll molecules 
for which the binding energies are "too close to call," 
[i.e., |eb| < 15 meV, based on Eq. (fBTj) ] include ethane, 
propane, cyclopropane, methanol, methyl fluoride, and 
water (cf. Table lIXp. and this is con sistent with the ob- 
servations (|Young and Sur k<J l2008ch . There are some 
disagreements in cases where the binding energy is pre- 
dicted and/or observed to be reasonably small. From 
Table ITXl the only molecule that exceeds the |£b| > 15 
meV criterion is NO2 with a Z^g jZ ratio of 47 and yet a 
predicted binding energy of —26 meV. Like benzene, NO2 
has two resonant electronic states involving 7r bonds, but 
the geometry of its valence orbitals differs greatly from 
that of the aromatics. Additional experiments are needed 
to better distinguish the effects of the different types of 
electronic bonds, especially for the smallest molecules. 

As shown in Table IIX1 the ratios of Z$f jZ for perfluo- 
roalkane molecules are much smaller than those for alka- 
nes. This has been interpreted as evidence that resonant 
annihilation is "switched off" f or them, possibly due to 
lack of binding ([Gribakinl . l2000h . However, Eq. (|61l) pre- 
dicts that perfluorocarbons other than perfluoromcthanc 
and perfluoroethane can bind positrons. This warrants 
further investigation. 



VIII. OTHER TOPICS 

A. Gamma-ray Doppler-broadening measurements 

As discussed in Sec. III. Bl positron annihilation usually 
results in the production of two, approximately back-to- 
back gamma rays, each with an energy E 1 ~ 511 keV. 
However, there are small shifts in the energies of the 
two gammas that can be used to obtain microscopic 
information about the annihilating pair. According to 
Eq. (fT3"l) , these energies are Doppler shifted by the energy 
e = ±cP z /2, where P z is the component of the center- 
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of-mass momentum of the electron-positron pair along 
the direction of the gamma rays. For the low-energy 
positrons that we deal with here (i.e., e < 0.5 eV), the 
center-of-mass momentum P is dominated by the mo- 
mentum distribution of the electron orbitals. As a bench- 
mark, a 4 eV electron traveling along the direction of the 
gamma rays produces a 1 keV Doppler shift. 

Extensive measurements have been made of the 
Doppler broadening of annihilation gamma rays for the 
case of a thermal distribution of positrons at 300 K in- 
teract i ng with a w i de var iety of molecules ( Iwata et all 
Il997at iTang et all Il992l ). The apparatus and proce- 
dures for these measurements are described in Sec. IIII.Fl 
Shown in Fig. [35] is the spectrum for hexane. The de- 
tector response is modeled by a Gaussian lineshape with 
a FWHM of 1.16 keV and an error function to account 
for the effect of Compton scattering of the gamma rays 
in the detector. The measured gamma spectra could be 
fit reasonably well using a single Gaussian to approxi- 
mate the Doppler spectrum, convolved with the detector 
response function. 

A more accurate fit was obtained using two Gaussians 
to model the Doppler spectrum. The second Gaussian 
was required to describe a smaller, higher-momentum 
component, which was typically about 1-3% of the larger 
component. The resulting linewidths (FWHM), with the 
detector resolution de-convolved, ranged from 1.7 keV for 
H2 to ~ 3.1 keV for fluorocarbons. The smaller, higher- 
momentum component had widths from 4 to 9 keV, with 
most molecules in the narrower range from 5 to 7 keV. 
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that the linewidths are only weakly correlated with the 
values of either £& or Z^. For example, comparing the 
ethane and dodecane in Table IXl £& increases from near 
zero to 220 meV, and Z*^ increases by 3 orders of magni- 
tude, while the gamma-ray line width increases by only 
5%. Note that there is, however, a significant increase 
in the linewidth when a hydrocarbon is partially or fully 
fluorinated. 



TABLE X Gamma-ray linewidths AE (FWHM, single Gaus- 
sian fit), annihilation rates for a thermal distribution of 
positrons at 300 K, and positron-molecule binding energies 
Eb, for selected molecules. 



Name 


Formula 


AE a 
(keV) 


rrth b 
^eff 


e b c 
(meV) 


Methane 


CH4 


2.09 


142 


< 


Ethane 


C2H6 


2.18 


1780 


> 


Propane 


C3H8 


2.21 


3500 


10 


Pentane 


C5H12 


2.24 


40200 


60 


Nonane 


C 7 Hi 6 


2.32 


643000 


145 


Dodecane 


C12H26 


2.29 


1780000 


220 


Benzene 


CeHg 


2.23 


15000 


150 


Naphthalene 


CioHs 


2.29 


494000 


~ 300 


1-fluoroethane 


C 2 H 5 F 


2.62 


3030 


> 


Hexafluoroethane 


C 2 F 6 


3.04 


149 


< 


1-fluorohexane 


C 6 H 13 F 


2.46 


269000 


80 



"From Iwata et al. (1997a). 



b See 
c See 



Iwatal il997h:llwata et all dl995ft and references therein. 
Young and Surko ( 2008b. c) and Sec. EH 
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FIG. 35 7-ray spectrum of hexane (C6H14): (o) observed 
spectrum; solid curve, single-Gaussian fit; dotted curve, fit us- 
ing a hydrogenic functional form; dot-dashed curve, fit using 
a noninteracting hydrogenic form convolved with a Gaussian; 
dashed curve (indistinguishable from the data), two-Gaussian 
fit. The statistical error bars are c omparable to or sm aller 
than the size of the data points. See llwata et all (|l997af ) for 
details. 

Shown in Table [X] are typical values of the linewidths 
AE (FWHM) using the single Gaussian fit, the annihi- 
lation rates Z e R for thermal positrons, and the positron- 
molecule binding energies £5, for a selection of molecules. 
A key conclusion of these Doppler-broadening studies is 



The linewidths AE for alkanes, shown in Table [Xj in- 
crease monotonically with increasing molecular size. This 
increase has been ascribed to the change in the relative 
number of electrons in C-C vs C-H bonds, since the elec- 
trons in these bonds are expected to be characterized by 
different Doppler widths. An analysis was made to de- 
termine the fraction of annihilation events i nvolving elec- 
trons f rom these types of orbitals in alkanes ( Iwata et all 
Il997al) . It showed that there is a linear increase in the 
linewidth as a function of the fraction of electrons in C-C 
orbitals. This is consistent with the assumption that the 
positron annihilates statistically on any of the valence 
electrons in the molecule (see below). 

A similar study was carried out for the linewidths asso- 
ciated with annihilation in partially and fully fluorinated 
alkanes. In this case, the observed line was decomposed 
into fluorinated and hydrogenated components by fitting 
to a sum of two measured line shapes, one for the fully hy- 
drogenated compound, a nd one for the fully fluorinated 
compound ( Iwata et all Il997al ). As shown in Fig. |3T>1 
there is a smooth, linear increase in the linewidth over the 
full range of the degree of fluorination, from to 100%. 
This also supports the hypothesis that the positrons an- 
nihilate statistically on any valence electron. 

The results of these gamma-ray studies provide con- 
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FIG. 36 Normalized fraction of positrons annihilating on flu- 
orine atoms in partially and fully fiuorinated alkanes, plot- 
ted against the fraction of valence electrons in the fluorine 
atoms. Open symbols are methane (circle), ethane (square), 
propane (triangle) , and hexane-based molecules (inverted tri- 
angle); filled symbols are six-carbon benzene-based molecules: 
1,2-difluorobenzene (square), 1,3-difluorobenzene (triangle), 
1,4-difluorobenzene (inverted triangle), and othe r six-carbon 
benzen e- based molecules (circle). See text and llwata et al\ 
l|l997al ) for details. 



firmation of a model by ICrawfordl ( 19941 ) for positron 
annihilation on molecules. It was developed to explain 
the degree of ion fragmentation that is observed follow- 
ing positron annihilation on molecules (see Sec. IVHI.BI) . 
Crawford argued that a low-energy positron should anni- 
hilate statistically on any valence electron. The gamma- 
ray linewidth measurements shown in Fig. are consis- 
tent with this hypothesis. Such uniform spreading of the 
positron density over the molecule is related to the fact 
that at low energies the positron de Broglie wavelength 
AdeB is larger than the size of the target. This prevents 
the positron from being localized strongly on any par- 
ticular site in the molecule. This is true both for direct 
annihilation, for which AdeB = 27r/fc (e.g., k — 0.05 a.u. 
for thermal positrons at 300 K), and for resonant anni- 
hilation, for which AdeB ~ 2tt/k, where k = v2~£& (see, 
e.g., Fig. [TU] for model bound-state wave functions). 

In some cases, such as strongly polar molecules, local- 
ization of the positron near specific sites in the molecule 
might be expected. For example, in LiH the positron den- 
sity is strongly localized at the negatively charged hy dro- 
genic end of the molecule ( Strasburgerl [l999, 2001); cf. 



Fig. [SJ However, annihilation gamma-ray spectra from 
this class of molecules have not yet been investigated. 

An analysis is in progress to understand in more de- 
tail th e implications of t he Doppler-broadening experi- 
ments (jWang et all I2009T) . It relies on modern quantum- 
chemistry methods, such as the density-functional theory 
B3LYP/TZVP, which give electron momentum de nsities 
in good agreement with experiment (Wang, 2003). The 
linewidths calculated for methane and fluoromethanes, 
ethane, propane, butane and benzene, by taking only 



the electron momentum distribution into account, are 
about 30% greater than the values measured. The main 
source of this discrepancy appears to be the neglect of the 
Coulo mb repulsion betwe en the positron and the atomic 
cores 1 Wang et all I2009T) . This repulsion suppresses the 
positron wave function at small positron-nuclear sepa- 
rations, which effectively reduces the high-Doppler-shift 
components in the annihilation spectra. 

A more complete theory of the annihilation gamma-ray 
spectra for molecules should include the full electron and 
positron dynamics and account for the positron-nuclear 
repulsion and elect ron-positron correl ation effects. How- 
ever, the results of lWang et al\ ( 20091) suggest that these 
effects in positron annihilation spectra can be modeled 
by a relatively simple scaling factor. At the next level of 
analysis, one should examine the high-momentum com- 
pon ent in the spectra [e .g., as described by two-Gaussian 
fits (jlwata et qj.l . ll997bh ]. which is likely due to annihila- 
tion on inner-shell electrons. Further work on this topic 
is in progress. 



B. Annihilation-induced fragmentation of molecules 

Two-body collisions between positrons and molecules 
can produce positive ions by two mechanisms. For in- 
cident positron energies greater than the threshold for 
Ps formation, ionization proceeds via this channel. If 
the positron energy is below the Ps formation threshold 
(i.e., the principal regime of interest in this review), the 
positron can annihilate with a molecular electron, also 
producing a positive molecular ion. 

The initial studies of annihilation-induced ion for- 
mation in molecules below the positronium formation 
threshold were conducted using positrons confined in a 
buffer-gas trap in the presence of low pressure gases of 
alkane molecules ( Passner et all Il989l) . In spite of com- 
plications due to the presen ce of molecular nitrogen in 
the trap ([Glish et all , [l994), these experiments estab- 
lished that sub-Ps-threshold ionization can produce sig- 
nificant amounts of fragmentation. Subsequently, exten- 
sive studies of positron-induced ionization of molecules, 
both below and above the positronium formation thresh- 
old, were conducted by L D. Hulett and collabor a- 
tors (iDonohue et all Il99di: iHulett et all Il993lll996allbl : 



iMoxom et all . 120001; IXu" et all I1994L Il993l I1995L Il997j) 
The experimental procedures are described in Sec. IIII.GI 
These studies produced several important results. One 
is that the extent of fragmentation depends in a non- 
trivial way on the energy of the incident positron. For 
example, fragmentation is a minimum at energies near 
the Ps-formation threshold, and increases toward smaller 
and greater positron energies. Figure [37] shows the mass 
spectrum of 1-dodecene (C12H24) at an incident energy 
of 1.0 eV (i.e., ~2 eV below the Ps formation threshold). 
It illustrat es the broad spe ctrum of product ions that is 
produced (|Xu et al. , 1994). Figure 135] shows the cross 
sections for producing these fragments as a function of 
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positron energy. 
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FIG. 37 Time-of-flight mass spectra of ion fragments from 
1-dodecene at 1.0 eV, wh ich is 2.1 eV belo w the Ps formation 
threshold. Adapted from lXu et all (|l994f) . 




POSITRON ENERGY (eV) 

FIG. 38 Cross sections for fragmentation of 1-dodecene 
(C12H24) by positron annihilation near and below the Ps- 
formation threshold at 3.1 eV. From lXu et~cd\ (|l994h . 

The fact that in many molecules the degree of fragmen- 
tation has a minimum close to the Ps-formation threshold 
may offer the possibility of using positrons to advantage 
in ion mass spectroscopy. Another interesting effect is 
that, below the Ps-formation threshold, double and triple 
bonds can stabilize the species with respect to fragmen- 
tation. For example, in the series decane (C10H22, single 
bonds only), 1-decene (C10H20, one double bond) and 
1,9-decadiene (Ci H 18 , two double bonds), the yield of 
i on fragments de creases with each additional double bond 
(IXu et ■ a*.l . ll995h . This same effect was also observed in 
other molecul es [e.g., tetravyn ilsilane vs tetraethylsilane 
(|Hulett et all Il993t IXu et all ll995|V[ 



Motivated by these experiments. ICrawfordl (1994) con- 
structed a simple and insightful model of the fragmen- 
tation process. He argued that annihilation will occur 
with comparable probability on electrons in any valence 
molecular orbital i, not just the highest occupied molec- 
ular orbital (HOMO). As a result, the molecular ion is 



typically left in an electronically excited state, with the 
excitation energy provided by the energy difference be- 
tween the HOMO and the orbital i. Figure [5^1 shows the 
probability that this energy e xceeds a certain value for 
propane, hexane, and decane ( Crawford! . Il994f ). Due to 
vibronic coupling, this energy flows from the electronic to 
the vibrational degrees of freedom. Given time and suf- 
ficient energy (e.g., a few electron volts), the molecular 
ion can break up in the process known as unimolecular 
dissociation. 

Fragmentation can occur when the incident positron 
annihilates, either during direct annihilation or following 
capture into a VFR. The resonant process likely domi- 
nates in larger polyatomics at low positron energies and 
is central to this review. In this case, both the incident 
positron energy and the binding energy contribute to the 
total excitation energy of the molecular ion. The gamma- 
ray annihilation studies described above provided vali- 
dation of Crawford's model, showing that positrons do 
indeed annihilate with approximately equal probabilities 
on all valence electrons. 




1 2 3 4 5 

EXCESS ENERGY IN PARENT ION (eV) 

FIG. 39 Probability P(E) that the energy deposited into a 
molecular ion due to positron annihilation with an electron 
in a valence orbital below the HOMO is greater than a given 
excess energy for p ropane, hexane, a nd decane. For details of 
the calculation, see ICrawfordl l| 19941 ). 



The experiments of Hulett and coworkers used a beam 
with an energy resolution ~0.5 eV FWHM. This raises 
a question as to the extent to which ionization could 
be controlled (e.g., to produce only parent ions) using 
the much higher resolution beams that are now available 
(e.g., 40 meV FWHM). Another possibility would be to 
exploit ionization via the second bound states such as 
those observed in larger alkane molecules (cf. Sec. IV. Bl) . 
In this case, the positron wave function has a nodal plane 
at the center of the linear carbon chain (see Fig. ITO)) . 
which would be expected to produce a decrease in anni- 
hilation near this location. An experiment to test this 
might also provide information as to whether the ex- 
cess electronic energy deposited in the molecule in the 
annihilation process could diffuse away from the anni- 
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hilation site and then break a bond. Investigations to 
date have only scratched the surface of this rich area 
of matter-antimatter chemistry. We still know relatively 
little about the chemical specificity of annihilation and 
annihilation-induced ion production. 



C. Nonlinear dependence of annihilation on molecular 
density 

Equations ((2]) and (p~0|) imply that the annihilation rate 
A is a linear function of the gas density n. This was 
shown to be correct at low densities, such as those used 
in the positron-trap experiments described in Sec. MI. CI 
There are, however, a number of effects that can make 
Z e g density dependent. For example, in dense gases 
and at lower temperatures, positrons can cause a local 
phase transition and become self-trapped in clusters that 
gives rise to a strongly nonlinear behavior of the an- 
nihilation rate ( Iakubov and Khrapakl [l982). Also, in 
low-density gases one typically observes Z e g 3> Z, while 
at solid or liquid densities the value of Z c g for almost 
any material are clos e to the number of valence electrons 
(|Pomeranchuklll949l ). 

For diatomic molecular gases such as H 2 , N 2 , CO 
and O2 at densities n < 50 amagat, the Z e g values re- 
main practically independent of n. In contrast, other 
species, such as CO2, CH4, and SF« displ ay consider- 
able variation with n ( Hevland et «L[l982|) . For larger 
molecules (e.g., C2H 6 , CCI2F2, CaHg, and C4H10) a 
stron g density dependence is observed even at n ~ 1 am- 
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agat ( Hevland et aZ.L [l982h ■ as shown in Fig. 1301 In these 
cases, the density effect can be expressed as a quadratic 
correction to Eq. (1111)) . 



A = an + bn 



(62) 



Here a — r Kr^cZ c g represents annihilation in binary col- 
lisions, and the b term accounts for three-body annihi- 
lation events involving a pos itron and two gas molecules 
(|Charlton et all I2002L ,12006). A key finding is that the 



coefficient b for molecules with large Z c g, for which res- 
onant annihilation is observed (e.g., C2H4 or CsHg), 
is approximately proportional to Z c g. This suggests 
that three-body collisions affect the resonant annihila- 
tion, rather than being an additional independent reac- 
tion pathway. 

There is currently no information about the micro- 
scopic nature of this process, and so one can only men- 
tion some possibilities. The annihilation in the resonant 
positron-molecule complex cou ld be enhanced through 
collisions with other molecules. iGribakinl ( 20091) recently 
proposed that these collisions may stabilize the resonant 
complexes. If this stabilization is complete (i.e., prevent- 
ing the positron from escaping the molecule), then the 
corresponding collision rate T c = a c v m n must be added 
to the annihilation rate r" in the numerator of Eq. (T3"T1) . 
Here v m is the mean molecular collision velocity and o~ c 
is the stabilizing collision cross section. The values of o~ c 
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FIG. 40 Density dependence of the room-temperature 
positron annihilation rates for CCl2F2(D), C2H6 (o), and 
C 3 H s (A). Solid lines are fits of the form Eq. (J62J). Adapted 
from I Charlton et al\ (|2006r h 



required to explain the observations can be estimated us- 
ing t he experimental values of a and b ( Charlton et all 
120061 ) . and estimating T* from Eqs. <j23ji and ([36]) using 



the binding energies inferred from the energy-resolved 
Z e g measurements. The values range from a c ~ 1 a.u. 
in ethane and ethylene to ~10 a.u. in butane. 

There is also the possibility that a collision with a 
second molecule enhances the annihilation rate by an 
amount pT c , where p < I is the probability that the colli- 
sionally stabilized positron-molecule complex survives to 
annihilation rather than autodetaches. In this case, to 
achieve the same effect, the cross sections a c would need 
to be enhanced by a factor p^ 1 above that for complete 
stabilization. 

The actual mechanism of collisional stabilization is also 
unclear. It could be due to the transfer of the bound 
positron from the vibrationally excited host molecule to 
a vibrationally colder molecule of the same species. This 
process would be facilitated by the fact that the positron 
binding energies in the two molecules would be similar, 
making it a resonant charge transfer. The vibrational 
de-excitation of the positron-molecule complex could also 
be due to intermolecular vibrational energy transfer in a 
collision with a colder molecule. However, in both cases, 
for the species studied, the thermal vibrational energies 
at room temperature are larger than the positron binding 
energies. Hence these processes might not completely 
stabilize the complex. 

Finally, if a VFR is excited when two molecules are 
in close proximity, the positron binding energy is likely 
to be greater than that for a binary positron-molecule 
interaction. This too might help to explain the observa- 
tions. At a minimum, further research is warranted to 
understand this collisional enhancement effect. 
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IX. SUMMARY AND A LOOK TO THE FUTURE 

This review focuses on the interaction of positrons 
with molecules in the range of energies below the thresh- 
olds for electronic excitation and Ps formation. Exper- 
imental studies of annihilation resolved as a function of 
positron energy show that positrons bind to many molec- 
ular species. This enables the formation of long-lived 
resonant states (i.e., vibrational Feshbach resonances) in 
two-body positron-molecule collisions, thus altering the 
dynamics in a fundamental way for the incident positron 
energies in the range of the molecular vibrations. 

These VFR are responsible for the large annihilation 
rates observed in many polyatomic molecules. They ex- 
ceed by orders of magnitude the contribution of direct, 
"in-flight" annihilation. This enhancement distinguishes 
annihilation from conventional scattering processes. For 
example, the resonant contribution would be hard (if at 
all possible) to observe in elastic scattering, where much 
larger contributions come from potential scattering. 

Theoretically, resonant annihilation can be described 
using a Breit-Wigner approach. The key quantities re- 
quired to make predictions are energies of the VFR and 
their decay rates (i.e., the annihilation rate, positron cap- 
ture and autodctachment rates, etc.). Such a calcula- 
tion can be done for small molecules (e.g., the methyl 
halides and similar molecules with IR-active modes) us- 
ing a theory which makes use of the small parameter in 
the problem, namely, the positron binding energy. This 
enables the use of long-range dipole coupling to evaluate 
the positron elastic capture rate r e and the (universal) 
y/el scaling to estimate the positron annihilation rate T a . 
The theory also makes use of the fact that the vibrational 
spectra are sufficiently simple and that the IR strengths 
of the modes are known. The only free parameter is the 
positron-molecule binding energy, which can be deter- 
mined experimentally. In the case of methyl halides, the 
theoretical predictions are in excellent agreement with 
the experimentally measured annihilation rates as a func- 
tion of positron energy. Application to the deuterated 
methyl halides provides a complete test of the theory with 
no adjustable parameters. The theory has also been suc- 
cessfully extended to include combination and overtone 
vibrations. 

A general result of the theory is that the resonances for 
all modes with coupling strengths T e 3> r a produce an- 
nihilation resonances of the same magnitude, modulated 
only by the factor g = \J Sb/s. There are no known cases 
where dipole coupling strengths greater than T a fail to 
produce resonances. On the other hand, the theory is 
incomplete for molecules such as ethylene and acetylene, 
where explaining the measured annihilation spectra re- 
quires the inclusion of VFR associated with nominally 
IR-inactive modes and overtones and combination vibra- 
tions. In these cases, the couplings are more difficult to 
evaluate. 

As the number of atoms N in the molecule is increased, 
the magnitudes of the annihilation resonances exceed (of- 



ten greatly) those explicable on the basis of individual 
single-mode resonances. The annihilation rates in these 
species scale as Z a g oc gN iA . While a quantitative ex- 
planation for this N 4 scaling is lacking, it is likely that 
the process of intramolecular vibrational energy redistri- 
bution is responsible for this enhancement. This IVR 
process couples the single-mode resonant doorway states 
to baths of dark states, namely, states that are not di- 
rectly coupled to the positron continuum. In this pic- 
ture, N likely reflects the number of vibrational degrees 
of freedom of the molecule. 

For attachment and annihilation in large molecules, 
there is a related key piece of the theory that is missing. 
If the positrons were able to populate the VFRs associ- 
ated with all possible vibrational excitations, Z e ff would 
be expected to increase much faster with molecular size 
than is observed. While vibrationally inelastic escape 
channels could moderate such growth, they appear to be 
generally inoperative. Hence, it is not known presently 
what subset of multimode dark states and subsequent 
positron escape channels are accessible. By inference 
from the experimental results, detachment from these 
dark states likely takes place via a quasielastic escape 
channel or channels, including the original and nearby 
doorway states (e.g., dipole-allowed modes close in en- 
ergy to the doorway state through which the positron 
entered). 

The energy-resolved annihilation data have provided 
measurements of positron-molecule binding energies Eb 
for 30 molecular species. The molecular dipole polariz- 
ability plays an important role in fixing the magnitude of 
Eb- Binding is further enhanced in aromatic molecules by 
an amount that increases with the number of electronic tt 
bonds. In small molecules, the permanent dipole moment 
also increases Eb- These results indicate that most large 
molecules will bind positrons, and they provide some in- 
sight into promising candidates for further experimen- 
tal and theoretical study. The goal of finding molecular 
species for which the binding energy can be both cal- 
culated and measured is closer to realization with the 
recent discovery o f small molecule s with relatively large 
binding energies (|Danielson et all l2010t> . Other inter- 
esting topics for further research include study of very 
large molecules, such as polycyclic ar omatic hydrocar- 
bons that are of astroph ysical interest ( Guessoum et all 
2010; llwata et aLl.ll996f). and cagelik e structures, such as 
Ceo (|Gianturco and Lucchesd . [1999). 

Following positron annihilation, VFR-mediated or not, 
the remaining molecular ion is frequently found to frag- 
ment. There are many open questions in this area, in- 
cluding what determines the degree of fragmentation. A 
practical question of interest is whether positron-induced 
annihilation might be a way to produce unfragmented 
ions from large molecular species for applications such as 
mass spectrometry. 

Finally, the resonant enhancement of positron an- 
nihilation described here has an electron analog. In 
electron-molecule scattering, resonances are known 
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to drive processes such as dissociative attachment 
( Christophorou et aLl . ll984T) . a process that involves the 
slow motion of the heavy atomic nuclei and that would 
be weak in a direct electron scattering process. In many 
electron-driven reactions, shape resonances are promi- 
nent. While this is not the case for positrons, both 
positrons and electrons can populate VFR directly. In 
the positron case, even a weak positron- vibrational cou- 
pling is sufficient to "turn on" the resonant annihilation 
mechanism. It is possible that a better understanding 
of positron attachment to polyatomic molecules, gained 
through such annihilation studies, can be useful in un- 
derstanding similar electron-molecule processes. 
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